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INTRODUCTION 


Ever  since  the  first  man-made  Earth  satellite  was 
launched  in  1957,  the  scientific  community  has  undertaken 
extensive  and  intensive  studies  of  the  atmospheric  density 
at  satellite  altitudes.  The  critical  impact  of  the  neutral 
density  upon  orbital  motion  motivated  research  at  various 
institutions  and  resulted  in  diverse  models,  developed  for  a 
wide  range  of  space  systems  and  mission  applications. 

Atmospheric  neutral  density  research  and  neutral  density 
model  development  continue  at  numerous  government  and 
private  facilities.  Each  is  attempting  to  resolve  its  own 
problems,  but  the  issues  to  be  addressed  throughout  the 
community  remain  closely  related. 

Many  analytical  and  empirical  atmospheric  neutral 
density  models  have  been  developed  over  the  last  25  years, 
and  implemented  in  various  operational,  research,  and 
development  systems.  These  models  are  qualitatively 
equivalent  but  numerically  different.  In  applications  for 
example,  the  associated  ballistic  coefficients  for  a  near- 
Earth  satellite  can  be  significantly  different  if  the 
density  models  are  different.  Thus  the  orbital  data  can  not 
be  interchanged  without  experiencing  loss  of  accuracy.  This 
is  a  typical  example  of  astrodynamic  compatibility  problems. 
This  proliferation  of  incompatible  models  has  created 
serious  communication  and  data  sharing  difficulties; 
sophisticated  conversion  algorithms  may  be  required  to  bring 
about  compatibility  between  the  space  centers  and  their  data 
users . 

In  recent  years,  the  Air  Force  Space  Command  Directorate 
of  Operations  Analysis  (formerly  the  Directorate  of 
Astrodynamics  )  and  a  number  of  other  professional 
organizations  have  conducted  independent  studies  regarding 
the  state-of-the-art  atmospheric  neutral  density  models  and 
their  respective  performance  in  space  systems.  These 
studies  have  all  concluded  that,  despite  extensive  researcn 
efforts,  the  scientific  community  has  produced  no 
significant  improvement  in  either  neutral  density  model 
accuracy  or  computational  efficiency  in  the  last  two 
decades.  The  general  agreement  is  that  all  available 
empirical  neutral  density  models,  those  that  input  the 
measured  10.7  cm  solar  flux  and  a  geomagnetic  index,  still 
exhibit  an  average  error  in  the  neighborhood  of  15  percent 
and  that  the  newer  models  tend  to  be  more  complex  and  thus 
are  less  efficient.  In  addition,  the  solar  flux  and 
qeomaqnetic  indices  are  extremely  difficult  to  predict. 
Therefore,  the  prediction  of  satellite  orbits  and  positions 
remain  inaccurate  and  unreliable. 


\ri 


* 


< 


Many  government  agencies,  research  institutions,  and 
laboratories  are  expending  tremendous  effort  to  address  the 
atmospheric  neutral  density  development  and  are  making 
encouraging  progress.  However,  the  problem  is  extremely 
complex,  the  focus  of  interest  remains  diverse  throughout 
the  scientific  community,  and  research  and  devlopment  funds 
are  becoming  ever  more  limited.  Therefore  Air  Force  Space 
Command  decided  to  sponsor  this  conference  to  offer  a  forum 
to  coordinate  and  promote  atmospheric  neutral  density 
research  and  development  that  is  ultimately  targeted  at 
supporting  space  operations  requirements. 

Operational  neutral  density  model  users  were  invited  to 
discuss  their  operational  requirements  and  the  difficulties 
that  they  have  encountered,  and  atmospheric  scientists  were 
asked  to  present  recent  advances  and  prospects  for 
improvement  in  the  neutral  density  models.  Cooperation  and 
dialogue  among  members  of  the  atmospheric  neutral  density 
operational  and  research  communities  is  mutally  beneficial 
and,  in  fact,  synergistic.  The  scientific  community  needs 
relevance  to  justify  continued  research  during  this  era  of 
tight  fiscal  restraint,  and  they  need  data  to  develop  and 
validate  the  advanced  algorithms  that  will  ultimately 
comprise  improved  models.  The  operational  community  has 
cited  and  validated  requirements  to  specify  and  to  predict 
the  atmospheric  neutral  density  at  satellite  altitudes  with 
specific  accuracy  requirements  as  a  function  of  altitude. 
Progress  will  continue  only  so  long  as  the  various  agencies 
continue  to  provide  mutual  support. 

The  operators  were  invited  from  4th  Weather  Wing  ( 4WW ) , 
Air  Force  Space  Command  ( AFSPACECOM ) ,  Air  Weather  Service 
( AWS ) ,  Air  Force  Global  Weather  Central  ( AFGWC ) ,  Naval  Space 
Surveillance  Center  (NAVSPASUR),  and  United  States  Space 
Command  ( USSPACECOM  ) ;  atmospheric  scientists  were  invited 
from  Air  Force  Geophysics  Laboratory  (AFGL),  NASA  Goddard 
Space  Flight  Center  (GSFC),  NASA  Marshall  Space  Flight 
Center  (MSFC),  National  Oceanic  and  Atmospheric 
Administration  ( NOAA ) ,  National  Center  for  Atmospheric 
Research  ( NCAR ) ,  MIT  Lincoln  Laboratory  ( LL ) ,  and  the 
University  of  Michigan  ( UM  ) . 

Topics  for  the  presentations  were  selected  to  provide 
information  and  to  generate  discussion  among  the  conference 
participants.  The  specific  session  topics  included: 

o  Examine  the  status  of  the  current  operational 
requirements  and  operational  difficulties, 

o  Assess  the  state-of-the-art  atmospheric  neutral 
density  models  and  operational  feasibilities, 
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o  Identify  the  available  databases  and  additional  data 
necessary  for  further  improvements, 

o  Address  the  operational  standardization/compatibility 
issues , 

o  Consolidate  the  near-term  and  long-term  research  and 
development  activities,  inter-agency  cooperation/ 
collaboration  and  funding  support,  and 

o  Generate  and  document  recommendations  and  action 
items. 

The  conference  was  held  at  the  Raintree  Inn,  Airport,  in 
Colorado  Springs,  Colorado,  on  22  and  23  March  1988,  and  was 
well  attended. 


JOSEPH  J.  F.  LIU 
GEORGE  R.  DAVENPORT 
Co-Chairmen 
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GENERAL  EARL  VAN  INWEGEN'S  REMARKS 
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TO 

NEUTRAL  DENSITY  SPECIALITY  CONFERENCE 


22  MARCH  1988 


WELCOME  ...  I'M  DELIGHTED  TO  HAVE  THE  OPPORTUNITY  TO  WELCOME  YOU 
TO  COLORADO  SPRINGS  AND  TO  SPACE  COMMAND. 


MAKING  OPENING  REMARKS  ON  THE  SECOND  DAY  OF  A  CONFERENCE  IS 
INDICATIVE  OF  THE  SIGN  OF  THE  TIMES.  THE  PRESS  OF  THE  CURRENT 
BUDGET  PROBLEM  ...  OF  HAVING  TO  MINIMIZE  THE  OPERATIONAL  IMPACT 
OF  DOING  LESS  WITH  LESS  . . . 


AT  LEAST  I'M  HERE  ON  THE  SECOND  DAY!!! 


ACCEPTING  THE  FACT  THAT  WE  WILL  HAVE  TO  DO  LESS  WITH  LESS  IS  A 
HARD  BUT  NECESSARY  PILL  TO  SWALLOW.  IT  IS  UP  TO  US  ...  ALL  OF 
US  ...  TO  MINIMIZE  THE  OPERATIONAL  IMPACT  OF  HAVING  LESS. 


IN  ASTRODYNAMICS  WE  HAVE  AN  OPPORTUNITY  TO  MAYBE  NOT  HAVE  TO 
DECIDE  HOW  WE  WILL  DO  LESS  ...  EVEN  THOUGH  WE  GET  LESS  ...  BECAUSE 
I  BELIEVE  THAT  IF  WE  WORK  TOGETHER  AND  FOCUS  OUR  EFFORTS,  WE 
MIGHT  BE  ABLE  TO  DO  MORE  ...  EVEN  THOUGH  OUR  BUDGETS  ARE 
SHRINKING  IN  REAL  TERMS. 


IN  A  VERY  REAL  SENSE,  THAT  IS  THE  PURPOSE  OF  THIS  CONFERENCE  ... 


WHAT  CAN  WE  DO  BETTER  BY  WORKING  TOGETHER? 
WHAT  CAN  BE  ACHIEVED  IF  WE  COOPERATE? 


THIS  IS  THE  FIRST  CONFERENCE  WE,  AIR  FORCE  SPACE  COMMAND,  HAVE 
SPONSORED  AS  PART  OF  THE  NEW  DIRECTION  WE  ARE  TAKING  IN  DEALING 
WITH  OPERATIONAL  ASTRODYNAMICS. 


WE  HAVE  THREE  MAJOR  UNDERTAKINGS  IN  OUR  SPACE  COMMAND  OPERATIONS, 
DIRECTORATE  OF  OPERATIONS  ANALYSIS  ...  WHICH  MANY  OF  YOU  MIGHT 
REMEMBER  AS  OUR  DIRECTORATE  OF  ASTRODYNAMICS  ...  FOR  1988. 
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THE  FIRST  IS  A  COMPREHENSIVE  REVIEW  OF  THE  ERRORS  IN  OUR 
ASTRODYNAMIC  ORBIT  DETERMINATION  AND  PROPAGATION  METHODOLOGIES. 
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IF  WE  CAN  ACCURATELY  CHARACTERIZE  THE  ERROR  SOURCES  AND  THE 
IMPACT  OF  THOSE  ERRORS  ON  OUR  OPERATIONS  ...  WE  HAVE  THE  KEY  FOR 
GETTING  STRONG  SUPPORT  FOR  RESEARCH  AND  DEVELOPMENT  IN 
ASTRODYNAMICS . 


WITHOUT  THAT  FOUNDATION  I  BELIEVE  IT  WILL  BE  VIRTUALLY  IMPOSSIBLE 
TO  WIN  IN  THE  BUDGET  BATTLE. 


OUR  SECOND  MAJOR  INITIATIVE  IS  DEVELOPMENT  OF  AN  ASTRODYNAMIC 
STANDARD  FOR  MILITARY  SPACE  OPERATIONS. 


WE  ALL  STARTED  WTTH  THE  SAME  PHYSICS,  THE  SAME  DATA  BASES,  AND 
THE  SAME  MATH  A  NUMBER  OF  YEARS  AGO.  THINGS  WERE  SO  MUCH  SIMPLER 
THEN. 


WHEN  THERE  WERE  ONLY  A  FEW  SATELLITE  PROGRAMS  AND  THE  TOTAL 
NUMBER  OF  MAN  MADE  OBJECTS  IN  ORBIT  WERE  NUMBERED  IN  THE  HUNDREDS 
INSTEAD  OF  THE  THOUSANDS  AMD  TENS  OF  THOUSANDS,  WE  COULD  TALK  TO 
EACH  OTHER. 


SPACE  HAS,  HOWEVER,  GROWN  UP  AND  NOW  IT  IT  COMPLEX.  IN  A  NUMER 
OF  INSTANCES,  THAT  COMPLEXITY  IS  NOT  AN  ESSENTIAL  CHARACTERISTIC 
.  .  .  IT  IS  MAN  MADE . 


WE  ARE  LOOKING  IN  OUR  OPERATIONS  ANALYSIS  SHOP  AT  HOW  TO  STEP 
AWAY  FROM  SOME  OF  THAT  COMPLEXITY.  IN  SO  DOING  WE  HOPE  TL 
IMPROVE  PERFORMANCE  AND  REDUCE  COST  AT  THE  SAME  TIME. 


OUR  THIRD  MAJOR  INITIATIVE  IS  A  SERIES  OF  PROGRAMS,  OF  WHICH  THIS 
CONFERENCE  IS  THE  FIRST,  TO  TRANSFORM  THE  GOOD  IDEAS  FROM  THE 
ERROR  ANALYSIS  AND  STANDARDS  DEFINITION  EFFORTS  INTO  REALITY. 


AS  YOU  KNOW,  WE  FIX  THE  POSITION  OF  AN  INCREDIBLE  NUMBER  OF  NON- 
COOPERATIVE  ORBITAL  OBJECTS  EACH  DAY  USING  A  VARIETY  OF  SENSORS 
OF  VARYING  OUALITY.  THIS  HAS  LED  US  AND  OTHERS  IN  OUR  BUSINESS 
TO  MAKE  A  NUMBER  OF  APPROXIMATIONS  TO  BOOST  COMPUTATIONAL  SPEED, 
TO  COMPENSATE  FOR  THE  QUALITY  OF  OUR  OBSERVATIONS,  AND  TO 
COMPENSATE  FOR  POOR  INPUT  DATA. 
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ONE  INPUT  DATA  THAT  HAS  BEEN  AND  CONTINUES  TO  BE  WEAK  IS  NEUTRAL 
DENSITY. 


WE  RECOGNIZE  THAT  CONSIDERABLE  GOOD,  HIGH  QUALITY  RESEARCH  HAS 
BEEN  PERFORMED  ON  THE  NEUTRAL  DENSITY  PROBLEM  OVER  THE  LAST 
TWENTY  YEARS.  THE  CONTRIBUTIONS  YOU  HAVE  MADE  ON  SPECIFIC 
PROBLEMS  ARE  PARTICULARILY  NOTEWORTHY  BECAUSE  OF  THE  DIFFICULTIES 
ASSOCIATED  WITH  HIGH  ATMOSPHERE  RESEARCH.  I  APPRECIATE  THE  FACT 
THAT  WE  HAVE  NOT  FOUND  A  WAY  TO  KEEP  A  SENSOR  IN  THE  50  TO  250 
KILOMETER  ALTITUDE  RANGE  FOR  ANY  EXTENDED  PERIOD  OF  TIME.  THAT 
FORCES  YOU  TO  WORK  A  TOUGH  PROBLEM  WITH  VERY  LIMITED  DATA. 


OUR  PROBLEM  IS  THAT  FOR  MOST  OF  THAT  TWENTY  YEARS  THERE  HAS  NOT 
BEEN  A  SIGNIFICANT  IMPROVEMENT  IN  NEUTRAL  DENSITY  MODELING  FOR 
MILITARY  SPACE  OPERATIONS  ...  AND  YET,  IN  THE  LOW  ALTITUDE 
DOMAIN,  NEUTRAL  DENSITY  IS  ONE  OF  THE  PRINCIPAL  ...  IF  NOT  THE 
PRINCIPAL  ...  SOURCE  OF  ORBIT  PROPAGATION  ERROR. 


WE,  THE  SPACE  OPERATORS  ...  AIR  FORCE  SPACE  COMMAND  AND  THE 
UNITED  STATES  SPACE  COMMAND  NEED  BETTER  NEUTRAL  DENSITY  MODELS  TO 
DO  TODAY'S  MISSIONS  BETTER  AND  MAKE  TOMORROW'S  MISSIONS  POSSIBLE. 


SO  . . .  WHAT  DO  WE  DO  IN  AN  ENVIRONMENT  WHERE  MONIES  FOR  RESEARCH 
AND  DEVELOPMENT  ARE  GOING  AWAY? 


HOW  DO  WE  DO  MORE? 

HOW  DO  WE  DO  MORE  WITH  LESS? 

HOW  DO  WE  PREPARE  FOR  THE  FUTURE? 


I  BELIEVE  THE  KEY  TO  DOING  MORE  WITH  LESS  IS  CLOSE  COOPERATION 
AND  COORDINATION  BETWEEN  THE  RESEARCH  AND  DEVELOPMENT  COMMUNITIES 
AND  THE  CUSTOMERS  ...  THAT'S  US,  THAT'S  THE  SPACE  OPERATORS. 


WITH  SUCH  CLOSE  COOPERATION  THE  RESEARCH  COMMUNITIES  WILL  BETTER 
UNDERSTAND  THE  OPERATIONAL  PROBLEMS  WE  FACE: 


HOW  DO  YOU  ACHIEVE  THE  REQUIRED  ACCURACY  WITH  LESS  THAN 
PERFECT  INPUT  DATA? 


HOW  DO  YOU  HANDLE  NON-COOPERATIVE  OBJECTS? 
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HOW  DO  YOU  minimize;  the  computational  process  so  we  can  live  with 
AVAILABLE  COMPUTERS? 


AT  THE  SAME  TIME,  WE  THE  OPERATORS  WILL  BETTER  UNDERSTAND  HOW  TO 
USE  THE  TOOLS  YOU  PROVIDE  US. 


WE  WILL  ALSO  BE  IN  A  STRONG  POSITION  TO  PROMOTE,  PARTICIPATE  IN, 
AND  SUPPORT  CRITICAL  RESEARCH  AND  DEVELOPMENT  EFFORTS.  WE  MAY 
EVEN  BE  IN  A  POSITION  TO  PROMOTE  NEW,  DEDICATED  HIGH  ATMOSPHERE 
SENSORS  TO  ELIMINATE  SOME  OF  THE  PROBLEMS  WE  HAVE  WITH  THE  DATA 
TODAY. 


ONLY  TIME  WILL  TELL.  THIS  CONFERENCE  IS,  HOWEVER,  A  CRITICAL 
FIRST  STEP.  IT  WILL  BE  FOLLOWED  BY  A  MAJOR  INITIATIVE  ON 
ASTRODYNAMIC  STANDARDS  ...  SO  WE  ALL  CAN  WORK  EFFICIENTLY  AND 
EFFECTIVELY  TOGETHER. 


WE  CAM  NOT  AFFORD  AN  "ASTRO  TOWER  OF  BABLE"  ...  AND  YET  THAT  IS 
ALMOST  WHAT  WE  HAVE  TODAY. 


SOMETIME  SOON  WE  WILL  PUBLISH  OUR  ASTRODYNAM ICS  STANDARDS  PACKAGE 
FOR  YOUR  REVIEW  AND  COMMENT.  THAT  WILL  BE  FOLLOWED  BY  A  DETAILED 
REPORT  ON  ERROR  BUDGETS  IN  ASTRODYNAM I CS . 


Win-!  THOSE  INITIATIVES  IN  PLACE,  WE  SHOULD  COLLECTIVELY  BE  ABLE 
TO  DEFINE  AND  DEFEND  IN  THE  BUDGET  PROCESS  THE  CRITICAL 
AS  TRODYMAN ICS  RESEARCH  AND  DEVELOPMENT.  AND,  IN  SO  DOING, 
SIGNIFICANTLY  ENHANCE  FUTURE  MILITARY  SPACE  OPERATIONS. 


AND  ...  AS  I  SAID  IN  THE  BEGINNING  ...  I  AM  EXTREMELY  PLEASED  TO 
SEE  YOU  HERE  TODAY.  I  SINCERELY  WISH  YOU  SUCCESS  AS  YOU  DRAW 
THIS  CONFERENCE  TO  A  CLOSE  ...  BECAUSE  YOUR  SUCCESS  IS  OUR 
SUCCESS. 


BEFORE  YOU  LEAVE,  PLEASE  TAKE  A  FEW  MOMENTS  TO  REFLECT  ON  THIS 
EFFORT  AND  THEN  GTVE  US  FEED-BACK.  IF  YOU  ARE  NOT  SURE  WE  ARE  ON 
THE  RIGHT  TRACK,  LET  US  KNOW.  WE  NEED  YOUR  GOOD  IDEAS  AND  YOUR 
SUPPORT  ...  WE  NEED  YOUR  COOPERATION. 


CLOSE  COOPERATION  IS  OUR  ONLY  HOPE  TO  DO  MORE  WITH  LESS. 

NOW,  I  WILL  RE  GLAD  TO  TAKE  A  FEW  MINUTES  IF  YOU  HAVE  ANY  GENERAL 
OFFICER  TYPE  QUESTIONS  ... 
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1.  Introduction 


It  became  apparent  at  the  Neutral  Density  Workshop  held  at  the  Air  Force 
Geophysics  Lab  (AFGL)  during  October  1987,  that  many  members  of  the  research 
community  were  not  familiar  with  the  methods  used  by  Air  Weather  Service  (AWS) 
and  AFGL  to  evaluate  space  environmental  support  requirements,  seek  solutions 
to  operational  shortfalls,  and  finally  provide  operational  support.  This  paper 
provides  the  overview  of  AWS  space  environmental  support  and  how  the  neutral 
modeling  requirements  established  by  the  space  system  operators  that  AWS 
supports  fit  into  this  system  (Figure  1).  After  providing  an  understanding  of 
AWS  space  environmental  support,  the  operational  support  requirements  and  our 
neutral  modeling  plan  are  presented. 

2.  AWS  Space  Environmental  Support 

AWS  is  one  of  only  two  US  government  agencies  with  the  responsibility  to 
provide  basic,  operational  space  environmental  support  (Figure  2).  The  other 
organization  Is  the  National  Oceanographic  and  Atmospheric  Administration 
(NOAA) .  AWS  and  NOAA  try  to  complement  each  other's  observational  and  modeling 
capabilities  so  that  unnecessary  duplication  Is  minimized.  For  example, 
observations  of  the  space  environment  made  by  AWS  and  NOAA  systems  are  shared. 
Furthermore,  the  output  of  models  run  by  NOAA  is  available  to  AWS,  and  plans 
are  being  made  to  make  the  output  of  models  currently  under  development  for  AWS 
available  to  NOAA.  This  spirit  of  cooperation  has  served  both  organizations 
wel  1 . 

Using  the  shared  data  base  mentioned  above,  NOAA  has  the  ability  to  provide 
operational  space  environmental  support  to  civilian  agencies  and  individuals. 
One  such  civilian  agency  which  may  require  NOAA  support  is  the  National 
Aeronautics  and  Space  Administration  (NASA).  Meanwhile,  AWS  uses  this  same 
shared  data  base  to  provide  basic  space  environmental  support  to  all  Department 
of  Defense  (DOD)  organizations.  This  does  not  mean  that  other  DOD 
organizations  in  the  Army,  Navy,  or  Air  Force  do  not  tailor  AWS  space 
environmental  support  products  to  their  particular  space  system;  it  simply 
means  that  the  tailored  support  is  based  on  a  forecast,  alert,  or  analysis 
which  originated  from  AWS. 

In  order  for  AWS  to  provide  support  to  a  space  system  operator,  operational 
space  environmental  support  requirements  must  be  stated.  Then,  based  on  these 
stated  requirements,  AWS  can  determine  whether  an  operational  shortfall  exists. 
If  the  answer  is  yes  and  research  is  required  to  address  the  shortfall,  then 
AWS  generates  a  geophysical  requirement  (GR)  (Fiqure  3).  In  response  to  the 
GR,  AFGL  will  evaluate  the  relevant  research  (6.1)  and  exploratory  development 
(6.2)  that  has  either  been  completed  or  is  in  progress.  Applicable  6.1 


research  and  6.2  development  can  either  be  monitored  or  financially  supported 
by  AFGL  with  the  goal  of  getting  this  research  to  the  stage  where  it  has 
matured  and  is  ready  to  be  transitioned  into  an  operational  algorithm.  If  the 
research  has  matured  to  a  point  where  operational  shortfalls  can  be  satisfied, 
then  AWS  begins  the  program  acquisition  process.  The  first,  step  is  to  .jeneral- 
a  statement  of  operational  need  (SON).  The  Space  Env i ronmental  Technology 
Transition  (SETT)  SON  defines  AWS  space  modeling  needs.  The  methqd  jsed  to 
develop  an  operational  algorithm  1  tne  technology  transition  (TL)  program 
designed  to  use  6.3  advanced  development  followed  by  3400  operations  and 
maintenance  operational  software  development  to  make  the  bridge  between  6.2 
exploratory  development  and  an  AWS  operational  capability.  The  technology 
transition  program  is  managed 2Joi ntly  by  AWS  and  AFGL,  with  AFGL  maintaining 
technical  interaction  with  Tc  model  developers  (Figure  4)  and  evaluating 
progress  on  development  of  the  required  model.  Headquarters  AWS  monitors  the 
progress  AFGL  is  making  on  model  development  and  ensures  that  operational  needs 
are  addressed  during  development. 

At  this  time  there  are  two  organizations  providing  space  envi ronmental 
support:  the  Air  Force  Global  Weather  Central  (AFGWC)  located  at  Offutt  AFB 
NE,  and  4th  Weather  Wing  (4  WW)  located  at  Peterson  AFB  CO  (Figure  4).  AFGWC 
currently  provides  all  centralized  space  environmental  support.  All  available 
ground-based  and  space-based  measurements  are  used  by  AFGWC,  as  well  as 
computer  models  of  the  near-Earth  space  environment.  The  space  forecast 
products  of  the  AFGWC  support  many  space  system  customers,  in  addition  to  the 
US  Space  Command  (USSPACECOM)  community,  whose  needs  we  are  addressing  here. 

In  1984,  4  WW  became  the  main  AWS  ogani cation  with  space  env i ronmenta ' 
support  In  its  charter.  Currently,  4  WW  tailors  centralized  support  from  AFGWC 
to  the  units  it  directly  supports,  such  as  USSPACECOM,  Air  Force  Space  Commind 
( AFSPACECOM) ,  and  North  American  Aerospace  Defense  Command  (NORAD).  In  ;he 
near  future,  this  will  change.  A  new  centralized  forecast  facility  dedicated 
to  space  environmental  support  is  being  built  at  Falcon  AFB  CO.  This  unit  will 
be  called  the  Space  Forecast  Center  (SFC)  and  will  be  subordinate  to  4  WW 
(Figure  5).  Once  SFC  is  fully  operational,  all  the  centralized  space 
environmental  support  currently  done  at  the  AFGWC  will  be  transferred  to  the 
SFC.  In  addition,  all  the  new  space  environmental  support  moc  is  being 
developed  will  also  run  at  the  SFC.  Moving  all  the  centralized  spa.  e  support 
to  Colorado  has  several  advantages.  Since  US  space  operations  are  managed  in 
Colorado,  management  of  space  environmental  support  (by  4  WW)  and  the  actual 
centralized  support  (by  the  SFC)  should  be  dose  to  the  management  of  snare 
operations,  in  order  to  be  responsive. 

The  major  areas  where  AFGWC  currently  provide  space  environmental  support 
are  geomagnetic,  ionospheric,  and  spacecraft  (Figure  6).  The  basic  spare 
environmental  data  and  model -produced  fields  are  then  tailored  to  meet  the 
operational  requirements  of  different  space  systems.  In  the  ionospheric 
support  area,  for  example,  three-dimensiona1  electron  density  fields  may  be 
converted  to  maps  of  critical  frequencies  for  high-frequency  communications 
systems  operators.  While  the  areas  of  current  support  are  indicated  by  solid 
lines,  the  area  where  we  intend  to  provide  future  support  is  the  neutral 
thermosphere,  indicated  by  dashed  lines. 

There  are  two  types  of  neutral  density  support  required  by  AWS  customers. 
The  first  type  involves  forecasting  the  effects  of  the  atmosphere  on  spacecraft 


orbits,  particularly  for  satellites  in  low-Earth  orbit.  The  other  requirement 
is  to  provide  neutral  variable  specification  at  a  point.  At  this  time,  AWS 
does  not  provide  either  type  of  neutral  atmospheric  support.  Instead,  AWS 
provides  forecasts  of  both  geomagnetic  indices  and  solar  radio  flux  at  10.7  cm. 
The  F10.7  flux  is  a  proxy  indicator  of  the  solar  extreme  ultraviolet  (EUV) 
radiation,  which  results  in  dayside  thermospheric  heating.  These  indices  are 
then  used  by  space  system  operators  as  input  into  neutral  models  which  they  run 
for  themselves.  For  NORAD  this  consists  of  Jacchia's  1964  model. 

Data  used  ty  AWS  to  provide  space  environmental  support  is  collected  from 
ground-based  and  space-based  sources  (Figure  7).  The  ground-based  data  sources 
include  five  different  data  measurement  networks.  Likewise,  the  space-based 
measurement  systems  are  shown  in  Figure  8.  Note  that  at  this  time,  there  are 
no  measurement  systems  capable  of  providing  either  remote  or  in-situ  neutral 
measurements  to  AWS.  That  shortfall  has  several  implications.  First,  there  is 
no  way  for  AWS  to  use  neutral  measurements  either  to  initialize  or  to  verify 
the  planned  neutral  model.  Further,  as  we  attempt  to  build  a  comprehensive 
understanding  of  the  near-Earth  space  environment  through  both  measurements  and 
model,  there  is  at  least  99  percent  of  this  noar-Earth  space  environment  that 
we  currently  are  not  observing. 

AWS  is  acquiring  a  series  of  space  environmental  models  (Figure  9).  While 
each  model  can  function  in  a  stand-alone  mode,  it  is  logical  to  couple  the 
input  and  output  fields  of  the  models  in  a  manner  similar  to  the  way  these 
regions  are  coupled  in  nature.  Currently,  the  AWS  ionospheric  model 
operational  but  undergoing  evolutionary  upgrades.  The  magnetospheri c  model  is 
under  development  and  will  become  operational  in  fiscal  year  1990  (FY90). 
Further,  we  expect  to  begin  advanced  development  of  the  proposed  neutral  model 
in  FY89.  Work  on  integration  of  the  models  and  development  of  a 
sol  ar/ i nterpl anetary  model  is  expected  to  begin  in  FY90.  When  ail  the 
individual  models  are  completed  and  integrated,  AWS  will  have  the  first  known 
space  model  that  couples  the  physics  of  these  different  regions  in  a  dynamic 
manner  all  the  way  back  to  the  driving  source,  the  Sun. 

To  look  now  in  more  detail  at  the  area  of  interest,  a  block  diagram  of  the 
neutral  atmospheric  model  is  shown  in  Figure  10.  The  main  data  sources  for  the 
model  are  shown  at  the  top.  These  input  data  are  reduced  to  a  form  usable  by 
the  model.  At  the  same  time,  the  neutral  model  is  capable  of  accepting  input 
from  both  the  magnetospheric  and  ionospheric  models.  The  neutral  model  will 
produce  output  fields  of  density,  winds,  temperature,  and  composition. 
User-speci f ic  application  software  will  then  be  developed  to  address  specific 
customer  support  requirements. 

3.  AWS  Neutral  Modeling  Requirements 

Now  that  we  have  seen  an  overview  of  how  AWS  provides  space  environmental 
support,  it's  appropriate  to  look  at  the  specific  neutral  atmospheric  support 
requirements.  Figure  11  provides  these  details  by  first  showing  the  two  groups 
of  customers  requiring  support. 

These  requirements  are  quite  simple.  The  USSPACECOM/AFSPACECOM/NORAD 
community  requires  forecasts  of  the  neutral  atmosphere  out  to  72  hours  while 
other  00D  space  system  operators  require  neutral  atmospheric  speci f ication  at 
particular  locations  or  points. 


» 


The  output  variables  for  both  customer  groups  are  the  same:  neutral 
density,  neutral  winds,  neutral  temperatures,  and  neutral  composition.  The 
consolidated  altitude  range  and  accuracy  goals  of  interest  are  also  shown. 

4.  Neutral  Modeling  Status 

An  appropriate  question  following  exposure  to  the  neutral  modeling 
requirements  could  be,  "Where  do  we  stand  on  this  program?"  These  details  are 
outlined  in  figure  12.  Frank  Marcos,  AFGL/LIS,  is  the  manager  for  model 
development.  To  make  maximum  use  of  this  position,  he  organized  and  ran  the 
Neutral  Density  Workshop  held  at  AFGL  in  October  1987.  During  this  meeting, 
the  state-of-the-science  was  reviewed  and  the  science  team  formed.  A 
subsequent  meeting  of  the  neutral  modeling  science  team  was  held  in  Boulder  CO 
in  January  1988.  The  next  logical  step  was  to  present  the  joint  AWS-AFGL 
neutral  modeling  proposal  to  one  of  the  operational  groups  stating  neutral 
requirements,  the  USSPACECOM/AFSPACECOM/NORAD  community.  That's  the  purpose  of 
this  meeting. 

Simultaneously,  we  intended  to  check  for  the  existence  of  other  neutral 
requirements  or  refinements  to  current  requirements.  Since  we  intend  to  begin 
development  on  this  model  early  in  FY89,  we  do  not  intend  to  incorporate  major 
changes  to  the  already  stated  operational  requirements  after  May  1988  (Figure 
13).  For  that  reason,  if  any  DOD  space  system  operators  see  the  need  to  revise 
or  establish  additional  neutral  model  requirements,  now  is  the  time. 

Since  the  proposed  neutral  model  has  been  extensively  compared  with 
available  data  in  the  250-500  km  altitude  range,  this  is  the  region  where 
confidence  is  highest.  Therefore,  when  the  Phase  I  initial  operating 
capability  (IOC)  is  declared  in  FY91  (Figure  12),  the  proposed  neutral  model 
will  supply  reliable  data  in  this  altitude  range.  The  proposed  model  will  also 
be  used  to  supply  neutral  data  in  the  500-1500  km  altitude  range,  but 
confidence  in  the  output  variable  will  be  lower.  Values  of  density, 
temperature,  and  composition  below  250  km  during  Phase  I  will  be  supplied  by 
the  mass  spectrometer  and  incoherent  scatter  (MSIS)  model. 

Follow-on  phases  of  neutral  model  development  will  extend  the  proposed 
neutral  model  down  to  90  km  as  databases  for  verifying  the  model  performance  in 
the  90-250  km  altitude  range  mature. 

5.  Conclusion 


In  summary,  we  have  presented  an  overview  of  Air  Weather  Service  space 
envi ronmental  support  responsibilities.  In  addition,  we  have  seen  how 
operational  shortfalls  are  addressed,  the  relationship  between  technical  and 
operational  space  environmental  support  organi zations ,  and  the  kinds  of  space 
environmental  support  provided  by  the  Air  Force  Global  Weather  Central  and 
later  the  Space  Forecast  Center.  We  have  also  reviewed  the  sources  of  data  and 
overall  modeling  plan.  In  more  detail,  we  looked  at  the  neutral  atmospheric 
customers,  their  requirements,  and  how  AWS/AFGL  jointly  plan  to  satisfy  these 
requirements . 

In  conclusion,  we  know  that  an  accurate  neutral  atmospheric  model  is  needed 
to  satisfy  stated  requirements  and  believe  we've  chosen  the  best  approach 
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possible  based  on  anticipated  resources  (Figure  14).  However,  if  there  are  any 
00D  space  system  operators  who  need  to  refine  these  neutral  requirements  or 
levy  additional  requirements,  now  is  the  time  to  do  so. 
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Figure  14 


Air  Force  Space  Command  Requirements  For  Neutral  Atmospheric 
Density  Specification  and  Prediction 

Lt  Col  George  R.  Davenport 
Headquarters  Fourth  Weather  Wing 


The  Air  Force  Space  Command  Deputy  Director  of  Astro- 
dynamics  last  stated  the  Ai r  Force  Space  Command  "Requirements 
for  Atmospheric  Neutral  Density  Predictions"  in  a  26  flay  1986 
letter  to  the  4th  Weather  Wing  Aerospace  Sciences  Division.  The 
requirement  was  for  predictions  of  the  neutral  atmospheric  den¬ 
sity  for  24  hours  out  to  ?2  hours,  presumably  to  be  updated 
every  24  hours,  at  ail  latitudes  and  longitudes  for  the  100  km 
to  lOOO  km  altitude  range.  The  upper  limit  of  the  altitude 
range  is  likely  to  be  raised,  possibly  to  1  500  km.  The  accuracy 
requirement  is  in  the  range  of  5%. 

Air  Force  Space  Command  also  made  several  specific  recom¬ 
mendations  in  the  26  May  1986  letter.  The  first  recommendation 
is  to  "Create  brand  new  models  which  utilize  additional  environ 
mental  input  parameters. "  Additional  model  inputs  might  include 
the  solar  x-ray  and  solar  extreme  ultraviolet  fluxes,  auroral 
energy  inputs,  the  solar  wind,  and  measurements  of  the  inter¬ 
planetary  magnetic  field  The  second  recommendation  is  to  "Get 
instruments  in  p 1 aco  for  measuring  and  predicting  new  environ 
mental  parameters  that  may  be  nee>ied .  "  Such  instruments  might 
fly  on  the  Defense  tleteorological  Satellite  Program  (DMSP) 
spacecraft,  or  chey  could  include  data  from  the  NCAA  Solar  > 

Ray  Imager,  the  proposed  NASA  WlND/Sl'IPL  libration  point 
spacecraft,  or  a  follow-on  to  the  Navy's  experimental  Remote- 
Atmospheric  and  Ionospheric  Detection  System  (RAIDS).  The  thi rd 
recommendation  is  to  "improve  measurements  and  forecasts  for 
currently  used  environmental  parameters. "  These  parameters  are 
the  geomagnetic  indicies  and  the  2800  MHz,  10.  ?  centii  ’ter, 
solar  radio  flux  measured  daily  from  the  Algonquin  Obs i  vatory, 
the  "Ottawa  Flux"  or  "FIO. " 

The  Air  Weather  Service  established  Geophysical  Requirement 
2-8?,  "Neutral  Atmospheric  Specification  and  Prediction, "  to 
help  the  research  community  to  understand  shortcomings  in 
generating  products  that  describe  the  current  and  the  future 
conditions  of  the  neutral  atmosphere  to  operational  customers. 
Geophysical  Requirement  2-8?  includes  Air  Force  Space  Command 
requirements  and  recommendations  as  stated  in  the  21  May  1986 
letter. 

Currently  Air  Weather  Service  is  just  providing  observed 
and  predicted  surrogate  environmental  parameters  for  each  user 
to  input  into  their  own  neutral  atmosphere  model.  The  Ottawa 
Flux  and  the  geomagnetic  indices  are  all  that  is  available,  so 
providing  these  parameters  constitutes  the  current  level  of 


* 


ft 


support.  This  is  analogous  to  providing  an  aircraft  pilot  with 
raw  weather  data  for  him  to  use  to  prepare  his  own  flight  route 
forecast . 

The  goai  for  the  next  few  months  is  to  update  the  Air  Force 
Space  Command  requirement  that  was  last  stated  almost  two  years 
ago.  If  appropriate,  the  U.S.  Space  Command,  Air  Force  Space 
Command,  and  Navy  Space  Command  should  publish  a  formal  Required 
Operational  Capability  (ROC)  with  supporting  Statements  of 
Operational  Need  (SONs).  In  the  process  of  preparing  and 
coordinating  these  requirements  documents,  the  commands  should 
stimulate  similar  stated  requirements  from  other  DOD  agencies. 
For  example.  Strategic  Air  Command  identified  impacts  of 
increased  atmospheric  density  in  their  Statement  of  Operational 
Need  for  Military  Space  Flight  Capability,  which  was  validated 
on  22  September  1980.  Although  an  old  SON,  it  indicates  that 
there  is  at  least  hope  for  support  from  other  Air  Force  major 
commands. 
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Naval  Space  Surveillance  Center  Atmospheric  Drag  Requirements 


Stephen  H.  Knowles 
Naval  Space  Surveillance  Center 
Dahlgren,  Va.  22448 


The  Naval  Space  Surveillance  Center,  with  headquarters 
located  in  Dahlgren,  Va.,  maintains  a  complete  catalog  of 
positions  of  earth-orbiting  space  objects  as  part  of  its 
operational  responsibility  as  a  component  of  the  U.S.  Space 
Command.  As  such,  it  has  a  need  to  know  the  atmospheric  drag  to 
an  accuracy  comparable  to  the  nominal  accuracy  with  which  orbits 
are  maintained.  The  accuracy  of  individual  orbits  in  the 
NAVSPASUR  catalog  varies  considerably  in  a  way  that  depends  on  a 
number  of  parameters,  including  the  altitude  of  the  orbit,  the 
tasking  of  the  sensors,  the  stability  of  the  thermosphere,  etc. 
A  nominal  desired  accuracy  for  near-earth  orbits,  however,  is  of 
the  order  of  +/-  1  to  3  miles.  This  accuracy  is  relatively 
coarse  compared  to  precision  orbit  determinations  that  can 
produce  accuracies  of  a  few  meters.  Nevertheless,  the  inaccuracy 
associated  with  predicting  the  effect  of  atmospheric  drag  is 
usually  the  most  important  term  contributing  to  orbit  prediction 
degradation  with  time.  Thus,  improvements  in  drag  compensation 
are  of  high  importance  in  improving  NAVSPASUR  orbital  prediction 
ability.  This  is  of  even  more  importance  for  the  limited  but 


increasing  number,  of  orbits  for  which  precision  orbit  prediction 
is  required.  A  special  important  need  exists  for  more  accurate 
upper  atmospheric  modeling  for  the  prediction  of  atmospheric  re¬ 
entry.  Present  drag  calculation  accuracy  requires  near-real-time 
tracking,  and  accurate  advance  re-entry  predictions  arc  not 
possible 

Atmospheric  drag  is  a  dissipative  force.  Its  effects  on 
different  near-earth  orbits,  while  quantitatively  different,  are 
qualitatively  the  same,  namely,  a  decrease  in  the  orbital  energy. 
This  energy  decrease  is  exhibited  most  noticeable  as  a  constant 
rate  of  decrease  in  the  mean  motion.  This  results  in  a  change  in 
the  along-track  position  that  is  quadratic  because  of  the 
integration  involved.  The  mean  motion  effect  results  in  the 
typically  5  times  worse  accuracy  in  determining  along-track 
position  which  is  the  largest  single  uncertainty  in  most  near- 
earth  orbits.  A  second  important  effect  is  the  reduction  of 
eccentricity,-  orbits  tend  to  circularize  as  a  result  of  drag. 
Other  elements  describing  the  motion  of  the  satellite  within  its 
plane  are  also  perturbed,  including  a  secular  reduction  in  the 
semi-major  axis. 

Satellite  drag  is  largely  an  affect  of  that  region  of  the 
upper  atmosphere  known  as  the  thermosphere,  with  a  small 
contribution  from  the  lowest  part  of  the  exosphere.  The  density 
and  other  parameters  vary  greatly  over  this  region  (Figure  1). 
The  amount  of  along-track  atmospheric  effect  is  therefore  highly 
dependent  on  perigee  height.  Figure  2  gives  an  indication  of 
this  relationship;  above  1000  km.  atmospheric  drag  can  be 


35 


neglected  for  most  purposes.  Drag  can  be  several  tens  of  seconds 
per  day  for  near-reentry  situations,  and  more  typically  appears 
as  an  error  in  predicted  arrival  time  of  about  +/-  1  second. 

At  the  present  time,  NAVSPASUR  makes  no  correction  for 
atmospheric  drag  from  a  model  in  computing  our  production  orbits. 
The  effect  of  drag  is  included  empirically.,  however,  by  the 
inclusion  of  the  first  and  second  derivatives  of  the  mean  motion 
( ndot  and  n  double  dot)  as  orbital  elements  during  the 
differential  correction  fitting  process.  Our  differential 
correction  is  performed  once  per  day;  the  orbital  elements, 
including  the  drag  terms,  however,  are  the  average  over  the  data 
fitting  span  used  in  the  differential  correction.  This  can  vary 
from  3  to  10  days  for  near-earth  orbits,  but  is  typically  3  dayi 
for  most  near-earth  orbits  that  exhibit  appreciable  drag.  The 
lack  of  explicit  inclusion  of  a  model  has  several  effects,  both 
advantageous  and  disadvantageous.  In  the  first  place,  the 
correct  mean  atmospheric  drag  is  automatically  fitted  for  each 
orbit's  parameters,  without  any  necessity  to  include  an  error 
term  from  model  inaccuracy.  However,  there  is  no  ability  to 
include  dynamic  behavior  of  the  thermosphere.  During  the 
occurrence  of  major  solar  storms,  this  has  resulted  in  situations 
wherein  large  numbers  of  orbits  in  our  operation  center  have 
exceeded  the  tolerances  that  made  them  appear  as  'UCTs' 
(uncorrelated  targets),  necessitating  a  large  amount  of  labor  by 
analysis  staff  to  make  them  fit.  Some  of  this  problem  can  and 
should  clearly  be  overcome  by  use  of  an  atmospheric  drag  model  to 
make  the  orbits  measurable  in  terms  of  some  physical  parameter. 


Hopefully  this  can  be  done  in  a  way  which  will  retain  the 
empirical  advantages  of  the  present  NAVSPASUR  system  while  still 
improving  the  ability  to  respond  to  geophysical  changes  in  the 
thermosphere  density. 

Currently  available  upper  atmospheric  models  were  largely 
developed  in  the  early  1970' s  by  Jacchia  and  others.  The  general 
consensus  is  that  these  models  are  capable  of  predicting  density 
in  the  100  to  1000  km.  region  to  an  accuracy  of  about  15%. 
Attempts  to  produce  better  accuracy  than  this  have  not  so  far  met 
with  success.  There  are  several  problems  with  attempting  to 
produce  improved  accuracy.  The  present  models  are  parametrized 
in  terms  of  global  parameters  having  to  do  with  solar  and 
magnetic  field  activity,  specifically  solar  10  cm.  flux  level  and 
the  Kp  magnetic  field  index.  These  parameters  as  distributed  are 
after-the-fact  by  up  to  several  days  and  are  averaged  to  some 
extent  rather  than  instantaneous.  In  addition,  they  are  global 
in  the  sense  of  having  only  one  parameter  to  represent  the 
overall  atmospheric  performance.  The  solar/geophysical 
indicators  serve  basically  as  indicators  of  energy  inputs  which 
increase  the  thermospheric/exospheric  temperature  and  thus  the 
atmospheric  density.  The  relation  between  exospheric  temperature 
and  density  is  not  simple  (Figure  3). 

As  an  operational  orbit  determiner/predictor,  NAVSPASUR 's 
most  pressing  need  for  improved  treatment  of  drag  is  to  account 
better  for  the  effect  of  solar/geomagnetic  disturbances 
(’storms')  on  the  catalog.  In  order  to  compensate  for  this 
effect,  real-time  data  inputs  must  be  available,  as  solar  storms 
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cannot  be  predicted  accurately  in  advance.  Conventional 
predictor  indices  do  not  serve  for  this  purpose;  an  additional 
problem  is  that  their  inherent  accuracy  is  considerably  poorer 
during  solar-geophysical  storms  than  it  is  during  quasi-steady- 
state  conditions.  At  present  there  is  no  fast-response 
monitoring  instrument  to  provide  immediate  notification  of  major 
solar  events.  Probably  the  best  method  to  do  this  is  a  solar 
ultraviolet  or  X-ray  detector;  one  is  planned  for  circa  1990  as 
an  adjunct  to  the  next  GOES  launch  by  NOAA.  There  is  also  a  case 
to  be  made  for  real-time  monitoring  of  the  upper  atmospheric 
density  at  various  locations  around  the  earth.  To  a  limited 
extent  this  can  be  accomplished  by  the  tracking  network  itself, 
but  not  as  presently  constituted  in  real  time.  Various  satellite 
systems  have  been  suggested,  including  a  series  of  'dragsats' 
which  would  instantaneously  measure  the  deceleration;  this  latter 
would  provide  the  best  current  reading  on  the  density  of  the 
thermosphere,  although  not  providing  a  thorough  sample  of  density 
at  various  altitudes.  There  is  presently  a  fairly  large  amount 
of  research  underway  to  improve  knowledge  of  the  composition  of 
the  thermosphere;  translating  this  into  improved  predictions  of 
orbiting  objects,  is  not  straightforward.  Some  of  the 
complexities  involved  are  illustrated  in  Figure  4,  which  provides 
an  indication  of  some  of  the  processes  controlling  thermospheric 
structure  and  dynamics.  It  seems  clear  that,  whatever 
improvement  in  thermospheric  density  modeling  is  decided  upon,  it 
will  include  provision  for  real-time  data  input  and  will  be  built 


upon  semi-empirical  formulas  rather  than  the 
elegant  series  expansions  of  orbital 


theoretically 
mechanics . 
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FIGURE  2 

Relation  between  atmospheric  density  and  temperature  for  various  heights 
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VELOCITY  REQUIREMENT  -  FEET/SECOND/ OAY 


FIGURE  3 

Velocity  required  to  maintain  orbital  height  as  a 
function  of  altitude 
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Principal  orocesses  and  energy  sources  governing 
thermospheric  structure  and  dynamic. 
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Space  Control  Operations  -  Meeting  the  Challenge 

by 

Capt  Colette  M.  de  la  Barre 

The  mission  of  the  USSPACECOM  Space  Control  Operations 
Centers  at  the  Cheyenne  Mountain  Air  Force  Station  is  space 
surveillance,  satellite  protection  and  negation.  In  recent 
years,  the  protection  aspect  of  this  mission  has  been  expanding. 
The  Strategic  Defense  Initiative  Organization,  Missile  Test 
Ranges,  and  NASA  require  more  accurate  space  surveillance  support 
for  many  varied  applications;  however,  the  one  that  has  received 
the  most  attention  recently  is  prediction  of  close  conjunctions 
between  satellites  and  launch  vehicles.  In  order  to  provide  more 
accurate  space  surveillance  products,  USSPACECOM  is  demanding  the 
most  from  its  personnel  and  the  spacetrack  network.  Since  these 
resources  are  limited,  USSPACECOM  must  look  at  other  ways  to  meet 
the  demands . 

The  Challenge 

One  of  the  challenges  facing  USSPACECOM  is  to  provide  more 
accurate  predictions  of  close  conjunctions  between  launch 
vehicles  and  orbiting  space  objects.  According  to  the  Outer 
Space  Treaty  of  1967,  launch  vehicle  owner/operators  are  held 
liable  if  their  launch  vehicle  collides  with  an  orbiting  payload. 
Although  launch  agency  personnel  are  also  concerned  with  avoiding 
orbiting  debris  to  protect  their  launch  vehicles,  they  usually 
resolve  to  avoid  only  payloads  in  order  to  minimize  the  impact  on 
their  launch  opportunities.  Since  the  prediction  of  the  location 
of  orbiting  debris  is  much  less  accurate  than  payloads,  the 
quality  of  USSPACECOM  generated  safe  launch  windows  are 
significantly  reduced  when  debris  must  be  avoided.  Moreover,  the 
NASA  Space  Shuttle  program  management  is  not  content  to  avoid 
only  payloads  due  to  the  threat  of  orbiting  debris  to  human  life. 
NASA's  requirements  provide  USSPACECOM  with  a  tremendous 
challenge  -  to  predict  all  catalogued  space  objects  with 
signif icantly  more  accuracy. 

Meeting  the  Challenge 

Many  factors  enter  into  the  accuracy  of  USSPACECOM  space 
surveillance  products  including  accuracy  of  spacetrack  sensors, 
frequency  of  spacetrack  observations,  accuracy  of  orbit 
propagation  models,  and  Space  Control  Operations  Center  personnel 
and  computer  work  loads.  An  improvement  in  any  of  these  areas 
would  improve  space  surveillance  products.  An  improvement  in 
accuracy  of  orbit  propagation  models  would  have  the  greatest 
impact,  because  it  would  result  in  more  accurate  space 
surveillance  products  while  reducing  the  time  spent  by  USSPACECOM 
personnel  and  the  spacetrack  network  on  routine  satellite 
maintenance.  The  resources  saved  could  be  put  to  use  in  other 
aspects  of  the  Space  Control  Operations  mission.  Since  the  Space 
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Control  Operations  Centers  don't  control  the  accuracy  of 
spacetrack  sensors  and  satellite  propagation  models,  but  have 
control  over  sensor  tasking  and  how  often  satellite  orbits  and 
space  surveillance  products  are  updated,  USSPACECOM  develops 
procedures  to  maximize  space  surveillance  product  accuracy  while 
not  stressing  their  computer  system.  These  procedures  result  in 
increased  personnel  work  load  and  increased  tasking  of  the 
spacetrack  network. 

The  accuracy  of  all  space  surveillance  products  depends  on 
how  accurately  the  Space  Control  Operations  Centers  can  predict 
the  location  of  a  satellite.  An  accurate  satellite  prediction  is 
generated  by  a  good  propagation  model  using  good  initial 
conditions.  Good  initial  conditions  are  generated  by  obtaining 
frequent  observations  on  the  satellite  from  the  spacetrack 
network  and  updating  the  satellite's  orbit  with  those 
observations  more  often.  Using  this  updated  satellite  orbital 
element  set,  a  more  accurate  space  surveillance  product  is 
generated  by  ensuring  the  time  span  over  which  the  product  is 
valid  is  not  excessive.  The  following  example  shows  how  the 
USSPACECOM  develops  procedures  and  uses  its  personnel  and  the 
spacetrack  network  to  generate  the  mcst  accurate  space 
surveillance  product.  To  ensure  the  best  accuracy  of  the  impact 
location  of  a  decaying  satellite,  the  spacetrack  sensors  are 
tasked  for  observations  on  the  satellite  every  time  it  comes 
within  their  field  of  view  and  the  satellite's  orbital  element 
set  is  manually  updated  about  every  eight  hours  using  these 
observations.  The  Tracking  and  Impact  Prediction  (TIP)  for  the 
satellite  is  updated  at  a  minimum  of  eight  times  in  15  days  with 
the  final  prediction  being  made  within  two  hours  of  impact.  If 
the  model  used  to  predict  the  satellite's  position  were  more 
accurate,  observations  could  be  required  less  frequently  by  the 
spacetrack  network,  the  satellite  orbital  element  set  would  not 
need  to  be  manually  updated  as  often,  and  the  TIP  for  the 
satellite  would  be  accurate  for  a  longer  period  of  time.  T"»e 
result  would  be  a  significant  decrease  in  resources  required  for 
TIPs. 


USSPACECOM  Atmospheric  Density  Models 

USSPACECOM  predicts  the  location  of  near  Earth  space  objects 
using  an  analytical  propagation  model  which  includes  a  power 
density  function  (GPPF)  and  a  numerical  integration  model  which 
includes  the  Jacchia-SCC  64  (JSCC)  empirical  atmospheric  density 
model  (SPJSCC)  [1].  Since  USSPACECOM  has  catalogued  over  18,900 
orbiting  space  objects,  and  currently  keeps  track  of  over  7100  of 
which  over  6300  are  near  Earth*,  computation  time  is  a  real 
consideration  in  selecting  which  model  to  use.  Since  the  SPJSCC 
model  requires  significantly  more  computer  time  than  the  GPPF,  it 
is  only  used  for  special  interest  space  objects  including 

*  A  near  Earth  satellite  has  an  orbital  period  less  than  or  equal 
to  225  minutes. 
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decaying  satellites  or  less  than  one  percent  of  all  near  Earth 
space  objects.  One  of  the  most  important  differences  between  the 
models  is  the  SPJSCC  model  includes  the  sun's  effect  on  the 
Earth's  atmospheric  density  and  the  GPPF  doesn't  [1].  Even 
though  the-  JSCC  model  is  more  sophisticated  than  the  power 
density  function,  it's  predictions  for  density  of  the  atmosphere 
are  still  reported  to  be  16  percent  off  from  atmospheric 
measurements.  Both  models  use  a  constant  satellite  ballistic 
coefficient  [ 1 ] . 

Performance 

The  USSPACECOM  Lost  Satellite  List  and  TIPs  for  decaying 
satellites  were  chosen  as  indicators  of  the  performance  of  the 
USSPACECOM  atmospheric  density  models.  A  near  Earth  satellite  is 
considered  lost  by  USSPACECOM  when  it  has  not  been  seen  by  the 
spacetrack  network  for  more  than  two  days.  There  are  many 
factors  that  affect  the  number  of  satellites  on  the  lost  list  and 
the  TIPs  for  decaying  satellites;  however,  a  significant 
contributor  is  the  geomagnetic  disturbance  level.  Figure  1  shows 
the  number  of  lost  satellites  from  Oct  86  to  Mar  88.  Figure  2 
shows  the  daily  geomagnetic  index  (Ap)  from  Jan  86  to  Mar  88. 
Analysis  of  the  lost  list  during  the  two  largest  geomagnetic 
storms  in  1988  (15  Jan  88  and  22  Feb  88)  shows  that  the  number  of 
lost  satellites  doubled  within  six  days  after  the  storms  (Fig.  3 
and  4.)  Since  the  satellites  on  the  lost  list  are  processed  with 
the  GPPF,  and  since  the  GPPF  includes  no  solar  effects,  it  is  not 
surprising  that  the  number  of  lost  satellites  increases 
significantly  during  high  geomagnetic  activity.  It  should  also 
be  noted  that  the  average  number  of  lost  satellites  over  the  17 
months  studied  is  178  which  is  not  consistent  with  protecting  the 
Space  Shuttle  and  other  launch  vehicles  from  collision  with 
orbiting  space  objects.  An  analysis  of  TIPs  also  shows  that  decay 
times  for  satellites  are  significantly  affected  by  geomagnetic 
disturbances . 


A  TIP  study  was  performed  that  analyzed  how  far  off  from  actual 
decay  time  the  two  hour  and  24  hour  TIP  times  were  for  a  five 
month  period  having  no  geomagnetic  storms  (Jan  86,  Mar  86,  Apr 
86,  Apr  87,  and  Jun  87)  and  a  four  month  period  having  at  least 
one  geomagnetic  storm  with  24  hour  geomagnetic  disturbance  index 
(Ap)  greater  than  80  (Feb  86,  Sep  86,  Jan  88,  and  Feb  88.)  Table 
1  shows  the  results  of  the  TIP  study.  The  two  hour  TIP  times  for 
the  25  satellites  that  decayed  during  the  months  with  no 
geomagnetic  storms  were  an  average  of  14  minutes  off  from  their 
actual  decay  times  while  the  24  hour  TIP  showed  a  one  hour  and  46 
minute  difference.  The  25  satellites  that  decayed  during  the 
months  with  large  geomagnetic  storms  were  found  to  be  an  average 
of  36  minutes  off  from  their  actual  decay  times  for  the  two  hour 
prediction  and  3  hours  and  42  minutes  off  for  the  24  hour 
prediction.  Since  a  decaying  satellite  can  travel  more  than 
6500  km  in  14  minutes,  even  the  two  hour  TIP  during  low 
geomagnetic  storm  levels  is  not  very  accurate.  The  TIP  made  24 
hours  prior  to  impact  is  of  little  value. 
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Recommendations 

USSPACECOM  will  continue  to  look  at  new  ways  to  use 
personnel  and  the  spacetrack  network  to  improve  space 
surveillance  support.  However,  improvements  to  USSPACECOM 
atmospheric  density  models  will  do  the  most  to  improve  space 
surveillance  support.  USSPACECOM  needs  a  satellite  propagation 
model  that  very  accurately  models  atmospheric  density  during  high 
and  low  solar  activity  without  requiring  more  computation  time 
than  the  GPPF.  This  model  could  be  used  for  all  near  Earth  space 
objects  except  special  interest  satellites.  An  improved  version 
of  SPJSCC  is  needed  to  process  special  interest  satellites 
including  decaying  satellites.  One  may  hypothesize  that 
improvements  in  predicting  the  indicators  of  geomagnetic  activity 
will  significantly  improve  the  products  generated  by  the  SPJSCC 
during  high  geomagnetic  activity.  However,  improvements  in 
modelling  atmospheric  density  during  low  solar  activity  are  also 
needed . 

New  Developments 

Dr  Joseph  Liu  and  Dr  Felix  Hoots  have  developed  some  new 
semianalytic  methods  for  propagating  satellite  orbits.  These 
methods  may  be  able  to  handle  improvements  in  analytical  and 
empirical  atmospheric  density  models  [2],  Mr  Peter  Kolodzie jczyk 
is  investigating  ways  to  predict  how  satellite  ballistic 
coefficient  varies  with  satellite  altitude  [3].  He  hypothesizes 
that  an  improved  prediction  model  could  absorb  inaccuracies  in 
atmospheric  density  models. 

Conclusions 

USSPACECOM  requires  improvements  in  near  Earth  atmospheric 
density  models  to  meet  the  demands  of  future  space  surveillance 
support.  without  improvements,  USSPACECOM  will  continue  to 
increase  the  burden  on  personnel  and  the  spacetrack  network. 

Since  USSPACECOM  resources  are  not  without  limit,  the  time  will 
come  when  USSPACECOM  will  have  difficulty  meeting  new  space 
surveillance  demands.  Modelling  atmospheric  density  simply 
enough  to  be  used  in  a  satellite  propagation  model  which  uses 
minimum  computer  time  is  no  trivial  problem.  Improving 
atmospheric  density  predictions  seems  to  increase  the  complexity 
of  the  mathematical  model  which  increases  the  computer  time 
required  to  propagate  a  satellite's  orbit.  A  "Newton-size" 
solution  is  needed  -  Are  th~re  any  Newton's  out  there? 
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FIGURE  1  -  Near  Earth  l.obt. 


FIGURE  2  -  Geomagnetic  Disturbances 


FIGURE  3  -  Near  Earth  Lost  Satellites  during  Jan  88 
Geomagnetic  Storm  - - - - - 
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Tabli  1  ■  Predicted  Decay  Tiae 

Versus  Actual  Decay  Tlae 


B 

24  Hr 

No  Geomag  Storms 
(Jan  86 ,  Mar  86 , 
Apr  86 ,  Apr  87 , 

J u  n  87 ) 

1+  Hi  n 

1  Hr  46  Min 

Geomag  Storms 
(Feb  86 ,  Sep  86 , 
Jan  88 ,  Feb  88 ) 

36  M  i  n 

3  Hr  42  Mm 

50  satellites 


14  Min  >  6500  k  m 
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AFSPACECOM  -  ATMOSPHERIC  NEUTRAL  DENSITY  SPECIALIST  CONFERENCE 

MARCH  22-23,  1988 
SESSION  2:  OPS  SUPPORT 

"MSFC'S  MODELING  AND  USE  OF  ORBITAL  DENSITY  RESULTS" 
by  DALE  JOHNSON,  NASA-MSFC/ED44/ (205)  544-1665 


The  Earth  Science  and  Applications  Division  (ESAD)  of  NASA/MSFC 
provides  orbital  atmospheric  support  for  the  design  studies  of 
various  planned  NASA  launch  and  space  vehicles.  Very  little 
operational  support  is  provided.  Only  as  needed.  Atmospheric 
models  such  as  the  MSFC/J7 0 ,  now  called  the  MSFC/Thermospheric 
Engineering  Atmosphere  Model  (TEAM) ,  are  utilized  in  response  to 
guidance,  attitude  and  control  studies  which  are  needed  by  design 
engineers  for  vehicle  response,  precision  orbit  selection  and 
positioning,  attitude  gyro  saturation  problems,  etc.  Atmospheric 
density  is  the  key  driver  in  these  studies.  A  second  area  of 
density  model  support  goes  to  the  NASA  lifetime 

prediction/orbital  decay/satellite  tracking  program  areas.  Also 
atmospheric  constituent  input  is  used  by  those  interested  in 
material  degradation  and  spacecraft  glow  problems.  These 
atmospheric  inputs  and  models  are  usually  documented  or  baselined 
as  the  design  criteria  to  use  on  selected  NASA  projects. 

Some  of  the  current  orbital  projects  supported  are:  Space 
Station  (document,  JSC  30425) ,  Space  Telescope,  OMV,  and  Tethered 
Satellite.  Current  orbital/middle  atmosphere  projects  include: 
Space  Shuttle,  AOTV,  AFE  and  NASP.  ESAD  generally  supports 
projects  on  the  drawing  board,  in  the  design  stage,  which  are 
usually  many  years  into  the  future.  Comparisons  of  various 
atmospheric  orbital  models  are  also  conducted  by  ESAD  in  order  to 
determine  and  use  the  proper  model  for  the  proper  application. 

See  Figure  1. 

A  similar  modeler/user  density  workshop  was  held  at  MSFC  in 
November  1985,  and  has  been  documented  as  NASA  CP-2460,  entitled, 
"Upper  and  Middle  Atmospheric  Density  Modeling  Requirements  For 
Spacecraft  Design  and  Operations."  I  intend  to  highlight  the  key 
conclusions  and  recommendations  as  determined  from  this  workshop. 
The  workshop  was  called  for  because  of  a  communications  gap 
between  the  modelers  and  the  engineering  user  communities.  Once 
models  or  criteria  were  released,  the  modelers  were  unaware  if 
the  data  and  variations  within  the  models  were  being  used 
correctly  or  not.  Likewise,  the  users  could  be  unaware  of 
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certain  limitations  and  restrictions  on  the  data  or  model.  This 
could  result  in  certain  design  and  operational  penalties.  See 
Figure  2.  There  were  eleven  issues/user  needs  which  surfaced  and 
were  concluded  as  a  result  of  this  workshop.  These  are  presented 
in  Figure  3. 

The  first  five  issues  represent  general  recommendations  and  those 
dealing  with  middle  atmospheric  needs.  The  remaining  six  issues 
deal  with  thermospheric  orbital  density  needs.  #6  Better 
prediction  of  the  input  parameters  used  (i.e.  solar  flux  and 
geomagnetic)  in  density  models  is  needed  in  the  short-or 
long-term.  #7  The  small-scale  density  structures  are  not 
accurately  modeled.  Perhaps  some  of  the  15%  density  model 
uncertainty  can  be  resolved  using  gravity  wave  theory,  etc.  #8  A 
need  for  more  accurate  thermospheric  wind  fields,  especially  at 
high  latitudes,  is  needed.  #9  Improvement  of  the  theoretical 
Thermospheric  General  Circulation  Models  (TGCM)  needs  to  be 
continued,  and  comparative  studies  done.  Modelers  need  to  be 
guided  when  using  TGCM's.  #10  For  satellite  drag  studies  the 
older  Jacchia  type  model  are  generally  preferred  -  since  drag 
data  went  into  these  models.  But  with  erosion  and  glow  studies 
today,  a  need  for  a  more  accurate  representation  of  the 
atmospheric  constituents  is  needed  (like  the  MSIS  and  Jacchia 
1977  models) .  #11  The  continuation  of  the  routine  middle  and 

orbital  atmospheric,  solar  and  geomagnetic  type  measurements  are 
still  needed  today  for  model  improvements.  In  summary:  The 
complexity  of  any  density  model  does  depend  on  the  models 
application  but  the  model  also  needs  to  be  relatively  simple  and 
usable.  Model  accuracy  needs  improvement,  and  more  measurements 
will  help. 

Since  the  November  1985  meeting  there  have  been  a  few  steps  taken 
toward  accomplishing  some  of  these  goals.  An  MSFC  clearing  house 
has  been  established  concerning  the  GRAM  model,  which  improves 
communications.  Likewise,  NASA  and  the  USAF  are  more  aware  of 
what  each  is  doing  so  that  duplication  of  effort  will  be  at  a 
minimum.  Use  of  the  recent  middle  atmospheric  data  sources  (i.e. 
satellite  radiometry,  limb  scanners,  lidar,  MST  radar,  shuttle 
drag  measurements)  are  currently  being  worked  into  the  GRAM 
atmosphere.  Funding  for  middle  and  orbital  atmospheric 
measurements/research  is  still  difficult  to  obtain  today.  NASA 
Headquarters  has  recently  reorganized  and  has  a  Space  Physics 
Division,  which  indicates  some  potential  for  orbital  studies. 

The  USAF/NASA  NASP  program  is  showing  some  receptiveness  toward 
stratospheric  and  middle  atmospheric  programs  which  would  aid 
NASP.  The  ENVIRONET  (NASA-GSFC)  currently  has  both  the  MSFC/TEAM 
and  the  MSIS-CIRA-86  models  for  its  users.  NASA/MSFC  is 
currently  evaluating  the  new  geomagnetic  storm  model  for  TEAM, 


54 


t 


and  will  be  receiving  a  thermospheric  wave  model  product  from  SAO 
by  June.  An  abbreviated  TGCM  will  be  tested  against  the  TEAM, 
MSIS,  and  J77.  The  AFGL  Jac  70/Tidal  model  is  also  new,  and  will 
be  evaluated.  The  MSFC  orbital  mechanics  people  have  recently 
indicated  a  willingness  to  help  compare  density  models  with 
satellite  data.  The  NASP  Natural  Environment  Committee  on  May 
11-12,  1988  will  evaluate  some  density  measurement  programs  like 
the  mega  lidar,  and  theoretical  density  perturbation  studies, 
such  as  the  thermospheric  air  glow  data  program.  The  up-coming 
CEDAR  measurement  program  could  use  results  from  this  theoretical 
study,  for  evaluation  of  emission  data.  Shuttle  re-entry 
accelerometry  data  is  also  currently  being  backed  out  into 
density  values  along  a  160-60  km  altitude  trajectory.  The  data 
sample  will  also  aid  GRAM  in  the  middle  atmospheric  region  by 
providing  some  quasi-horizontal  gradient  measurements  of  density. 

Shown  on  Figure  4  is  a  schematic  of  the  past  history  of  the  MSFC 
Global  Reference  Atmosphere  Model  (GRAM)  and  what  future  middle 
atmosphere  and  thermospheric  plans  are  currently  being 
accomplished.  A  GRAM-89  model  is  foreseen. 

We  feel  that  since  orbital  density  models  may  be  difficult  to 
improve  much  beyond  the  current  -15%  uncertainty  limit,  that 
efforts  in  enhancing/increasing  the  computer  capability  and  more 
research  in  understanding  more  regarding  errors  in  the  drag 
coefficient,  will  also  help  this  USAF  operational  problem. 
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Lnvironet  uses  MSFC/TEAM  5  MSIS-CIRA-86 

MSFC  Thermospheric  Geomagnetic  Storm,  and  Wave  Model  Studies  (in  time) 
Use  of  Rees  abrcviated  TGCM  (using  coefficients) 

AFGL  .Tac70/Forbes  Tide  Model 

InterModel  Comparisons  and  Thermospheric  Theoretical  Airglow  Studies 
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The  Space  Environment  Services  Center  (SESC) ,  located  in  Boulder, 
Colorado,  is  operated  jointly  by  NOAA  and  the  Air  Weather  Service 
to  provide  space  environment  services  to  meet  national  needs.  The 
SESC  provides  forecasts,  data,  and  guidance  to  the  Department  of 
Defense  through  the  Global  Weather  Central  and  serves  the  civil¬ 
ian  community  directly.  Among  the  services  provided  are  fore¬ 
casts  and  current  values  of  indices  required  by  civilian  users  to 
run  operational  and  developmental  neutral  density  models.  Table  1 
summarizes  the  users  who  are  operating  neutral  density  models  and 
who  utilize  products  from  the  SESC. 
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Table  1.  Neutral  density  modelers  who  use  SESC  indices,  fore¬ 
casts  or  guioance 


User 

Type  Model 

Indices 

Appl ication 

Altitude 

Naval 

Astronautics 

Group 

Modi f ied 
Jacch i a 

10  cm 

KP 

Navigation 

satellite 

positioning 

1000  km 

Martin 

Marietta 

MS  IS 

l  o  cm 

AP 

Thermal 
loading 
and  payload 
fairing  release 
during  launches 

80  km 

NASA  Marshall 
Space  Flight 
Center 

Jacch ia 

10  cm 

Ap 

Projection  Various 

of  re-entry 
time  and 

attitude  control 
heights 

North 

American 
Rockwel 1- 
Johnson 

Space 

Center 

Mod i 1 ied 

J acch : a 

10  cm 

{p 

KP 

Shuttle 

Drag 

Coef  f i cient 
refinement 

105-270 

nm 

Lincoln 

Laboratories 

Var' ous 

1  0  cm 

tJ?d 

Development 
and  evaluation 
of  models  and 
inputs 

Various 

Notes : 

10  cm  is 

solar  radio 

flux  at  10  cm 

wavelength  (2800  mHz) . 

Ap  and  Kp  are  standard  geomagnetic  indices. 

TED  is  the  total  energy  deposited  in  the  upper  atmosphere 
by  precipitating  particles  as  measured  by  the  NOAA 
Satellite . 
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Paul  W.  Schumacher,  Jr. 
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NAVSPASUR,  Dahlgren , Vi rginia 
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Abstract 


The  main  features  of  NAVSPASUR 's  catalog  maintenance  procedures 
are  described.  It  is  noted  that  the  general  theory  used  to 
propagate  orbital  elements  forward  in  time  does  not  at  present 
include  drag  perturbations.  Rather,  observed  changes  in  mean 
motion  (inferred  by  differential  correction)  are  used  to  estimate 
mean  decay  rates  in  semimajor  axis  and  eccentricity .  While 
producing  the  correct  element  trends  for  predictions  over  2-3 
days,  this  extremely  simple  "reduced-order"  decay  model  has 
quantitative  deficiencies  which  can  lead  to  temporary  mis- 
identif ication  of  satellites.  Probably  the  best  long-term  remedy 
is  to  implement  a  new  general  theory  which  includes  drag 
perturbations.  However,  this  would  require  a  long  development 
and  test/evaluation  time.  The  present  discussion  examines  a 
near-term  alternative  which  can  be  pursued  with  modest  effort, 
namely,  refining  the  reduced-order  decay  model.  On  this  approach 
one  hopes  merely  to  obtain  more  accurate  mean  decay  rates  than 
are  now  used.  Some  general  analytical  results  are  reviewed  for 
the  case  of  a  locally  exponential  density  profile.  For  more 
complicated  profiles  one  can  consider  using  numerical  quadrature 
rather  than  literal  expansions,  even  for  catalog  maintenance 
work,  because  the  reduced-order  equations  are  fairly  simple.  It 
does,  of  course,  remain  to  be  seen  whether  simple  techniques  of 
this  type  can  be  pushed  far  enough  to  effect  any  noticeable 
improvement  in  day-to-day  operations. 
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Present  NAVSPASUR  Procedures 

The  catalog  maintenance  process  at  NAVSPASUR  can  be  described  very 
generally  as  follows.  Upon  receipt  of  observations  and 
identification  of  the  satellites,  differential  correction  (D/C)  is 
initiated  to  update  the  orbital  elements  to  a  new  epoch.  The 
parameters  used  in  the  D/C  are  nonsingular  combinations  of  the 
classical  elements  designed  to  avoid  geometrical  indeterminancies 
at  low  eccentricity  and  low  inclination.  These  parameters  are 
akin  to  the  so-called  "equinoctial  elements",  but  have  the 
additional  property  of  being  valid  even  for  retrograde  equatorial 
orbits.  The  D/C  also  produces  values  for  the  parameters  M2  and 
M3,  which  are  essentially  rate  of  change  and  acceleration  of 
mean  motion  respectively.  These  extra  parameters  are  quite 
important  at  NAVSPASUR  since  they  are  the  primary  means  by  which 
orbital  decay  is  detected  and  accounted  for.  Indeed,  the  D/C 
process  itself  is  affected  critically  by  these  parameters.  3y  an 
empirically  derived  formula,  the  values  of  M2  and  M3  are  used 
to  determine  the  time-span  of  data  over  which  the  next  day's  D/C 
will  be  performed  for  that  satellite.  This  so-called  "fit  span" 
can  fall  between  3  and  30  days,  the  shorter  spans  corresponding  to 
higher  decay  rates-  For  exceptional  satellites,  operations 
analysts  can  manually  specify  fit  spans  outside  this  range.  The 
importance  of  selecting  a  suitable  fit  span  should  be  appreciated 
in  light  of  the  fact  that  drag  per turbat ions  are  not  well  modelled 
at  NAVSPASUR.  Without  accurate  drag  modelling,  observations 
remote  in  time  from  the  present  are  not  significant  for  near-term 
predictions  and  should  be  ignored.  The  empirical  fit  span  formula 
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provides  this  adjustment  in  a  systematic  way  tuned  to  the 
characteristics  and  requirements  of  the  NAVSPASUR  system. 

Once  the  orbit  has  been  differentially  corrected  and  a  new  epoch 
established,  the  elements  are  propagated  forward  in  time  using  the 
general  theory  of  Brouwer  [5]  as  modified  by  Lyddane  [7]  to 
handle  low-eccentricity  cases.  This  theory  is  notable  for  its 
accuracy  in  following  gravitationally  perturbed  orbits,  and  has 
been  implemented  in  a  concise  and  efficient  FORTRAN  program  at 
NAVSPASUR.  However,  the  theory  does  not  include  any  atmospheric 
perturbations . 

With  Brouwer  elements  available  at  a  given  time,  conversion  is 
made  to  the  classical  elements  which  in  turn  are  used  for  a 
variety  of  purposes,  including  position  predictions  to  be  compared 
with  future  observations.  Typically,  predictions  are  needed  for 
2  to  3  days  in  the  future,  long  enough  that  decay  effects  must 
somehow  be  included  even  though  the  general  theory  does  not 
include  them.  The  decay  corrections  now  being  used  are  applied  to 
the  classical  elements  themselves.  In  particular,  mean  motion  is 
augmented  with  a  term  in  M2-  The  quantity  M3  is  available  and 
could  be  included  in  the  mean-anomaly  polynomial,  but  it  has  been 
found  that  generally  better  predictions  ensue  when  M3  is  not 
included,  at  least  over  2  to  3  days.  This  effect  has  not  been 
scrutinized;  it  seems  to  be  due  to  better  smoothing  of  unmodelled 
effects.  In  virtually  all  cases  M3  is  very  small.  The  observed 
mean  motion  rate  M2  is  also  used  to  estimate  mean  decay  rates  in 
semimajor  axis  a  and  eccentricity  e. 
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The  term  a  is  obtained  by  di f Lerent iating  Kepler's  Third  Law  in 
the  form  n^a^=  constant.  A  minor  refinement  is  to  use  the 
Brouwer  "mean"  mean  motion  M ^  rather  than  the  osculating  mean 
motion  n  in  the  calcuation  of  a.  Subsequently,  the  mean 
eccentricity  rate  e  is  computed  in  terms  of  a  by  a  formula  to  be 
explained  later.  These  equations  are  examples  of  what  might  be 
called  a  "reduced-order"  decay  model:  less  than  the  full  set  of 
elements  is  included  in  the  drag  calculation,  based  (hopefully)  on 
reasonable  approximations. 

The  association  of  a  single-station  observation  with  a  prediction 
is  accepted  based  on  a  set  of  empirical  tolerances  (observed  minus 
expected)  in  the  values  of  time,  direction  cosine  and  Doppler 
shift  for  the  satellite.  Of  particular  interest  is  the  tolerance 
in  the  time,  which  can  vary  from  2  to  30  seconds.  This  time 
tolerance  is  an  empirical  function  of  semimajor  axis  a,  observed 
mean  motion  rate  M2  and  time  from  epoch  expressed  as  elapsed 
mean  anomaly  M-Mq .  Low,  rapidly  decaying  satellites  far  from 
the  epoch  call  for  the  largest  tolerance.  As  in  the  case  of  fit 
span,  the  importance  of  selecting  a  suitable  time  tolerance  for 
association  purposes  should  be  appreciated.  Since  drag  is  poorly 
modelled  in  the  system,  too  tight  a  tolerance  would  result  in 
failure  to  associate  the  observation  with  any  prediction.  On  the 
other  hand,  with  an  average  of  45,000  observations  being  collected 
over  a  24-nour  span,  it  is  likely  that  more  than  one  prediction 
could  be  associated  with  an  observation  when  the  time  tolerance  is 
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relaxed.  The  dilemma  is  unavoidable  as  long  as  decay  effects  are 
poorly  predicted.  The  empirical  answer  at  NAVSPASUR  has  been  to 
try  to  make  sure  that  at  least  one  association  can  always  be  made. 
In  case  of  more  than  one  possible  association,  the  candidate  with 
the  smallest  time  difference  is  accepted. 

There  are  obvious  shortcomings  in  this  strategy.  A  typical 
difficulty  is  that  a  sequence  of  closely  spaced  satellites  may 
change  their  "parade  order"  under  drag  perturbations  of  the  mean 
motions.  This  results  in  either  unidentified  or  mi s ident i f ied 
observations.  The  former  call  for  some  manual  intervention  to 
identify  the  satellite;  the  latter  lead  to  manual  D./C  in  which  the 
analyst  attempts  to  find  and  omit  the  offending  observations  from 
the  input  data.  Some  idea  of  the  magnitude  of  this  effort  can  be 
gained  by  noting  that  from  15  February  to  15  March  1988  about  1500 
manual  identifications  and  about  4500  manual  D/C's  were  required. 
Of  course,  these  figures  count  occasions  due  to  all  causes,  not 
just  those  due  to  decay  uncertainties.  Another  effect  of  poor 
drag  modelling  is  that  the  point  of  breakup  in  a  separation  event 
is  poorly  determined  due  to  differential  decay  of  the  pieces. 

This  sometimes  makes  it  impossible  to  choose  with  reasonable 
certainty  among  many  candidate  "blast  points".  One  might  well 
anticipate  an  increasing  future  demand  for  precise  locations  of 
breakups,  such  as  might  be  needed  to  distinguish  accidents  from 
deliberate  and  hostile  actions.  A  third  and  very  fundamental 
difficulty  caused  by  poor  decay  modelling  is  that  the  final  few 
orbits  in  the  satellite's  lifetime  are  poorly  followed.  Needless 
to  say,  this  final  phase  is  of  great  interest  for  many  satellites. 
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Improving  the  Drag  Calculations 

The  options  for  including  accurate  drag  perturbations  in  NAVSPASUR 
calculations  are  just  beginning  to  be  examined  afresh.  In  the 
past,  when  only  some  hundreds  of  satellites  were  being  cataloged, 
there  was  no  operational  imperative  to  implement  elaborate  drag 
predictions.  Even  now,  with  more  than  6500  objects  in  the 
catalog,  there  is  still  no  operational  imperative  to  include  drag 
more  accurately;  that  is,  NAVSPASUR  is  presently  quite  successful 
in  its  operations.  However,  the  analyst  workload  due  to  decay- 
related  anomalies  is  significant.  Moreover,  that  workload  is  an 
obvious,  though  as  yet  undocumented,  function  of  atmospheric 
activity  as  measured  by  geomagnetic  and  solar  flux  (10. 7cm) 
indices.  (These  indices  are  available  daily  at  NAVSPASUR.)  The 
nature  of  the  satellite  identification  problem  is  such  that  it 
must  get  worse  geometrically  as  the  number  of  low-altitude 
satellites  increases.  Hence,  it  is  only  prudent  to  explore 
possible  improvements  before  the  matter  does  become  imperative. 

In  the  event  that  high  accuracy  in  position  were  required,  say,  on 
the  order  of  10  meters,  NAVSPASUR  would  have  to  resort  to  a 
special  theory  and  simply  bear  the  computational  burden.  The 
relatively  few  requests  to  date  for  arcs  of  such  accuracy  have 
always  involved  rather  high-altitude  satellites  for  which  drag 
perturbations  are  either  insignificant  or  not  the  limiting  factor 
in  the  overall  accuracy.  Thus,  no  really  accurate  atmospheric 
model  has  ever  been  .needed.  But  the  choice  of  atmospheric  model 
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would  instantly  become  a  crucial  issue  if  high-accuracy  low- 
altitude  arcs  were  requested.  With  the  choice  made,  though, 
including  drag  effects  in  the  special  theory  is  straightforward. 
For  catalog  maintenance,  one  can  usually  be  satisfied  with 
position  accuracy  on  the  order  of  1  kilometer,  for  which  a  general 
theory  will  suffice.  This  is  a  fortunate  circumstance;  numerical 
integration,  even  at  this  accuracy,  would  be  prohibitively  time- 
consuming  for  catalog  maintenance.  Hence,  the  best  prospect  for 
improving  NAVSPASUR  catalog  maintenance  in  the  long  run  is  to 
implement  a  general  theory  which  accounts  for  drag.  The  theory  of 
Lane  and  Cranford  [4],  currently  used  at  NORAD/SSC,  is  certainly 
one  candidate.  The  theory  of  Brouwer  and  Hori  [1,6]  was  studied 
at  NAVSPASUR  more  than  twenty  years  ago,  but  finally  was  not 
developed  for  implementation  due  to  lack  (at  that  time)  of  a 
compelling  need  for  it.  Whatever  general  theory  is  selected,  one 
can  anticipate  a  long  development  and  test/evaluation  time  for 
such  a  critical  software  component.  If  near-term  results  are 
desired  and  only  modest  effort  can  be  allocatd  to  the  task,  there 
is  still  an  alternative  worth  exploring,  namely,  to  refine  the 
reduced-order  decay  model. 

Drag  models  consisting  of  coupled  equations  for  the  elements  a  and 
e  have  been  developed  extensively  by  Sterne  [3],  King-Hele  [2]  ar 
others.  Typically,  eccentric  anomaly  is  used  as  the  independent 
variable,  and  some  averaging  operation  is  used  to  obtain  mean 
rates  for  these  two  elements.  Considering  only  a  and  e  is 
equivalent  to  assuming  spherical  symmetry  of  the  atmosphere  and 
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neglecting  atmospheric  rotation.  Under  these  approximations,  the 
orientation  elements  are  not  affected  by  drag,  giving  rise  to  the 
simple  reduced-order  system  of  two  equations.  (It  will  be  assumed 
in  this  discussion  that  the  perturbation  of  mean  motion  is  an 
observed  quantity  and  is  available  for  input  to  the  model.) 

Usually  tne  two  equations  are  solved  to  first-order.  Then  at  the 
second  order  changes  in  the  orientation  elements  can  be 
considered,  which  arise  due  to  atmospheric  rotation,  atmospheric 
oblateness  and  coupling  with  gravitationally  -  induced 
precessions  of  the  orbit. 

It  may  be  found  that  numerical  quadrature  is  a  feasible  way  to 
obtain  the  first-order  change  in  a  and  e,  even  for  catalog 
maintenance  applications.  The  equations  are  fairly  simple  and  are 
valid  for  arbitrary  density/height  profiles.  This  approach  does 
require  a-priori  knowledge  of  the  ballistic  coefficient  for  the 
satellite.  Useful  comments  are  found  in  the  book  by  Sterne  [3]. 

If  the  dens i ty/he ight  profile  is  assumed  to  be  exponential  in  the 
vicinity  of  some  reference  height,  then  analytical  results  ire 
possible.  The  rate  equations  for  a  and  e  are  found  to  contain 
cosine  eccentric  anomaly  both  in  an  exponential  function  and  in 
algebraic  factors.  If  the  algebraic  factors  are  developed  in 
Fourier  series  and  the  equations  are  averaged  over  one  revolution, 
the  solution  appears  as  an  infinite  series  of  modified  Bessel 
functions.  This  is  the  approach  taken  by  King-Hele  C2]  in  his 
analysis  of  satellite  lifetimes.  There  are  a  variety  of  ways  to 
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evaluate  the  Bessel  functions  for  both  small  and  larye  values  of 
the  argument.  Another  way  of  treating  the  rate  equations  is  to 
include  the  exponential  function  in  the  Fourier  expansion.  The 
Fourier  coefficients  will  now  be  more  complicated  due  to  expanding 
the  exponential  function,  but  when  each  equation  is  averaged  over 
one  revolution  only  the  leading  coefficient  remains.  It  turns  out 
that  the  leading  coefficient  can  be  obtained  explicit'/  without 
developing  the  rest  of  the  expansion,  and  it  can  be  arranged  as  an 
infinite  series  in  the  eccent r ic i ty -  The  series  converges,  but 
whether  it  converges  fast  enough  to  be  useful  is  not  known  yet.  A 
little  inspection  of  this  solution  shows  that  the  formula  for  mean 
decay  rate  of  eccentricity  currently  used  at  NAVSPASLJR  is  really  a 
crudely  truncated  version  of  the  Fourier  series.  Evidently,  here 
is  the  place  to  begin  refining  the  current  model. 

When  fust-order  drag  variations  in  a  and  e  are  available,  one  can 
address  the  second-order  effects  mentioned  before.  Of  these,  the 
easiest  to  illustrate  is  the  coupling  with  gravitationally-induced 
mean  rates  in  the  node  and  perigee  angles.  The  latter  happen  to 
be  simple  functions  of  a  and  e  so  the  necessary  partial 
derivatives  are  easy  to  obtain.  Variations  in  node,  perigee  and 
inclination  due  to  rotation  of  the  atmosphere  and  atmospheric 
oblateness  are  more  complicated.  King-Hele  [2]  provides  extensive 
d i scuss ion . 
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There  is  an  altogether  different  type  of  analysis  which  may  lead 
to  a  useful  reduced-order  model  and  which  has  not  been  thoroughly 
investigated.  The  concept  is  presented  here  in  its  "half-baked" 
form  because  several  of  its  features  have  not  been  discussed  in 
the  literature.  The  analysis  begins  with  a  set  of  regularized 
equations  of  a  type  first  proposed  by  Burdet  [8].  The  independent 
variable  is  true  anomaly  elapsed  from  the  epoch,  and  the  dependent 
variables  consist  of  reciprocal  radius,  angular  momentum,  time, 
and  a  set  of  orthonormal  unit  vectors  defining  the  osculating 
plane.  In  terms  of  these  variables,  the  Cartesian  position  and 
velocity  vectors  have  simple  forms.  The  differential  equations 
have  the  interesting  property  that  they  are  rigorously  linear  for 
Keolerian  motion.  A  further  remarkable  property  appears  when  drag 
is  modelled  as  any  scalar  function  whatsoever  multiplying  the 
velocity  vector.  In  that  case  all  the  equations  except  the  one 
for  angular  momentum  revert  to  their  Kepler ian  forms-  In 
principle,  the  only  approximation  enters  in  the  solution  for 
angular  momentum.  The  rest  of  the  system  can  be  solved  exactly 
using  quadratures  at  most.  Mittleman  and  Jezewski  [9]  examined 
the  special,  somewhat  unrealistic,  case  in  which  angular  momentum 
can  also  be  obtained  exactly.  However,  they  did  not  discuss  the 
more  general  property  just  noted. 


* 


Other  Poss ibi 1 i t ies  Related  to  Drag  Modelling 

It  is  worth  mentioning  a  few  topics  that  could  be  addressed  if 
drag  perturbations  were  being  accounted  for  in  a  routine  and 
accurate  manner.  As  researchers  attempt  to  model  atmospheric 
phenomena  more  accurately,  there  is  always  the  problem  of 
observing  and  estimating  the  new  parameters  that  arise.  Satellite 
observations  have  been  a  prime  source  of  observations  for 
ionospheric  and  exospheric  phenomena.  Often  the  orbital  effects 
of  the  phenomena  would  be  detectable  even  though  they  are  not 
significant  for  satellite  operations.  In  this  case,  the 
possibility  exists  of  using  NAVSPASiJR  historical  data  to  validate 
density  models  over  extended  past  time.  Complete  observational 
archives  have  been  maintained  on  magnetic  tape  since  the  inception 
of  Naval  Space  Survei  Llance.  The  data  is  in  the  form  of  element 
sots  and  triangulated  positions  for  every  satellite  observed  by 
"JAVSPASUR.  Records  extend  over  more  than  two  complete  solar 
cycles  -  almost  thirty  years.  Although  NAV3PASUR  is  not  in  a 
position  to  participate  directly  in  such  work,  in  principle  the 
archives  could  be  made  available  to  researchers. 
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Another  possiblity  of  perennial  interest  is  to  estimate  the 
ballistic  coefficient  of  the  satellite  and  from  that  infer  values 
for  the  satellite  mass  or  cross-sectional  area.  In  using  the 
reduced-order  (a,e)  model,  the  ballistic  coefficient  would  be 
obtained  by  observing  the  decay  in  mean  motion  and  inserting  this 
value  into  the  differential  of  Kepler's  Third  Law,  producing  the 
decay  in  semimajor  axis.  Then,  given  a  suitable  density  profile, 
the  ballistic  coefficient  in  the  reduced-order  model  is  adjusted 
to  match  the  calculated  decay  in  semimajor  axis  to  its  known 
value.  Information  from  other  sources  would  then  be  needed  to 
isolate  any  of  the  three  parameters  comprising  the  ballistic 
coefficient. 

Another  problem  which  has  received  scant  attention  is  the  rigorous 
derivation  of  the  law  governing  the  drag  force  magnitude.  The  law 
universally  assumed  nowadays  is  really  just  the  definition  of  drag 
coefficient.  It  has  little  physical  basis  as  a  force  law.  It  was 
originally  derived  by  dimensional  analysis  of  the  variables  judged 
to  be  significant  in  the  inviscid  flow  of  an  incompressible  fluid. 
This  is  hardly  the  flow  regime  encountered  by  a  satellite.  Of 
course,  a  proper  critique  of  the  drag  law  cannot  be  made  until 
drag  perturbations  are  being  calculated  very  accurately. 

Finally,  if  drag  perturbations  are  being  predicted  well,  one  may 
hope  to  detect  smaller  effects  now  considered  intractable  to 
model.  Examples  might  be  differential  rotation  between  the  upper 
atmosphere  and  the  solid  Earth,  and  the  orbital  precession  due  to 
lift  of  a  large-area  satellite  which  has  been  attitude-stablized 
in  the  direction  of  the  sun.  Naturally,  it  is  doubtful  whether 
■mole  red  a  d~order  dr. an  mode  '  .  would  .  1  ■  w  one  to  n\.  ?.  rate 
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NAVSPASUR  ORBITAL  CALCULATIONS 
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SATELLITE  IDENTIFICATION 
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TOTAL  OBSERVATIONS  PER  DAY:  80,000*  RECEIVED 

45,000  PROCESSED 


TYPICAL  DIFFICULTIES  WHEN  DRAG  IS  SIGNIFICANT 
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OPTIONS  FOR  INCLUDING  DRAG  PERTURBATIONS  AT  NAVSPASUR 
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(3)  FOR  ASPADOC  SUPPORT 


TYPICAL  REDUCED-ORDER  DECAY  MODEL 


FOR  "LOCALLY  EXPONENTIAL"  ATMOSPHERE: 


SOLUTION  IN  MODIFIED  BESSEL  FUNCTIONS 

•  FIRST-ORDER  IN  (  Sfa  ) 

•  CHANGES  IN  ONE  REVOLUTION  (e=o  TO  e=z *) 
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RECURRENCE  FORMULAE 


SOLUTION  IN  FOURIER  COSINE  SERIES 

•  FIRST-ORDER  IN  (  &feF  ) 

♦  CHANGES  IN  ONE  REVOLUTION  (  E=<?  TO  E'iir) 
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MEAN  INCREMENTS  IN  NODE  AND  PERIGEE  RATES 
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OUGHT  TO  INCLUDE  ATMOSPHERIC  ROTATION 
AND  OBLATENESS  AT  THIS  ORDER 


A  POSSIBLE  REDUCED-ORDER  DECAY  MODEL  (BURDET  VARIABLES) 
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(MITTLEMAN  AND  JEZEWSKI,  1982) 


WHAT  IF  DRAG  IS  BEING  ACCOUNTED  FOR  "VERY  WELL"? 
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M.I.T.  Lincoln  Laboratory  Analysis  of  Satellite  Drag 


E.  M.  Gaposchkin  and  A.  J.  Coster 


ABSTRACT 


The  use  of  atmospheric  density  models  for  the 
calculation  of  the  drag  force  on  satellites  is 
evaluated.  First,  the  variety  of  different  roles 
that  drag  models  play  in  the  field  of  satellite 
tracking  and  space  surveillance  is  discussed.  Then, 
problems  associated  with  determining  the  atmospheric 
drag  on  a  satellite  are  presented.  This  is  followed 
by  a  description  of  the  current  suite  of  atmospheric 
models.  A  subset  of  these  models  is  selected  for 
further  quantitative  consideration.  This  subset 
includes:  Jacchia  1971  [CIRA  1972],  Jacchia  1977, 

Barlier  et  al.,  1977  [ DTM ] ,  and  Hedin  1983  [MSIS83]. 
These  models  are  evaluated  using  precision  tracking 
data  on  three  spherical  satellites.  The  satellites 
have  perigee  heights  of  270  km,  780  km,  and  1500  km. 
Cook's  1965  definition  of  Cjj  is  used  in  this 
analysis.  At  the  lower  altitude,  270  km,  all  models 
performed  equally  well.  At  the  higher  altitudes,  the 
models  do  not  perform  as  well.  Using  our  tracking 
data,  an  evaluation  is  also  made  of  the  atmospheric 
indices  currently  input  to  the  models,  the  Kp  and 
F10.7  cm  flux,  and  of  the  precipitation  index,  which 
is  not  yet  included  in  the  models.  Significant 
correlation  is  found  with  our  data  and  the 
precipitation  index.  Based  on  this  information,  a 
suggested  first  revision  to  a  version  of  the  Jacchia 
77  model  is  presented.  In  conclusion,  drag  data  can 
still  play  an  important  role  in  understanding  the 
thermosphere  and  in  contributing  unique  data  to 
monitor  this  next  solar  cycle. 


i .  INTRODUCTION 

Atmospheric  drag  affects  all  satellites  -  in  all  altitude 
regimes  -  from  low  altitudes  to  beyond  geosynchronous  altitudes. 
For  many  satellites,  the  modelling  of  atmospheric  drag  is  the 
largest  error  source  in  describing  the  force  on  the  satellite. 

The  focus  of  this  paper  is  the  use  of  atmospheric  density  models 
for  calculation  of  the  drag  force  on  satellites.  Atmospheric 
drag  is  viewed  from  the  vantage  point  of  satellite  tracking.  The 
data  analysis  presented  in  this  paper  comes  from  precision  radar 
tracking  data  beginning  in  1985  on  two  spherical  satellites  and 
from  laser  ranging  data  taken  in  1986  on  a  third  spherical 
satellite.  The  satellites  have  perigee  heights  of  270  km,  780 
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km,  and  1500  km.  These  data  were  used  to  evaluate  the  f c 1 . ew . n  i 
atmospheric  density  models:  Jacchia  1971  [C1RA  15  7  2  dec  -hi  a 
1977,  Barlier  et  al.,  1977  [ DTM ] ,  and  Hedin  1983  iMSIS83],  using 
Cook's  1965  definition  of  C^. 

As  introduction,  Figure  1  illustrates  several  aspects  or  the 
atmospheric  drag  problem.  First,  observe  the  drag  force 
equation.  The  drag  force  is  a  product  of  four  factors:  C'd  the 
ballistic  coefficient;  A/M,  the  a rea-to-mass  ratio;  p  ,  the 
atmospheric  density,  and  V  ,  the  speed  of  the  satellite  with 
respect  to  the  atmosphere.5  Each  of  these  quantities  will  be 
discussed  in  a  later  section  of  this  paper.  However,  note  nere 
that  none  of  these  quantities  is  known  precisely.  An  error  made 
in  the  calculation  of  atmospheric  drag  can  be  related  to  an  error 
made  in  determining  any  one  of  the  above  factors,  and  is  more 
likely  due  to  a  combination  of  errors  in  all  of  them. 

A  second  aspect  of  the  atmospheric  drag  problem  is  suggested 
by  the  two  separate  terms  for  density  listed  directly  below  the 
drag  force  equation.  Two  terms  are  used  because  none  of  the 
standard  atmospheric  models  predict  reliable  densities  above  27Cu 
km.  This  is  despite  the  fact  that  drag  effects  are  observed 
above  this  altitude.  The  reason  is  that,  until  recently,  the 
primary  interest  in  drag  has  teen  in  the  low  altitude  regimes. 

The  majority  of  data  used  in  building  the  atmospheric  models 
comes  from  regions  below  1000  km.  Additional  data  is  neeued 
above  1000  km.  Also,  it  is  probab’e  that  the  physical  process  of 
the  two  regions  (above  and  below  2000  km,  are  different,.  Tn  cur 
analysis,  we  use  an  integrated  atmospheric  density  model.  Below 
2000  km,  the  atmospheric  density  is  determined  from  one  of  the 
conventional  atmospheric  models  using  ir.pm such  as  the  '-'10.7  sr 
flux  and  the  geomagnetic  index,  P’p.  It  1 1so  assumed  runt  tno 
lower  atmosphere  is  corrotating  with  the  solid  earth.  Above 
km,  an  empirically  determined  density  1 1  wed  - n  inertial  r pace  is 
used . 

Finally,  in  Figure  1,  the  size  of  the  drag  ef'ect  is 
suggested  in  a  simplified  model  calculation.  Tn  ' his  fiqu  e  all 
satellites  are  assumed  to  be  spheres  with  th*-1  same  urea-tc-mass 
(A/M)  ratio.  The  A/M  is  chosen  to  be  0.1  cm-'/gram  which  is 
typical  for  satellite  payloads.  We  calculate  the  time  for  drv.g 
to  change  the  satellite  position  by  12  km.  This  is,  in  one 
sense,  the  orbit  error  committed  ignoring  drag  altogether.  It  3  ■ 
found  that  for  satellite  COS  11796,  which  is  ..  calibration  sphere 
with  perigee  at  about  300  km,  there  will  12  km  in  error  (along 
track)  in  0.92  days  or  13.6  revolutions.  For  LOS  4,  another 
calibration  sphere  with  perigee  at  about  S00  km,  there  v ; 1 1  be  a 
12  km  error  in  22.8  days  or  323  revolutions.  While  for  L.CS  1, 
calibration  sphere  with  perigee  at  23Cb  km  there  will  be  a  12  km 
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error  in  38.9  days  or  386  revolutions.  Clearly,  although  the 
drag  effect  is  far  more  noticeable  at  lower  altitudes,  there  is 
significant  drag  observed  at  all  altitudes. 

II.  USES  OF  DRAG  MODELS 

The  primary  uses  of  drag  models  are:  precision  orbit 
determination,  mass  determination  or  weighing  of  satellites,  and 
investigation  of  geophysical  phenomena.  The  first  two  uses, 
precision  orbit  determination  and  weighing  of  satellites,  are 
considered  time  critical.  Results  must  be  available  within  hours 
to  be  of  use.  The  last  one,  geophysical  investigations,  which 
includes  the  specific  area  of  atmospheric  physics,  is  less  time 
critical.  One  can  await  a_posteriori  analysis  to  obtain  an 
optimal  estimate  of  the  orbit.  Table  1  summarizes  the  uses  of 
atmospheric  drag  models  and  the  different  areas  in  satellite 
tracking  that  are  involved  in  drag  measurement. 

Table  1.  Uses  of  Drag  Models 

Precision  Orbit  Determination  (Time  Critical) 

Catalog  Maintenance 
SOI 

Predict ion/Forecasting 
ASAT  Targeting/Threat  Analysis 
Data  Screening  and  Calibration 
Collision  Avoidance 

Navigation  of  Satellites  (e.g.  Transit) 

Weighing  of  Satellites  (Time  Critical) 

SOI 

Damage  Assessment 
Decoying  of  Space  Assets 
Space  Debris  Characterization 

Geophysical  Investigations 
Atmospheric  Physics 

Model  Development 
Calibration  of  Other  Sensors 
Synthesis  with  Other  Data 
Other  Investigations 

Polar  Motion  and  Earth  Rotation 
Scientific  Satellites  Precision  Orbits 
(MAGSAT , SEAS AT , TOPEX) 

A.  Precision  Orbit  Determination 

The  main  use  of  drag  models  is  the  determination  of  precision 
orbits.  Because  atmospheric  drag  is  one  of  the  largest  forces  on 
a  satellite,  accurate  modelling  of  the  atmosphere  is  required  to 
obtain  precision  orbits.  The  ability  to  know  precise  satellite 
positions  contributes  to  a  variety  of  areas  in  satellite 
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tracking.  For  example,  precision  orbits  on  satellites  are  needed 
for  catalog  maintenance,  SOI  (Satellite  Orbit  Identification) , 
collision  avoidance,  satellite  navigation  (e.g.  Transit),  ASAT 
targeting  and  threat  analysis,  and  data  screening  and 
calibration. 

Catalog  maintenance  requires  knowledge  of  a  satellite's 
orbital  elements  within  a  given  accuracy.  For  data  of  a  given 
quality,  the  best  dynamical  model  of  the  orbit,  one  that  includes 
the  best  model  of  atmospheric  drag,  will  give  the  best,  or  most 
accurate,  element  set  for  a  satellite.  By  maintaining  an 
accurate  catalog,  one  has  available  precision  element  sets  for 
each  satellite.  A  critical  use  of  these  element  sets  is 
Satellite  Orbit  Identification  (SOI) .  Reliable  and  accurate 
orbits  provide  the  most  commonly  used  technique  for  the 
identification  of  satellites. 

Another  important  use  of  precision  elements  is  for  the 
prediction  and  forecasting  of  satellite  orbits.  Prediction  is 
needed  for  reasons  such  as  acquiring  more  tracking  data  from  a 
pencil  beam  radar,  or  determining  when  an  orbit  will  decay,  a 
requirement  which  stresses  the  atmospheric  model  at  low 
altitudes.  Predictions  are  also  necessary  for  the  manouvre 
planning  and  execution  needed  for  orbit  maintenance.  Highly 
accurate  and  timely  predictions  of  satellite  position  are  also 
required  to  target  within  a  given  volume  (ASAT  targeting) ,  or  to 
assess  the  threat  to  an  asset  of  another  vehicle  (threat 
analysis) .  In  addition,  orbit  predictions  play  an  important  role 
in  the  testing  of  new  surveillance  systems.  New  systems,  such  as 
the  space-based  radar  and  the  visible  optical  satellite  tracking 
system,  both  currently  under  development,  will  need  real-time 
precision  predictions  of  position,  up  to  24  hours  in  advance,  to 
interpret  the  sensor  data  properly.  Finally,  accurate 
predictions  are  needed  by  operational  satellites.  For  example, 
the  Navy  Transit  Navigation  satellites  are  limited  by  the 
accuracy  of  drag  models.  At  this  time  there  are  several  Transit 
satellites  which  are  surface  force  compensated  and  drag  free. 
These  satellites  are  expensive.  Better  drag  models  could 
alleviate  the  drag  issues  in  a  fundamental  way. 

Atmospheric  models  also  play  an  important  role  in  the 
calibration  of  satellite  tracking  data.  One  must  be  certain  that 
the  data  has  been  validated,  screened,  and  calibrated  in  order  to 
have  confidence  in  the  inferences  made  based  on  tracking  data. 
This  is  an  ongoing  data  analysis  function.  Reliable  and  accurate 
orbit  computation  is  needed  to  calibrate  the  data.  Calibration 
computations  are  no  more  accurate  than  the  precision  of  the 
orbit,  which,  in  low  altitudes,  is  limited  by  drag. 

A  final  application  for  precision  orbits  is  needed  in  the 
area  of  collision  avoidance.  There  is  a  growing  population  of 
cataloged  objects  in  space,  almost  7000  in  low  altitude.  The 
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possibility  of  collision,  although  small,  is  nevertheless  real. 
Despite  the  small  probability,  the  cost  of  a  collision  of  one  of 
our  satellites,  or  even  worse  of  our  debris,  with  a  foreign  asset 
is  so  high  that  every  step  must  be  taken  to  eliminate  such  an 
occurrence.  The  precise  and  timely  monitoring,  i.e., 
maintenance  of  the  catalog,  is  the  only  way  to  accomplish  this. 

B.  Mass  Determination  of  Satellites 

The  effects  of  all  non-gravitational  forces,  e.g.,  drag  and 
radiation  pressure,  are  proportional  to  A/M.  This  fact  is  often 
used  to  measure  A/M  from  the  observed  change  in  the  satellite 
orbit.  Given  that  information  is  known  about  the  size  (A)  of  a 
satellite  from  radar  cross  section  (RCS)  observations  then,  by 
measuring  the  drag  force,  the  mass  (M)  of  a  satellite  can  be 
determined.  This  has  been  done,  in  fact,  with  more  or  less 
accuracy  for  some  years.  It  hinges  on  knowing  the  density  and 
the  contribution  to  the  drag  force  equation.  Hence  it  places 
a  requirement  on  the  accurate  assessment  of  the  atmospheric 
density. 

The  determination  of  a  satellite's  mass  and  size  has 
important  implications  in  the  area  of  Satellite  Orbit 
Identification  (SOI) .  The  combination  of  information  concerning 
a  satellite’s  mass  and  size  (and  perhaps  its  shape,  attitude,  and 
spin  rate)  is  a  second  discrimination  tool  to  be  used  in  SOI, 
perhaps  equally  powerful  to  the  information  contained  in  the 
element  sets.  By  monitoring  this  information,  it  can  also  be 
used  to  indicate  if  a  change  has  taken  place  in  either  the  orbit 
or  in  the  operational  status  of  the  satellite.  Such  changes  can 
be  indicators  of  possible  damage,  in  the  event  of  an  ASAT  attack 
or  collision. 

Another  use  of  weighing  satellites  is  in  the  realm  of 
decoying  space  assets.  It  is  usual  to  view  the  decoy  problem  in 
terms  of  replicating  the  size,  RCS,  and  optical  characteristics 
of  the  asset  in  a  decoy.  However,  with  the  ability  to  weigh  the 
objects,  the  mass  will  also  have  to  be  replicated  in  order  to  be 
credible.  Since  the  cost  of  space  objects  depends  on  the  mass  in 
orbit,  such  a  requirement  sharply  reduces  the  attractiveness  of 
decoying.  Of  course,  issues  such  as  response  time  and  accuracy 
become  paramount,  since  decoy  strategies  are  manifold. 

Finally,  a  further  use  of  weighing  satellites  is  in  the  low 
altitude  regime  where  there  is  considerable  debris.  The  debris 
population  is  continually  changing.  Debris  is  depleted  by 
atmospheric  drag  and  replenished  by  a  variety  of  sources, 
including  breakups  and  new  rocket  launches.  An  organized 
observation  and  analysis  program  is  required  to  properly  assess 
the  hazard  and  the  statistics  of  the  population.  A  critical 
element  of  space  debris  characterization  will  be  the  measurement 
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of  each  particle's  mass,  which  can  be  done  with  a  good  drag  model 
and  accurate  tracking  data. 

C.  Geophysical  Investigations 

Drag  models  used  for  operational  tracking  can  provide 
significant  insight  into  the  basic  physics  of  the  Thermosphere. 
Excellent  models  have  been  derived  from  analysis  of  tracking 
data.  Such  tracking  data  can  iso  provide  basic  information  to 
test  models.  For  example,  drag  data  can  be  used  to  assess  the 
models'  performance  and  to  determine  basic  constants  within 
models.  Drag  data  can  also  be  used  to  calibrate  other 
atmospheric  sensors,  such  as  instrumented  satellites.  Over  the 
long  term,  one  can  hope  to  derive  a  complete  thermospheric  model 
based  on  physical  principles,  that  combines  observations  of 
constituents  from  satellite  instruments,  with  total  density 
obtained  from  analysis  of  satellite  drag. 

Other  geophysical  investigations  also  depend  on  precision 
orbit  computation.  Currently,  there  are  low  altitude  satellites 
used  for  measuring  the  earth's  geopotenr. iai ,  its  motion  in  mace, 
its  polar  motion,  and  its  rotation  rate.  Atmospheric  drag  is  a 
significant  error  source  in  the  analysis  of  these  low  altitude 
satellites.  In  addition,  there  are  a  number  of  low  altitude 
geophysical  sensors  (e.g,  MAG SAT,  SEASAT,  GEOSAT,  TOPEX)  that 
need  to  know  the  position  of  the  satellite  at  the  time  of 
measurement.  Again,  improvement  in  d-nag  modelling  will  lead 
directly  to  improvement  in  the  analysis  of  this  type  of  sensor 
data . 

III.  Outline  of  Drag  Problems 

The  primary  topic  of  this  conference  is  Neutral  Atmospheric 
Density  Models.  However,  this  emphasis  is  somewhat  misplaced  as 
the  main  interest  is  the  calculation  of  satellite  drag.  T’ 
calculation  of  neutral  density  is  only  one  piece  of  the  pr' i  vem. 

A  number  of  other  issues  have  to  be  addressed  in  order  to  acr.ieve 
the  desired  capabilities. 

As  shown,  the  drag  force  per  unit  mass  on  a  satellite,  which 
is  the  force  that  is  measured,  is  determined  from 

F  =  (l/2)Cd(A/M)PVs2  (1) 

where  Cd  is  the  ballistic  coefficient,  A/M  the  area-to-mass 
ratio,  p  is  the  atmospheric  density,  and  Vs  is  the  speed  of  the 
satellite  with  respect  to  the  atmosphere.  Vs  can  be  written  as 
the  speed  of  the  satellite,  Vsat,and  the  speed  of  the  atmosphere, 
Vatxns*  *n  our  analysis,  we  assumed  that  the  lower  atmosphere  co¬ 
rotates  with  the  earth. 

One  of  the  main  problems  in  determining  che  drag  on  a 
satellite  is  that  none  of  the  quantities  in  this  product  is  known 
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without  error.  For  example,  there  are  fundamental  unknowns  in 
the  physics  of  scattering.  This  issue  is  grouped  in  the 
definition  of  which  is  given  in  Figure  2.  has  to  be 

determined  for  each  satellite.  It  will  be  discussed  in  more 
detail  in  the  following  section. 

The  projected  cross-section,  A,  in  Equation  1  also  has  to  be 
determined  for  the  majority  of  satellites.  In  general,  the  cross 
section  depends  on  a  satellite's  aspect  angle.  This  is  often 
difficult  to  calculate  and  probably  will  have  to  be  observed.  In 
our  study  we  determined  the  atmospheric  drag  on  selected 
spherical  satellites  so  that  questions  of  aspect  did  not  arise. 

To  calculate  the  area-to-mass  ratio  (A/M)  in  Equation  1,  the  true 
mass  (M)  must  also  be  known.  This  is  essential  if  drag 
measurements  are  to  be  used  to  predict  absolute  density. 

The  next  term  in  Equation  1  is  P  ,  the  atmospheric  density. 

This  value  is  usually  predicted  by  the  atmospheric  models 
although  the  models  themselves  do  not  predict  the  density  to 
better  than  15%  at  low  altitudes  (less  than  300  km).  The  model 
predictions  get  even  worse  with  increasing  altitude. 

Finally,  the  last  term  in  Equation  1  is  the  speed  of  the 
satellite  with  respect  to  the  atmosphere.  The  speed  of  the 
satellite  with  respect  to  the  earth  is  known  with  fairly  high 
precision.  However,  this  term  also  depends  on  the  speed  of 
atmospheric  winds.  It  is  clear  from  satellite  experiments,  drag 
analysis,  and  theoretical  models  that  there  are  significant  winds 
and  gravity  waves  in  the  thermosphere.  These  winds  can  have  a 
speed  of  several  hundred  meters  per  second,  a  speed  which  is 
comparable  to  satellite  velocities.  Therefore,  any  advance  in 
the  modelling  of  satellite  drag  must  include  information  about 
winds . 

Other  problems  related  to  the  determination  of  atmo' pheric 
drag  include  the  abundant  evidence  that  sensible  drag  exists  at 
all  altitudes.  From  the  standpoint  of  space  surveillance,  high 
altitude  drag  is  not  as  serious  as  low  altitude  drag.  The 
cumulative  effects  take  months  to  be  operationally  significant. 
Yet  the  existence  of  this  drag,  not  predicted  by  any  of  the 
current  models,  indicates  a  flaw  in  our  understanding  of  the 
thermosphere.  We  cannot  accept  a  model  as  comprehensive  if  it 
does  not  satisfactorily  explain  all  the  observed  atmospheric 
phenomena. 

Further  difficulties  in  the  determination  of  atmospheric  drag 
can  be  found  in  the  geophysical  inputs  to  drag  models.  Recall 
that  in  the  previous  section  we  identified  a  number  of  important 
uses  that  were  time  critical,  i.e.,  that  need  results  within 
hours  to  be  useful.  Most  of  the  present  suite  of  models  use 
geophysical  parameters  such  as  F10.7  and  Kp  as  measures  of 
energy  input  to  the  atmosphere.  These  are  surrogate  parameters  to 
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begin  with.  The  reason  is  that  direct  measurements  of  the 
ionospheric  current  and  the  solar  EUV  flux  are  not  available. 

The  indices  that  are  used  have  certain  inherent  problems.  The  Kp 
and  Ap  are  planetary  indices  and  so  can  not  represent  localized 
disturbances.  The  F10.7  cm  flux,  on  the  other  hand,  is  used  to 
represent  the  EUV  solar  flux,  although  it  in  itself  has  no  direct 
influence  on  the  atmosphere.  The  models  also  assume  thac  we  use 
data  obtained  after  careful  calibration  and  reduction  by  the 
geophysical  agency,  NOAA/USAF.  Generally  the  final  data  values 
are  available  only  after  some  weeks  or  months,  and  are  not 
available  for  our  time  critical  mission.  We  must, 
therefore,  use  predictions  of  these  values  based  on  incomplete 
information.  A  certain  amount  of  error  is  associated  with  this 
approximation,  and  it  needs  to  be  quantified. 

Table  2  summarizes  the  problem  areas  associated  with 
determining  atmospheric  drag.  In  the  following  section  we  will 
explore  one  problem  area  in  more  detail. 


Table  2.  Summary  of  Proolems  Associated  with  Atmospheric  Drag 
Physics 

Satellite  Aspect 
Geophysical  Input  : 

Accuracy  and  Prediction 
Composition,  Temperature,  and  Density  Models 
Winds  and  Super-Rotation 
Gravity  Waves 
High  Altitude  Drag 


IV.  THEORY  OF  SCATTERING 

One  has  to  model  the  scattering  mechanism  to  determine  the 
ballistic  coefficient  C^  used  in  Equation  1.  However,  this 
process  is  not  well  understood.  For  example,  we  can  not  reliably 
predict  how  particles  in  free  molecular  flow  scatter  from  a 
surface  in  space,  nor  is  it  known  how  that  surface  and  the 
scattering  change  after  long  exposure  in  space. 

A  theory  of  scattering  is  summarized  in  Figure  2.  It  is  a 
theory  based  on  general  principles.  Verification  of  the  theory 
is  needed  along  with  definition  of  several  constants.  Figure  2 
gives  the  mathematical  formalism  for  solar  radiation  pressure, 
neutral  drag,  and  charge  drag.  These  three  effects  have  some 
common  elements.  It  is  instructive  to  consider  them  together. 


95 


In  all  three  cases,  the  force  on  a  general  surface  element 
depends  on  the  scattering  mechanism  and  the  flux,  the  flux  of 
neutral  particles,  photons,  or  charged  particles .  Two  scattering 
mechanisms  are  indicated,  specular  and  a  Lamberts  law  "diffuse" 
scattering.  Other  scattering  mechanisms  are  possible.  These  two 
are  considered  to  be  the  limiting  cases;  one  in  which  the 
scatterer  completely  "remembers"  the  information  about  the 
incoming  flux  vector_momentum  (specular) ,  and  the  other  in  which 
no  information  is  "  remembered"  (diffuse).  In  general,  surfaces 
do  not  behave  completely  as  specular  or  diffuse.  A  better  model 
would  be  to  assume  a  fractional  part  of  both  the  diffuse  and 
specular  components  on  the  satellite  surface. 

For  neutral  drag,  the  scattering  interaction  depends  on  the 
surface  material  of  the  satellite,  and  on  the  molecular  weight 
and  temperature,  or  thermal  velocity,  of  each  particle.  The 
thermal  dependence  is  theoretically  modelled  in  terms  of  the 
ballistic  coefficient.  The  molecular  weight  of  the  atmospheric 
constituent  and  that  of  the  surface  material  of  the  satellite  are 
modeled  through  an  accomodation  factor, 5  .  However,  even  if  the 
chemical  constituents  of  the  atmosphere  are  known  (using 
laboratory  measurements  from  Herraro,  1983),  the  value  of  can 
not  be  determined  to  better  than  5%  (Gaposchkin  &  Coster  1986) . 

V.  PRESENT  DRAG  MODELS 

Atmospheric  density  models  for  heights  above  120  km  have 
been  derived  from  the  analysis  of  satellite  drag  since  the 
launch  of  Sputnik  I.  Satellite  drag  measures  only  the  total 
density  and  contains  no  direct  information  about  composition. 
Early  models  identified  the  fundamental  dependence  of  the  upper 
atmosphere  on  solar  flux,  geomagnetic  index,  diurnal,  monthly  and 
seasonal  variations.  Atmospheric  models  based  on  satellite  dng 
data  are  typified  by  the  COSPAR  International  Reference 
Atmosphere  of  1972  (CIRA72),  which  is  fundamentally  based  on  the 
Jacchia  1971  model  (Jacchia,  1972),  and  by  the  DTM  1978  model 
(Barlier  et  al.,  1977). 

Atmospheric  composition  can  be  inferred  or  measured  using  both 
ground-based  incoherent  backscatter  radar  measurements  and 
satellites  instrumented  with  mass  spectrometers  and 
accelerometers.  This  type  of  data  has  been  used  to  construct  the 
so-called  Mass  Spectrometer  and  Incoherent  Backscatter  (MSIS) 
models,  typified  by  Hedin  et  al.,  1977  a,b,  1983. 

Jacchia  attempted  to  merge  drag  and  composition  data  into  a 
combined  model,  Jacchia  1977.  The  Jacchia  1977  model  has  two 
modifications:  a  1977  addendum,  and  a  revision  [Jacchia  & 

Slowly,  1982].  An  additional  version  of  the  Jacchia  1977  model 
known  as  the  Jacchia-Bass  model  [Bass,  1980  a,b]  was  developed  at 
the  Air  Force  Geophysical  Laboratory. 
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All  models  ?re  relatively  simple.  They  can  be  characterized 
as  static  diffusion  models  that  only  incorporate  dynamics 
implicitly.  The  processes  are,  as  yet,  too  complex  to  formulate 
a  model  based  on  physical  principles  alone. 

The  models  we  tested  inirially  were  the  Jacchia  1977  and  the 
MSIS  83  models.  We  also  tesred  the  Jacchia  1971  (CIRA72)  model 
since  it  is  a  'widely  used  standard,  and  the  DTM  1978  mode]  which 
was  designed  specifically  to  evaluate  satellite  drag.  However,, 
during  the  analysis,  a  number  of  issues  concerning  the  Jacchia 
model  arose  that  required  testing  of  several  variants.  Finally, 
the  results  on  seven  models  have  been  assembled.  The  development 
and  testing  of  these  models  led  to  additional  insights  into  the 
models  and  the  scattering  mechanisms  used  by  the  models. 

The  seven  models  are  listed  in  Table  4.  The  four  different 
versions  of  the  Jacchia  1977  model  evolved  because  of  two 
revisions  dealing  with  geomagnetic  effects.  The  first  was  an 
analytic  approximation  to  a  geomagnetic  effect  described  by  Eq. 

32  in  the  original  1977  Smithsonian  Astrophysical  Observatory 
Special  Report  No.  375  (Jacchia  1977) .  The  second  revision  was 
described  in  Jacchia  &  Slowey,  1981  resulted  from  combining  drag 
data  with  ESR04  satellite  data.  These  revisions  are  discussed  in 
more  detail  in  Gaposchkin  and  Coster  1986.  J771  includes  the 

analytic  approximation,  J77’’  includes  the  Jacchia  &  Slowey  1981 
revision,  and  J 77,,,  includes  both  sets  of  revisions. 

The  J71 ,  J77 ,  J77 ' ,  and  DTM  models  are  predominantly  based  on 
satellite  drag  measurements  obtained  in  the  region  from  250  to 
1000  km.  In  addition,  the  J77  •  '  and  J77'''  models  incorporate  a 
considerable  amount  of  satellite  measured  composition  data  from 
the  ESRO  4  satellite  at  altitudes  ranging  250  to  800  km  and  less 
than  40  degrees  in  latitude.  The  MSIS  model  is  completely  based 
on  satellite  mass  spectrometer  and  ground-based  incoherent 
scatter  data,  the  latter  used  primarily  to  measure  neutral 
temperatures.  The  majority  of  mass  spectrometer  data  was 
obtained  from  the  Atmospheric  Explorer  satellites  in  the  altitude 
regions  from  100  to  500  km. 


97 


» 


Table  4.  Seven  Density  Models  Tested 


J71 
DTM 
MSIS83 
J77 
J77  ' 
J77  '  ' 
J77  '  '  ' 


=  Jacchia  71,  aka  CIRA  1972 

=  Barlier  et  al.,  1978 

=  Hedin,  1983 

=  Jacchia  77,  as  defined  in  SAO  SR  375 
=  Jacchia  77  +  analytical  approximation 
=  Jacchia  77  +  1981  changes  in  AFGL  report 

=  Jacchia  77  +  1981  changes  +  analytical 

approx . 


A  timing  test  was  performed  on  the  different  models.  In  the 
test,  each  model  was  called  4000  times  in  a  variety  of  positions 
and  hour  angles.  The  results  are  presented  in  Table  5. 


Table  5.  Timing  Test  of  the  Atmospheric  Models 


J71 
DTM 
311 
311  ' 
311  '  ' 
311 '  '  ' 
MSIS83 


11.10  seconds 
16.22  seconds 
35.38  seconds 
41.98  seconds 
33.26  seconds 
36.70  seconds 
116.84  seconds  (86.74) 


The  two  times  listed  for  the  MSIS83  model  refer  to  the  how 
the  geomagnetic  Ap  parameters  were  input.  The  MSIS83  program  has 
the  option  of  using  either  the  daily  value  of  the  Ap  or  an  array 
of  7  Ap  indices,  including  the  daily  Ap  index  and  different  time 
averages  of  the  3-hr  Ap  indices.  The  shorter  time  listed  for 
MSIS83  refers  to  program  run  with  only  the  daily  Ap  value  input. 
No  significant  difference  was  found  between  these  two  options. 

In  our  standard  procedure,  the  array  of  Ap  values  as  inputs  to 
the  MSIS83  model  was  used. 


Our  timing  test  shows  that  the  MSIS  83  model  is  significantly 
slower  than  all  of  the  Jacchia  models.  This  result  is  in 
contradiction  with  the  general  consensus  in  the  community  which 
is  that  the  J77  model  is  much  slower  computationally  than  both 
the  J71  model  and  the  MSIS  model  (Liu,  et  al.,  1983,  Afonso,  et 
al.,  1983).  We  assume  this  is  because  our  implementation  of  the 
J77  model  used  a  look  up  table  similar  to  that  used  in  the  J71 
program.  Note  also  that  the  DTM  model  runs  faster  than  all  the 
models  except  for  the  J71  model. 

We  ran  a  computer  simulation  in  which  the  ratios  of  the 
densities  predicted  by  the  JACCHIA  71,  J77,  311' ,  311'''  ,  and 
MSIS  83  models  were  compared  to  those  predicted  by  the  J77 ' ' 
model  for  the  300  km  altitude.  This  simulation  was  done  to  see 
how  muc1-  variance  there  is  between  the  different  models.  The 
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J77 ' '  model  was  selected  as  the  standard  model  because  it  was 
found  to  be  the  best  overall  model  in  predicting  satellite  draq 
for  our  data  set.  All  hour  angles  were  sampled  with  latitude 
coverage  between  +60  to  -60  degrees.  The  daily  and  average 
solar  flux  values  were  set  at  74,  and  a  Julian  date  of  46144  was 
used.  The  results  for  two  different  Kp  values  corresponding  to  a 
typical  mean  Kp  value  (Kp  =  2.5)  and  to  a  high  Kp  value  (Kp  = 

5.0)  are  presented  in  Table  6.  In  each  case,  twenty-five  hundred 
data  values  were  averaged. 


Table  6.  Comparison  of  Densities  Derived  by  the  Various 
Atmospheric  Models  to  the  J77"  Model  at  300  Km 


KP  = 

2 . 5 

KP  = 

5.0 

Mean  Ratio 

Stnd.  Dev. 

Mean  Ratio 

Stnd.  Dev. 

(Model/J77 

") 

(%) 

(Model/J77" ) 

(%) 

1 . 

Jacchia 

71 

1.082 

6.5 

1.16 

9.1 

2  . 

Jacchia 

77 

J77 

0.886 

6.3 

0.81 

8 . 7 

J77  ' 

0.872 

7 . 3 

0.81 

9.4 

J77  " 

1.000 

0.0 

1 . 00 

0.0 

J77  '  '  ' 

0.880 

8.4 

0.78 

14 . 3 

3  . 

DTM  78 

0.867 

5.3 

0.81 

5.1 

4  . 

MS IS  83 

0.934 

9.3 

0.93 

10.2 

Clearly 

the 

results  show 

that 

at  300 

km  all  of  the 

models  are 

in  reasonable  agreement.  A  15%  difference  between  the  mean 
ratios  of  the  models  is  apparent.  It  is  known  that  atmospheric 
models  can  not  predict  densities  better  than  15%  (a  number 
which  increases  at  higher  altitudes).  Space  Command,  however, 
would  like  a  5%  prediction  capability  of  drag. 

In  Gaposchkin  &  Coster  (1986)  we  found  that  the  J77,  J77 '  and 
J77 '  '  '  were  clearly  inferior  to  the  J77 ' ’  in  all  cases.  These 
models  were  not  evaluated  further. 

VI.  MEASUREMENTS  OF  SATELLITE  DRAG 

The  rest  of  this  paper  is  concerned  with  measurements  of 
satellite  drag  and  with  what  can  be  done  with  these 
measurements.  In  particular,  we  describe  the  measurements  we 
made  of  atmospheric  drag  and  how  we  used  these  measurements  to 
evaluate  the  different  atmospheric  models.  We  were  also  able  to 
use  these  measurements  to  evaluate  different  atmospheric  indices, 
including  those  that  are  currently  used  as  inputs  and  a  new  index 
that  is  not  yet  being  used  by  the  models.  Following  this,  we 
will  discuss  a  suggested  first  revision  to  the  J77''  model  based 
on  a  least  squares  fit  of  our  measurements  to  the  following 
terms  :annual,  semi-annual,  seasonal-lat ifudinal ,  diurnal,  , 

PI .  p 


t  s 


Measurements  of  satellite  drag  were  acquired  by  taking  daily 
tracks  of  spherical  satellites  from  Lincoln  Laboratory  and  NASA 
facilities.  Two  satellite  tracking  radars  that  provide  data 
with  an  accuracy  of  1  meter  [Gaposchkin  1985,  1986  a,b]  are 
operated  by  Lincoln  Laboratory.  One  radar  is  the  Altair  radar 
located  on  the  Kwajalein  Atoll,  the  other  is  the  Millstone  L-band 
radar  located  in  Westford,  Mass.  A  network  of  laser  ranging 
stations  that  provide  data  on  satellites  equipped  with  cube 
corner  reflectors  is  operated  by  NASA.  The  accuracy  of  this  data 
is  better  than  2  cm  [Gaposchkin  et  al.,  1987]. 

Daily  tracks  of  data  on  two  satellites,  LCS4  and  COS1179, 
have  been  taken  since  the  beginning  of  1985.  Daily  tracks  on 
the  third  satellite,  EGS  aka  Agasii,  have  been  taken  since 
launch  in  1986.  Data  from  all  of  1985  have  been  analyzed  for  the 
former  satellites,  as  well  as  the  initial  two  months  of  data  on 
EGS  in  1986.  These  satellites  are  described  in  Table  7  where  the 
semi-major  axis  (a)  and  the  perigee  height  £q)  are  given  in  km, 
the  area-to-mass  ratio  (A/M)  is  given  in  cm  /gram,  and  the 
inclination  (I)  in  degrees.  The  eccentricity  is  referred  to  as 
e. 


Table  7.  Satellites  Used  for  Evaluation  of  Density  Models 


Cospar 

# 

NSSC# 

Name 

a  (km) 

e 

K°) 

q  (km) 

A/M 

1980 

37 

A 

11796 

COS1179 

7028 

0.051 

82.9 

270 

0.038 

1971 

67 

E 

5398 

LCS#4 

7207 

0.008 

87 . 6 

780 

0.285 

1986 

61 

A 

16908 

EGS 

7878 

0.001 

50.1 

1500 

0.045 

Atmospheric  density  models  were  used  in  the  orbit  computation 
of  these  data  sets  to  calculate  atmospheric  drag.  The  use  of 
these  models  included  the  final  data  on  the  10.7  cm  solar  flux 
and  the  geomagnetic  index  Kp  (Solar  Geophysical  Data,  1985, 

1986) .  A  scale  parameter  (S)  is  introduced  in  the  orbit 
determination  program.  S  is  a  least  squares  fit  parameter  that 
is  used  for  each  orbital  arc.  S  scales  the  entire  drag  model. 

It  can  be  interpreted  as  an  indication  of  the  adequacy  of  the 
drag,  and  thus  of  the  atmospheric  model  used  to  compute  the  drag. 
If  S  is  less/greater  than  unity,  then  the  atmospheric  density 
predicted  by  the  model  is  too  large/small. 

With  a  complete  force  model,  and  the  scale  factor,  the 
orbital  arcs  are  generally  fit  to  the  accuracy  of  the  data 
(several  meters) .  The  orbital  arcs  are  computed  with  one  day 
spacing,  using  between  two  and  four  days  of  data.  Any  given  pass 
of  data  will  be  in  at  least  two  orbit  fits.  This  is  used  for 
data  validation,  as  well  as  checking  orbital  consistency.  The 
computed  scale  factors  can  be  plotted  as  a  function  of  epoch. 
These  plots  can  then  be  used  to  give  an  indication  of  how  the 
density  model  involved  in  the  drag  calculation  may  be  deficient. 
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Figure  3  is  an  example  oi  the  scale  factor  data  1  or  satellite 
11796  and  for  the  J771'  model  aver  a  years  time. 


Scale  factors  were  computed  for  each  satellite  and  for  e^ch 
atmospheric  nsocel  .  In  addition,  scale  factors  were  computes  for 
both  the  specular  and  diffuse  cattering  case.  These  data  are 
summarized  in  Table  8.  The  first  term  under  each  satellite  is 
the  average  scale  factor,  S,  for  that  data  set.  The  second  term 
is  the  standard  deviation,  q  . 


Table  8.  Summary  of  Scale  Factors 


Diffuse 

Satellite 

Model 

J71 

J77" 

MSIS 

DTM 


11796 

s/a 

1. 106/0. 146 
1. 177/0. 138 


5398 

S/0 

0.715/0.321 

0.753/0.345 

0.878/0.255 


16908 

s/o 

1 . 293/0.691 
1. 389/0.617 
1.665/1.220 
1. 344/0.570 


Specular 
Satellite  11796 
Model  S/a 


5398 "C" 


5398 "G" 


16908 


s/a 


s/o 


s/a 


J71 

1.024/0.114 

1.024/0.374 

1.111/0.374 

311 '  ' 

1.108/0.153 

0.876/0.333 

0.950/0.333 

MSIS 

1.183/0.152 

0.923/0. 325 

1 . 001/0.325 

DTM 

1.152/0.169 

1.036/0.229 

1. 124/0.269 

1.919/0.839 
2 . 397/ 1.220 


The  results  for  each  satellite  will  now  be  discussed. 


COS 117 9 

COS 1179  (NORAD  Space  Surveillance  Center  [NSSC] #11796)  is 
known  to  be  a  sphere.  Its  mass,  however,  is  not  known  from 
independent  information.  Therefore,  an  area-to-mass  (A/M)  ratio 
based  on  an  average  drag  has  been  adopted.  Because  of  this  11"06 
can  only  be  used  for  a  relative  assessment  of  the  atmospheric 
models.  The  scale  factors  for  311''  were  computed  assuming  a 
diffuse  scattering  mechanism  and  are  presented  in  Figure  3. 

Throughout  the  time  period  analyzed,  11796  has  a  sufficiently 
large  eccentricity  that  drag  occurs  mainly  at  perigee,  wher^  the 
scale  height  is  about  7  km.  Therefore,  the  orbit  only  samples 
the  density  at  one  geographic  point  each  revolution.  Figure  4 
gives  the  latitude  and  solar  hour  angle  of  the  perigee  point  for 
the  interval  analyzed.  One  can  see  that  both  latitude  and  hour 


angle  are  completely  sampled,  and  there  is  no  simple  correlation 
with  the  variation  in  S  seen  in  Figure  3. 

Scale  factors  were  computed  again  for  this  satellite  using  the 
assumption  of  specular  scattering.  The  mean  values  of  the  scale 
factors  for  these  two  cases  are  fairly  similar  as  can  be  seen  in 
Table  8.  From  the  theoretical  standpoint,  a  significant 
difference  is  not  expected  between  the  specular  and  diffuse  cases 
because  the  drag  at  295  km  is  primarily  from  oxygen  and  nitrogen. 
The  mean  molecular  weight  is  approximately  16,  which  results  in  a 
value  for  the  accommodation  coefficient  of  nearly  1.  The 
accommodation  coefficient  is  the  term  used  in  the  model  for 
diffuse  scattering.  An  accomodation  coefficient  of  1  results  in 
the  same  value  of  as  that  predicted  by  the  model  for  specular 
scattering.  The  scale  factors  for  the  specular  scattering  are 
not  significantly  different. 

As  mentioned  above,  the  absolute  mean  scale  factors  are  not 
significant  in  this  satellite.  However,  the  fact  that  the  scale 
factors  differ  by  at  most  5%  indicates  that  all  the  models  are  in 
good  agreement.  The  variation  about  the  mean  is  smallest  for  the 
J77 1 '  model,  though  all  the  models  differ  by  2%  at  most.  On  this 
basis  one  could  marginally  chose  the  J77 ' '  as  the  best  model. 

This  difference  is  not  believed  to  be  significant.  All  models 
are  equally  good  (or  bad)  at  275  km  altitude  for  all  latitude, 
hour  angle,  and  geophysical  data. 

LCS4 

LCS4  is  a  satellite  built  by  Lincoln  Laboratory.  It  is  in  a 
circular  orbit  at  approximately  780  km  altitude.  Its  physical 
characteristics,  such  as  its  area-to-mass  ratio  (A/M)  and  the 
composition  of  its  outer  surface  (aluminum),  are  known.  Because 
of  its  nearly  circular  orbit,  the  satellite  does  not  provide  a 
clear  association  of  the  position  and  the  effect  of  dr.-  .. 
Nevertheless,  the  average  drag  can  be  discussed. 

At  this  altitude  the  scale  factors  for  all  models  are  quite 
different.  The  orbital  fits  have  been  carefully  scrutinized.  The 
observed  variability  in  S  is  believed  to  represent  real 
variations  in  the  atmospheric  density  at  780  km  altitude  that,  are 
not  predicted  by  any  of  the  models.  In  all  cases,  it  can  be  seen 
that  the  scale  factors  are  20%  to  30%  less  than  unity,  indicating 
that  the  predicted  drag  is  too  large  by  this  factor. 

At  780  km  altitude  the  principal  atmospheric  constituents  are 
hydrogen  and  helium.  For  these  constituents,  the  difference 
between  specular  and  diffuse  scattering  is  quite  large,  resulting 
in  diffuse  ballistic  coefficients  20%  to  30%  larger  than  specular 
ones.  The  mean  scale  factors,  assuming  the  specular  scattering 
mechanism,  are  approximately  0.9.  The  general  agreement  among 
the  various  models  persuades  us  that: 
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a)  the  specular  scattering  model  is  correct  for  the 
polished  aluminum  spherical  satellite,  and 

b)  there  may  be  a  systematic  overestimate  of  the 
density  by  as  much  as  12%. 

The  J71  model  has  the  best  mean  value,  and  the  DTM  has  the 
smallest  variability.  However,  the  differences  may  not  be 
significant  considering  that  the  uncertainty  in  may  be  as  much 
as  5%. 

A  further  issue  arises  due  to  an  uncertainty  about  LCS4 .  LCS4 
was  one  of  two  satellites  fabricated  by  Lincoln  Laboratory  at  the 
same  time.  Both  were  one  square  meter  radar  calibration  spheres 
made  of  polished  aluminum,  identical  in  all  but  one  respect.  The 
records  show  that  one  sphere,  the  "GH  sphere,  failed  on  launch 
and  never  made  it  to  orbit.  The  second,  "C"  sphere,  became  LCS4 
upon  launch.  Both  spheres  were  carefully  measured,  including  the 
mass,  which  was  given  in  the  records  as  38.186  kg  for  the  "G" 
sphere,  and  35.203  kg  for  the  "C"  sphere.  The  two  spheres  were 
unmarked,  and  virtually  indistinguishable.  One  cannot  rule  out 
the  possibility  that  they  were  interchanged  during  ground 
handling,  and  that  the  "G"  sphere  actually  made  it  to  orbit  as 
LCS4 .  In  this  case  all  the  recent  density  models  are  in  better 
agreement  with  the  data.  To  account  for  this,  the  scale  factors 
can  be  multiplied  by  1.085.  Then  assuming  specular  reflection 
the  MSIS83  has  a  mean  scale  factor  of  virtually  unity,  and  the 
J77 1 '  is  0.95.  Both  of  these  values  are  within  the  uncertainty 
of  the  C . .  In  any  event  the  recent  models  are  in  acceptable 
average  agreement  at  780  km,  although  there  are  large  variations 
in  the  density  that  are  not  modeled.  It  is  true  that  the  model 
with  smallest  variability  for  LCS4  is  the  DTM. 

EGS 


EGS  (Experimental  Geodetic  Satellite)  is  a  spherical 
satellite  in  a  nearly  circular  orbit  at  1500  km  altitude,  and 
equipped  with  laser  cube  corner  reflectors.  The  laser  ranging 
data,  provided  by  NASA  since  launch  in  July  1986,  has  allowed  the 
study  of  the  density  model  at  1500  km.  This  is  higher  in 
altitude  than  any  data  used  in  constructing  the  models,  and 
serves  to  measure  how  well  the  predicted  densities  of  the  models 
can  be  extrapolated.  The  EGS  satellite  has  a  known  A/M  that  is 
given  in  Table  8  and  a  surface  made  of  aluminum  with  holes  for 
the  laser  cube  corner  reflectors. 

Preliminary  calculations  on  two  months  of  data  have  been  done 
using  the  MSIS83  and  J77 ' '  models.  The  scale  factors  were 
computed  with  both  a  diffuse  and  specular  scattering  mechanism. 
Since  the  atmosphere  at  1500  km  is  mostly  hydrogen,  one  would 
expect  a  significant  difference  in  the  result.  In  Table  8,  we 
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see  that  for  EGS  the  specular  scattering  assumption  leads  to 
scale  factors  of  1.9  and  greater;  the  MSIS83  mean  value  is  2.39. 
The  conventional  diffuse  scattering  model  gives  an  average  scale 
factor  of  1.39  for  the  J77'*  and  1.67  for  the  MSIS83.  A-priori 
one  would  choose  the  diffuse  scattering  for  aluminum,  and  this  is 
adopted.  It  is  also  concluded  from  this  analysis  that  J771' 
underestimates  the  density  at  1500  km  by  about  30%,  whereas  the 
MSIS83  underestimates  the  density  by  about  60%.  This  is 
consistent  with  the  results  for  LCS4,  where  the  J77''  gives 
larger  densities  than  MSIS83  at  the  780  km  altitude. 

VII.  EVALUATION  OF  ATMOSPHERIC  INDICES 

In  addition  to  evaluating  the  atmospheric  models,  the  derived 
scale  factors  can  be  used  to  determine  how  well  the  effects  of 
the  different  indices  are  being  modeled  and  to  evaluate  whether  a 
new  index  is  of  value  to  future  modelling  efforts.  This  was  done 
by  correlating  the  series  of  derived  scale  factors  with  different 
series  of  atmospheric  indices,  including  the  daily  F10.7  cm  flux, 
the  Kp,  and  the  precipitation  index  (PI) .  The  precipitation 
index  is  the  only  parameter  we  studied  that  is  not  currently 
being  used  as  an  input  to  any  of  the  atmospheric  models.  It  is 
an  index  which  quantifies  the  intensity  and  spatial  extent  of 
high-latitude  particle  precipitation  based  on  observations  made 
along  individual  passes  of  the  NOAA/TIROS  weather  satellites.  It 
is  measured  in  near  real-time  by  the  NOAA/TIROS  weather 
satellites  which  are  in  circular  sun-synchronous  orbits  at  850 
km.  (Foster,  et  al.,  1986).  We  analyzed  the  power  levels  of  the 
precipitation  index. 

The  atmospheric  parameters  were  averaged  over  the  same  time 
interval  as  the  data  used  to  compute  the  orbital  arcs.  For 
COS 1179 ,  3  days  of  data  were  averaged  for  each  parameter,  while 
for  LCS4 ,  4  days  of  data  were  averaged.  Each  series  of 
atmospheric  data  was  then  correlated  against  the  series  of  scale 
factors  determined  for  both  satellites  using  each  of  the  four 
different  atmospheric  models.  In  certain  cases,  two  series  of 
scale  factor  data  existed  for  each  model  corresponding  to 
specular  and  diffuse  scattering.  This  issue  was  discussed 
earlier.  It  was  not  significant  when  correlating  the  series  of 
atmospheric  parameters  against  the  series  of  scale  factors. 

The  results  for  each  atmospheric  parameter  are  presented  in 
Tables  9  a,b,  and  c.  In  all  cases,  the  correlation  coefficients 
were  computed  from  scale  factors  and  geophysical  data  with  the 
average  subtracted,  i.e.  these  data  have  zero  mean.  A 
correlation  coefficient  was  determined  for  each  atmospheric 
model  and  each  satellite.  A  value  of  greater  than  0.2  was 
considered  to  be  significant  correlation.  The  data  discussed  in 
these  next  sections  represent  the  first  six  months  of  data  in 
1985  with  one  exception.  The  bracketed  value  in  the  J77"  columns 
for  11796  Diffuse  represent  the  entire  year  of  data  for  1985. 
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TABLE  9a.  Correlation  Coefficients  for  F10.7  cm  Flux 


11796 


5398 


Specular 

Diffuse 

Specular 

Diffuse 

J71 

0.18 

0.35 

J77" 

-0.14 

-0.15  (-0.02) 

-0.07 

-0.01 

DTM 

1978  0.02 

0.12 

0.12 

MSIS 

1983  0.04 

0.06 

0.11 

0.02 

TABLE  9b.  Correlation  Coefficients  for  Kp 

11796 

5398 

Specular 

Diffuse 

Specular 

Diffuse 

J71 

0.21 

-0.24 

J77 

0.40 

0.39  (0.27) 

0.27 

0.29 

DTM 

0.40 

0.09 

0.14 

MSIS83  0.28 

0.35 

0.28 

TABLE  9c.  Correlation 

Coefficients 

for  the  Precipitation  Index 

11796 

5398 

Specular 

Diffuse 

Specular 

Diffuse 

J71 

0.42 

-0.15 

J77  ' 

'  0.48 

0.46  (0.37) 

0.36 

0.40 

DTM 

0.53 

0.28 

0 . 19 

MSIS83  0.42 

0.42 

0.45 

0.36 

The  results  of  the  correlation  between  the  series  of  daily 
F10.7  cm  values  and  the  series  of  determined  scale  factors  .^re 
given  in  Table  9a.  It  appears  that  for  the  lower  satellite  the 
F10.7  cm  flux  is  modeled  fairly  accurately.  Problems  that 
existed  in  the  J71  model  seem  to  have  been  corrected  in  the  later 
models . 

The  results  of  the  correlation  between  the  Kp  indices  and  the 
scale  factors  are  given  in  Table  9b.  The  Kp  index  is  used  to 
model  the  influence  of  geomagnetic  fluctuations  in  the 
atmosphere.  There  appears  to  be  significant  correlation  between 

and  the  series  of  scale  factors  determined  for  both 


Figure  5  shows  the  correlation  between  the  Kp  data  and  the 
scale  factors  for  the  Jucchia  1977  S  and  the  MSIS  1983  S  models. 
S  indicates  specular  scattering.  The  mean  of  each  data  set  has 
been  determined  and  subtracted  from  the  series.  The  y-  axis 
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satellites . 
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therefore  represents  the  data,  with  the  mean  subtracted,  while 
the  x-axis  represents  time.  This  data  is  presented  here  to 
illustrate  the  features  in  the  data  sets  that  are  correlated. 

Based  on  the  data  for  the  lower  satellite,  the  MSIS  model 
shows  slightly  less  correlation  with  the  Kp  than  does  the 
Jacchia  1977  model,  indicating  that  it  better  models  the 
atmospheric  response  to  changes  in  this  index.  We  should  point 
out  that  the  MSIS  model  uses  the  Ap  index  instead  of  the  Kp 
index.  We  computed  the  correlation  between  the  Ap  index  and 
these  data  sets.  The  results  were  very  similar  to  the  results 
presented  in  Table  9b.  Finally,  the  DTM  model  correlation  with 
Kp  at  275  km  is  similar  to  the  others  (0.40),  but  at  780  km  it 
has  a  significantly  smaller  correlation  (0.09). 

The  results  of  the  correlation  between  the  precipitation  power 
index  (PI)  and  the  scale  factors  are  given  in  Table  9c.  As  is 
evidenced  in  Table  9c,  there  is  significant  correlation  between 
this  index  and  the  scale  factors  of  the  atmospheric  models  and 
both  satellites.  Of  the  parameters  we  investigated,  the  highest 
average  correlation  coefficients  were  seen  with  this  index. 

Figure  6  illustrates  in  greater  detail  the  correlation  between 
the  scale  factor  data  and  the  precipitation  index.  The 
precipitation  index  is  not  used  by  any  of  the  models  to  predict 
the  atmospheric  response  to  geomagnetic  activity.  Based  on  this 
data,  we  think  that  the  precipitation  index  should  be  included  in 
future  atmospheric  models.  One  should  also  note  that  there  is 
significant  correlation  between  the  precipitation  index  and  the 
Kp.  The  correlation  coefficient  between  the  Kp  and  prescription 
data  sets  is  0.47.  Some  of  the  correlation  observed  with  the  Kp 
could  actually  be  leakage  of  the  precipitation  index 
correlation.  Clearly  this  issue  needs  more  data  to  be  resolved. 

VIII.  REVISIONS  TO  THE  JACCHIA  77"  DENSITY  MODEL 

Model  corrections  were  suggested  by  the  fact  that 
significant  correlation  of  our  scale  factors  and  certain 
geophysical  information  were  found .  The  scale  factor  information 
we  determined  can  be  used  to  suggest  revisions  to  the  models.  A 
change  in  the  scale  factor  can  be  shown  to  be  directly  equal  to  a 
change  in  the  log  of  the  atmospheric  density.  What  we  present 
here  is  a  suggested  revision  to  the  J77"  model  based  on  a  least 
squares  fit  of  our  scale  factor  data  to  the  following  terms  : 
annual,  semi-annual,  seasonal-latitudinal,  diurnal,  Kp,  and  PI 
(precipitation  index) .  We  were  able  to  investigate  seasonal 
effects  because  we  had  a  full  year  of  data. 

To  determine  the  corrections  to  the  model  we  chose  the  scale 
factors  associated  with  the  lower  satellite  11796  computed  using 
the  J77"  model.  For  each  scale  factor,  an  averaging  time  and 
values  for  the  latitude,  hour  angle,  solar  declination,  and 
geophysical  index  (PI,  Kp,  F10.7)  were  associated.  The  averaging 
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time  was  chosen  so  that  it  corresponds  to  the  data  span  used  in 
determining  the  scale  factor,  S.  The  geophysical  indices  were 
averaged  over  the  same  time  interval.  The  values  of  latitude, 
hour  angle,  and  solar  declination  were  the  values  for  the 
satellite  at  perigee  during  that  time  period. 

Several  functional  forms  were  postulated  and  a  joint  least 
squares  solution  for  the  parameters  was  made.  The  formal 
statistics  of  each  parameter  were  then  assessed  for  significance. 
In  this  way  the  set  of  statistically  significant  parameters  was 
found.  For  this  series  of  S  the  results  are  given  in  Table  10. 

A  mean  value  and  slope  were  also  recovered  although  not  shown. 


Table  10.  11796  Scale  Factor  Analysis 

d  S  =  d  log  (  p  ) 

Cosine  Sine 

ANNUAL  (-0.0178  +-  0.0063)  (-0.0034  +-0.0107) 

SEMIANNUAL  (-0.0026  +  -  0.0047)  (-0.0629  +-0.0067) 

LATITUDE  (  0.0139  +-  0.0047)  (-0.0089  +-0.0040) 

SEASONAL  LATITUDINAL  (  0.0022  +-  0.0154)  (  0.0399  +-0.0087) 

DIURNAL  (-0.0129  +-  0.0111)  (  0.0035  +-0.0059) 

PRECIPITATION  INDEX 

( (0.0009+-0.0005)+(0.0017+-0, 0009) *COS(0) 2) ) *PI 


To  understand  the  annual  and  semi-annual  variations  in  S, 
the  F10.7  cm  flux,  the  90  day  mean  F10.7  cm  flux,  and  the  scale 
factor  S  data  were  also  analyzed  for  annual  and  semiannual 
variations.  These  coefficients  are  given  in  Table  11.  The 
annual  variations  in  the  F10.7,  the  mean  F10.7,  and  S  are  more 
or  less  in  phase.  The  annual  corrections  in  S  might  be  expected, 
as  there  are  known  to  be  long  period  variations  for  a  number  of 
geophysical  observations.  For  example,  there  are  decade  scale 
variations  in  weather,  earth  rotation,  the  ocean  circulation, 
upper  atmospheric  circulation  and  climate.  This  variation 
probably  cannot  be  predicted,  and  must  be  measured  for  each  year 
and  during  each  solar  cycle. 
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TABLE  11.  Annual  and  Semiannual  Variation  (Epoch  =  MJD  46325) 

ANNUAL  SEMIANNUAL 


AMP 

PHASE 

AMP 

PHASE 

SCALE  FACTOR 

0.0181 

-169.19 

0.063 

-92.37 

PI 

2.5820 

162.50 

1.519 

-30.55 

F10.7 

1.5970 

-120.40 

3.690 

147.60 

[F10.7] 

2.5200 

-143.40 

2.440 

158.10 

Kp 

0.1807 

51.05 

0.094 

-12.98 

It  is  possible  that  the  latitude,  seasonal-latitudinal,  and 
diurnal  variation  are  all  related  to  the  shape  of  the  diurnal 
variation  (the  diurnal  bulge) . 

The  most  interesting  result,  is  the  dependence  of  S  on  the 
precipitation  index,  PI.  The  PI  is  not  an  input  to  any  model  at 
present.  Our  analysis  indicates  that  it  shows  promise.  In  the 
modej.  form,  experiments  were  done  with  an  hour  angle  dependence 
of  the  PI,  but  none  was  found  significant.  Recall  that  there  is 
significant  correlation  between  the  PI  and  the  Kp  index.  The 
correlation  coefficient  between  these  two  data  sets  is  0.47. 

Both  the  Kp  and  PI  are  measures  of  the  solar  particle  input 
energy  to  the  atmosphere.  What  is  unclear  is  how  much  of  the 
information  is  common  to  both.  Note,  that  we  did  not  find  a 
significant  parametric  dependence  of  S  on  Kp.  However,  when  the 
dependence  on  PI  was  determined,  and  subtracted  from  S,  the 
residual  correlation  of  S'=S-dS  with  Kp  was  reduced  to  nearly 
zero.  Recall  that  the  correlation  of  S  with  Kp  had  been  0.27 
(see  Table  9b.)  From  this  preliminary  analysis  it  would  seem 
that  Kp  and  PI  have  somewhat  different  information  to  contribute. 
Furthermore  the  functional  dependence  of  PI  grows  toward  the 
equator  (0. 0017*cos (Q) 2)  which  is  counter  intuitive. 

Conventional  wisdom  holds  that  the  influx  of  energy  is  in  the 
auroral  zone,  and  the  effect  of  the  thermosphere  propagates 
equatorward  decreasing  in  size.  In  fact  the  Jacchia  77"  model 
has  a  sin(0)3  dependence  on  Kp,  which  is  zero  at  the  equator.  A 
recent  revision  (Slowey  1984)  gives  a  more  detailed  analysis  of 
the  geomagnetic  effect  and  this  is  no  longer  true.  Our  model 
dependence  may  only  indicate  a  flawed  model  assumption  in  the 
J77"  model.  The  fact  remains  that  there  are  still  fundamental 
questions  about  the  particle  flux  that  should  be  answered. 

With  the  seasonal  terms  removed,  the  correlation  of  scale 
factors  with  the  precipitation  index  is  reduced  from  0.37  to 
0.33.  The  correlation  of  the  scale  factors  with  the  other 
geophysical  parameters  does  not  change.  The  Kp  correlation,  for 
instance,  remained  at  0.27.  Apparently,  the  annual  variation  in 
-he  precipitation  index  is  in  phase  with  the  error  in  the  J77" 
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density  model  This  leads  us  to  conclude  that  long  term 
variation  (seasonal)  in  the  atmosphere  may  be  controlled  by  the 
same  phenomenon  that  influences  the  precipitation  index. 

The  overall  result  of  the  scale  factor  analysis  is  summarized 
in  Table  12.  There  are  significant  (10%)  changes  in  the  Jacchia 
77  model  suggested  by  this  preliminary  analysis.  Most  of  the 
changes  would  be  considered  revisions,  and  they  must  be  tested 
against  other  satellite  data  before  taking  them  too  seriously. 

The  precipitation  index  model  is  new,  and  could  be  quite 
important.  Here  too,  one  must  clearly  understand  the  relation 
between  Kp  and  PI  before  making  any  final  model.  We  believe  that 
we  have  only  scratched  the  surface  and  that  there  is  significant 
insight  to  be  gained  from  this  type  of  analysis. 


TABLE  12.  11796  SCALE  FACTOR  ANALYSIS 

1.  Correlation  of  S  with  Kp  is  through  cross  correlation  of 

Kp  with  PI.  There  is  no  direct  correlation  of  Kp  with 
scale  factor.  1 

2.  Significant  Model  Corrections  for  Semiannual  (6.3%), 
Latitude  (1.6%),  Seasonal  Latitudinal  (4.0%),  and  Diurnal 
(1.3  %)  are  appropriate  for  this  satellite.  The  latitude 
correction  coud  be  the  same  as  the  diurnal  correction. 

3.  Precipitation  index  predicts  as  much  as  10  %  change  in 
density. 

4.  Precipitation  index  effect: 

(O.OO17*cos(0) 2  +  0.0009)  *  PI  *  F ( z ) 

is  concentrated  at  equator  and  minimum  at  the  pole. 

This  is  opposite  to  the  geomagnetic  Kp  which  depends  on 
sin(0)3  in  the  J77"  model. 

5.  Precipitation  index  effect  seems  planetary,  not  depending 
on  hour  angle. 

IX.  SUMMARY 

We  have  tested  atmospheric  density  models  at  275  km,  780 
km,  and  1500  km.  At  the  two  lower  altitudes,  all  models  are  in 
good  agreement  in  the  sense  that  they  give  the  same  average 
performance.  They  all  exhibit  significant  departures  from  the 
actual  density,  and  it  is  necessary  to  include  "solve  for" 
parameters  in  orbit  determination  to  match  the  tracking  data.  As 
our  analysis  continues,  the  differences  between  models  becomes 
more  subtle,  and  a  selection  becomes  more  difficult.  For 
example,  the  finding  of  Gaposchkin  and  Coster  (1986)  of  a  large 


difference  between  the  J77 1 '  and  the  MSIS83  at  780  km  is 
significantly  reduced  here.  From  the  standpoint  of  computational 
speed  the  J71  is  clearly  superior,  and  the  MSIS83  is  deficient. 

We  are  continuing  work  in  five  areas: 

1.  To  extend  the  data  for  all  three  satellites,  to 
avoid  any  biasing  due  to  seasonal  effects,  and  to 
improve  the  sampling. 

2.  To  determine  the  mass  of  C0S1179  and  LCS4  using 
satellite  orbit  perturbations.  It  is  at  present 
routine  to  determine  the  mass  of  high  altitude 
satellites  from  perturbations  due  to  solar  radiation 
pressure.  This  has  not  been  possible  until  now  for 
lower  satellites  because  the  effects  of  geopotential 
model  errors  dominated  the  solution.  Recent 
geodetic  solutions  show  promise  of  changing  this 
situation.  In  this  case  we  can  hope  to  make  some 
statement  about  the  absolute  densities  at  275  km 
altitude,  and  clarify  the  ambiguity  about  the  mass  of 
LCS4 . 

3.  Explore  alternate  indices  or  variables  as  input  to 
the  density  models,  such  as  the  precipitation  index. 

4.  Expand  the  source  of  accurate  tracking  data  on  these 
satellites,  in  order  to  increase  the  space  and  time 
resolution  of  the  drag  determinations. 

5.  Examine  the  accuracy  of  f c"~  ;ting  drag,  for 
prediction  of  satellite  orbits. 

6.  Obtain  data  on  other  low  altitude  satellites. 

X.  CONCLUSIONS 

1.  No  model  does  an  adequate  job  of  modelling  the 
atmospheric  density. 

2.  There  is  no  agreement  on  what  is  the  best  model.  We 
find  that  the  differences  between  models,  though 
measurable,  are  less  than  the  agreement  between 

the  models. 

3.  There  are  real  physical  variations  in  the  atmosphere 
that  are  not  modeled  by  any  of  the  current  suite  of 
atmospheric  models.  New  model  parameters  are  needed 
(e.g.,  winds,  gravity  waves,  (Gross,  1985). 

4.  The  inclusion  of  the  precipitation  index  in  future 
atmospheric  models  should  be  investigated. 
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Significant  correlation  was  observed  between  the 
precipitation  index  and  the  scale  factors  at  275  and 
780  km. 

5.  Overall  these  models  have  at  most  18%  difference 
about  the  mean  when  averaged  of  all  latitudes  and 
hour  angles  below  800  km.  However,  the  variance  of 
the  model  differences  exceed  16%. 

6.  The  J71  continues  to  perform  exceptionally  well,  and 
is  the  fastest  overall.  The  MS1S83  is  by  a  large 
factor,  the  slowest. 

7.  For  our  overall  use,  balancing  accuracy,  computer 
speed,  and  range  of  height,  we  use  the  J77'1 

8.  The  final  conclusion  is  that  satellite  drag  data  can 
still  play  a  role  in  understanding  the  thermosphere, 
and  in  contributing  unique  data  to  monitor  the  next 
solar  cycle. 

XI .  RECOMMENDATIONS 

OBTAIN  NEW  DRAG  DATA  WITH  EXISTING  FACILITIES 

PRECISION  RADAR  DATA 

CALIBRATION  SPHERES  (14  SATELLITES  NOW  AVAILABLE) 

RECOVER  MEAN  DENSITIES 

MODEL  IMPROVEMENTS 

MONITOR  IN  NEAR  REAL  TIME  ATMOSPHERIC/GEOPHYSICAL  INDICES 

EXISTING  DATA  IS  AVAILBALE  (E.G.,  TRANSIT  SATELLITES) 

MAKE  NEW  MEASUREMENTS  OF  ENERGY  INPUTS 

CONTINUE  PRECIPITATION  INDEX  MEASUREMENTS 

OTHER  INDICES  NOW  MAINTAINED 

SATELLITE  MEASUREMENTS  OF  EUV  AND  CORRLE LAT I ON  WITH  F10.7 

CM  FLUX, ETC. 

WIND  AND  GRAVITY  WAVES  SHOULD  BE  INCORPORATED  IN  MODELS 

LAUNCH  DRAG  SATELLITES  (BALLOONS)  SUPPORTED  WITH  HIGH  ACCURACY 

TRACKING 

NEW  SATELLITE  EXPERIMENTS  ARE  NECESSARY,  SHOULD  BE  SIMULTANEOUSLY 

CALIBRATED  WITH  DRAG  MEASUREMENTS. 


Ill 
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DRAG  EFFECT  FOR  (A/M)  =  0.1  cm2/gram 
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Accelerometers  flown  on  low-altitude  satellites  have 
provided  our  most  extensive  density  data  base  in  the 
1 bO- 250  km  region.  These  data  have  been  obtained 
over  a  wide  region  of  solar  and  geomagnetic  condi¬ 
tions.  A  comprehensive  comparison  of  our  measure¬ 
ments,  obtained  1974-1982,  with  several  empirical 
models,  has  resulted  in  assessment  of  the  current 
status  and  improvement  prospects  of  our  capability 
to  specify  and  forecast  density  (dray).  The  model 
accuracies  are  compared  by  their  mean  values  and 
standard  deviations  relative  to  the  accelerometer 
data.  Results  indicate  a  lack  of  significant  model 
accuracy  improvements  during  the  past  two  decades. 

Mean  values  are  estimated  to  be  given  typically 
within  +_  10%  and  standard  deviations  are  approxi¬ 
mately  "+  15%.  Improved  representat i on  of  atmos¬ 
pheric  Variability  requires  a  program  including 
coordinated  measurements  of  atmospheric  structure 
and  dynamics,  development  of  more  accurate  indica¬ 
tors  of  solar  and  geomagnetic  activity  and  theore¬ 
tical  investigations  using  self-consistent 
numerical  models. 

INTRODUCTION 

Empirical  atmospheric  density  models  based  on  large  data  sets  have 
been  developed  to  describe  variations  of  the  upper  atmosphere.  In 
atmospheric  research  these  models  serve  as  a  reference  and  permit 
analysis  of  di screpanc i es  with  respect  to  new  measurements.  For 
practical  applications  models  are  also  required  for  problems  concerned 
with  the  effects  of  aerodynamic  drag  on  satellites.  Satellite  lifetime, 
desiyn,  control,  tracking,  on-board  fuel  requirements  and  reentry  are 
all  affected  by  atmospheric  dray  and  its  variability.  Two  versions 
of  models  are  currently  used.  Those  of  Jacchia  [1-4J  are  based 
mainly  on  total  density  data  derived  from  satellite  orbital  decay 
observations.  Models  based  on  in  situ  composition  data  and  temperatures 
inferred  from  ground-based  incoherent  scatter  radar  have  been  formulated 
by  Hedin  et  al  [5-6J  and  Hedin  |7-HJ.  A  common  feature  of  all  the 
models  is  the  use  of  relatively  simplified  physical  concepts.  Varia¬ 
tions  are  described  as  a  function  of  altitude,  solar  and  geomagnetic 
activity,  latitude,  longitude  [4,  b-8|,  local  time  and  day  of 
year . 


tThi  s  paper  also  presented  as  AAS  paper  No.  87-552,  AAS/AIAA  Astro- 
dynamics  Specialist  Conference,  Kalispell,  Montana,  Aug  10-13,  1987 
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Several  limitations  of  present  models  with  regard  to  solar  and  geo¬ 
magnetic  variations  are  recognized.  Heating  due  to  solar  EUV  radia¬ 
tion  activity  is  taken  into  account  by  relations  involving  ground- 
based^  measurements  of  the  10.7  cm  solar  flux,  F^o.7»  and  some  mean 
value,  F,  measured  over  several  solar  rotations,  this  index  does  not 
necessarily  represent  the  complex  mechanisms  of  interactions  between 
the  solar  extreme  ultraviolet  (EUV)  flux  and  the  thermosphere.  Since 
it  has  been  found  to  generally  relfect  variations  in  the  thermospheric 
energy  input,  it  is  routinely  used  as  a  readily  available,  but 
imperfect,  indicator  of  solar  activity.  Studies  to  determine  the 
relationship  between  F^q#7  and  satellite  measurements  of  EUV  [9]  reveal 
correlations  that  are  non-linear  and  both  wavelength  and  solar 
cycle  dependent.  Similarly,  the  planetary  3-hr  Kp  or  daily  Ap  index 
used  as  a  geomagnetic  activity  indicator  does  not  necessarily 
represent  the  physical  mechanisms  responsible  for  density  variations. 
The  indices  are  derived  from  ground-based  magnetic  field  fluctuations. 
They  are  used  as  indicators  of  heating  at  high  latitudes  caused  by 
the  interaction  of  solar  plasma  with  the  upper  atmosphere.  However, 
it  is  the  only  routinely  available  index  for  geomagnetic  activity. 

Based  on  a  network  of  midlatitude  stations,  it  is  particularly 
limited  when  describing  data  at  high  latitudes.  The  need  for 
improved  geomagnetic  indicators  is  particularly  borne  out  by  the 
data  presented  in  this  study. 

Lack  of  comprehensive  experimental  data  over  the  range  of  variables 
described  in  models  also  contributes  to  their  limitations:  (a)  The 
dominance  of  the  semidiurnal  tide  below  200  km  has  recently  been 
revealed  by  low  altitude  satellite  measurements.  [10-11 J  This  feature 
is  incorporated  in  the  MS  I S  models  but  not  in  those  of  Jacchia. 

However,  the  considerable  phase  and  amplitude  structure  of  the 
semidiurnal  tide  as  a  function  of  altitude  and  latitude  is  beyond 
the  scope  of  these  models,  (b)  Density  variations  with  annual  and 
sub-annual  periodicities  are  character i zed  as  a  "semiannual 
variation"  with  maxima  near  equinox  and  minima  near  solstice.  The 
cause  of  this  variation  has  not  been  unambiguously  resolved. 

Measured  data  show  irregu’ar  variations  from  year  to  year  in  maxima 
and  minima,  (c)  Satellite  composition  measurements  have  provided 
evidence  of  longitudinal  variations.  Density  maxima  for  heavy 
constituents  occur  at  the  longitude  of  the  north  and  south  magnetic 
poles  while  the  behavior  of  the  lighter  constituents  is  negatively 
correlated.  This  behavior  is  related  to  ionospheric  plasma  density 
and  motion  influenced  by  earth's  magnetic  field.  The  MSIS  models 
rely  on  geographic  coordinates  and  relatively  few  terms  to  describe 
the  broad  features  of  this  variation,  (d)  Large  scale  gravity  wave 
structures  also  occur.  These  have  periods  different  from  those  of 
tidal  waves  and  are  not  directly  related  to  the  earth-moon-sun 
geometry.  They  may  be  generated  either  locally  in  the  thermosphere 
or  at  lower  altitudes.  Gravity  waves  with  the  laryest  amplitudes 
(~  70%  peak-to-peak )  observed  at  high  altitudes  detected  by  accel¬ 
erometers  on  the  0V1-15  satellite,  [12]  were  related  to  high-latitude 
heating  as  indicated  by  the  auroral  electrojet  index.  [13]  Wave  struc- 
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tures  are  not  incorporated  into  present  models,  (e)  Latitude 
dependences,  associated  with  all  of  the  atmospheric  variations,  are 
not  well  understood.  A  particular  problem  for  the  yeomaynet i c 
activity  effect  is  that  the  time  delay  between  storm  onset  and  den¬ 
sity  variation  is  latitude-dependent. 

In  the  lower  thermosphere ,  below  about  250  km,  model  deficiencies 
are  particularly  acute  since  they  are  based  primarily  on  data  at 
higher  altitudes.  Physical  and  mathematical  difficulties  limit 
extrapolation  of  these  data  to  lower  altitudes:  the  assumption  of 
diffusive  equilibrium  is  not  always  valid,  and  the  lower  thermos¬ 
phere  boundary  conditions  are  poorly  known  and  variable.  AFGL  is 
accumulating  an  extensive  atmospheric  density  data  base,  from  satel¬ 
lite  accel erometer  measurements,  for  the  150-240  km  region  [14J# 

The  present  study  uses  these  data  to  determine  the  validity  of 
empirical  models  and  establishes  phenomena  not  contained  in  them. 
Definitive  low-altitude  satellite  measurement  programs,  to  guide 
development  of  computationally  efficient  dynamic  models  utilizing 
realistic  indicators  of  atmospheric  energy  sources,  are  recommended. 
The  effects  of  high  latitude  processes  on  lower  thermosphere 
dynamics  are  emphasized. 

DATA  DESCRIPTION 

Extensive  measurements  of  the  lower  thermosphere  neutral  density 
have  been  derived  using  the  satellite  accel erometer  experiment. 
Instrument  operating  principles  have  been  described  by  Marcos  and 
Swift  [15J. Density  is  derived  from  direct  measurements  of  aero¬ 
dynamic  dray  toyether  with  knowledge  of  the  satellite's  mass,  area, 
velocity  and  dray  coefficient.  Figure  1  shows  the  flight  history 
of  the  accelerometer  from  1974  to  the  present.  All  spacecraft  had 
near-polar  orbits  except  for  AE-C  (68°)  and  AE-E  (19°).  This  exten¬ 
sive  lower  thermosphere  data  set  has  been  obtained  over  a  wide  ranye 
of  latitudes  and  solar  and  yeomaynetic  conditions.  Data  from  the 
first  seven  flights  have  been  used  in  the  present  analysis  of 
empirical  models.  The  satellites  used  and  the  dates  of  data  acquisi¬ 
tion  are  yiven  in  Table  1.  The  altitude  ranye  of  the  data  is 
yenerally  150-240  km  for  AE-C,  -D,  and  -E  and  S3-1  and  170-240  km 
for  S3-4  and  SETA-1  and  -2 

Data  reduction  techniques  implemented  for  the  SETA  flight 
data  are  described  in  Reference  15.  Because  of  the  high  lenyth-to- 
diameter  ratio  of  the  host  vehicle,  the  equations  used  for  determin¬ 
ation  of  atmospheric  density  (  o  )  and  cross-track  winds  Lib]  are: 

A  ref  ) 

aD  = - Ci  p\T  (2) 

2M 

where  i  =  x,  y,  z 

M  =  satellite  mass 
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Aref  =  satellite  frontal  area 
V  =  -  vG  +  vA  +  Vw 

with  Vq,  V/\  and  Vw  representing  respectively  inertial  satel¬ 
lite  velocity,  the  atmospheric  rotation  velocity  (assumed 
equal  to  earth's  rotation  velocity)  and  the  neutral  wind 
velocity. 

Ci  =  drag  coefficients;  (Vx,  V y,  Vz). 

The  equation  for  the  component,  in  a  particular  direction,  of  the 
total  force  on  an  element  of  area  is  given  by: 

1  1  2C2 

dC  = - •  {(rk  +  yi  +  nt)  Ly  ■  (1  +  erf  yS)  +•  —  e-Y  s  J 

Aref  S2  ft 


l  ^r  y/n  1  2C2 

+  —  (1+erf  yS)  +  -  L —  (1+erf  yS)  +  --e“Y  s  J}  dA  , 
2S^  2  Ti  S  S2 

where 

2  dF 


ArefP ^ 

F  =  aerodynamic  force 

f,  ,y  ,n  =  direction  cosines  between  mass  velocity  vector  and  axes  of 

element  of  area,  respectively 

k,y,t  =  direction  cosines  between  the  axes  of  element  of  area,  and  the 

direction  in  which  force  is  desired  (accelerometer  axis) 

erf  YS  =  error  function  of  yS 

S  =  molecular  speed  ratio  =  (mass  velocity  of  gas)/(most  probable 

random  speed  of  molecules) 

=  temperature  of  incident  molecules 

Tr  =  temperature  of  reflected  molecules  (assumed  =  300°K) 

This  equation  is  integrated  over  the  entire  satellite  surface  area 
to  obtain  the  drag  coefficient  for  a  particular  axis.  Values 
appropriate  to  tnis  satellite  were  used  with  the  MS  I S  model  and  the 
above  equation  to  derive  drag  coefficients.  For  AE-C,  -D  and  -E 
and  S3-1,  which  had  1 ength-to-di ameter  ratios  near  unity,  a  constant 
value  of  Cq=2.2  was  used. 


123 


Fig.  1.  Satellite  accelerometer  flight  history  and  solar  activity 
v  9 •  time. 


MODEL  EVALUATIONS  USING  ACCELEROMETER  DATA 


The  models  considered  were  those  of  Jacchia,  designated  J65  L 1 J » 

J70,  [2]  J71,  [31  J73  (unpublished,  see  Ref  17)  and  J77,  [4J 
Jacchi a-Wal ker-Bruce ,  L17J  Densel  ,  [17J  Lockheed-NASA,  L 1 H  |  U.S. 
Standard  Atmosphere,  1952,  L19J  U.S.  standard  Atmosphere  Supplements, 
19b6  L 20 J  and  the  Mass  Spectrometer  Inco  herent  Scatter  models 
designated  MSIS/7,  L  b J  MSIS79  L6J,  MSIS83  L  7 J  and  MSIS86  L8J.  Two 
versions  of  the  1983  MSIS  model  were  tested:  MS1S83A  used  a  daily 
averaged  Ap  index  for  the  geomagnetic  activity  effect  and  MSIS83B 
used  a  weighted  average  of  the  3-hr  index  over  a  59  hour  interval. 

The  latter  version  of  MSIS86  was  used  in  this  study. 

Comparisons  of  model  predictions  with  measurements  not  used  in  the 
model  construction  are  valuable  for  testing  the  model  reliability. 
While  some  AE  mass  spectrometer  data  were  used  in  developing  the 
MSIS  models  the  density  data  were  not.  Also  the  S3-1  and  -4  and 
the  SETA-1  and  -2  data  were  not  used  in  construction  of  any  of  the 
models.  The  quantity  statistical ly  analyzed  is  the  ratio  (R) 
between  measured  density  and  model  density.  Results  are  given  for 
mean  ratios  (R)  and  standard  deviations  (  a  p)  in  Table  2.  Models 
are  listed  in  chronological  order.  This  Table  indicates  the  lack 
of  significant  improvement  in  model  accuracies  during  the  past  two 
decades . 

The  results  of  Table  2  indicate  some  of  the  difficulties  in  improv¬ 
ing  the  accuracy  of  empirical  models.  The  MSIS79  model  added 
longitude  variations,  based  on  direct  satellite  composition  measure¬ 
ments,  to  the  MSIS77  model.  This  was  done  without  changing  the 
zonal ly  averaged  values.  Incorporation  of  an  observed  phenomenon 
resulted  in  a  less  accurate  model  (by  0.1  to  0.5%  in  r  )  for  every 
data  set.  The  MSIS83  models  are  based  on  a  data  set  much  more 
extensive  than  that  used  for  MSIS77.  While  a  small  accuracy  improve¬ 
ment  results  in  MSIS83B  from  a  weighted  average  vs  a  daily  average 
for  MSIS83A,  neither  model  is  generally  more  accurate  than  MSIS77. 

As  part  of  a  revised  and  updated  J71  model,  a  more  complex  geomag¬ 
netic  variation  was  incorporated  into  the  J7/  model.  For  all  data 
sets,  the  standard  deviations  of  J77  are  0.4  to  2.8%  higher  than 
those  of  J71.  Revisions  made  to  the  J71  coefficients  and  formulas 
to  produce  the  J73  model  resulted  in  the  most  accurate  of  the 
Jacchia  models;  yet,  this  model  was  never  published.  MSIS86  is 
very  similar  to  MSI583  but  with  additional  terms  to  represent 
seasonal  differences  in  polar  composition  variations.  The 
modifications  are  based  mainly  on  Dynamics  Explorer  satellite 
composition  data  obtained  above  300  km.  The  model  provides  typically 
a  1%  reduction  in  standard  deviations  over  MSIS83  in  the  150-240  km 
region.  However,  the  MSIS86  errors  remain  comparable  to  those  in 
J 7 1  and  J73. 
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TABLE 


2.  ACCELEROMETER,  TOTAL  MASS  DENSITY  RATIOS  TO  MODELS 
(ALTITUDE  150-240  KM) 


AE 

-C 

AE 

-D 

AE 

-E 

S3 

-1 

R 

*R 

R 

*R 

R 

<*r 

R 

Or 

MSIS86B 

1.10 

14.5 

0.95 

15.1 

1.02 

13.0 

1.03 

14.4 

MSIS83B 

1.11 

15.0 

0.98 

15.6 

1.02 

13.0 

1.07 

14.6 

MSIS83A 

1.12 

15.2 

0.99 

16.2 

1.03 

13.6 

1 .08 

14.8 

MS1S79 

1.05 

14.2 

0.98 

16.5 

1.01 

13.6 

1.00 

14.7 

MSI S77 

1.05 

14.0 

0.98 

16.3 

1.01 

13.5 

1.00 

14.2 

377 

1.05 

15.9 

1.01 

15.3 

1.02 

15.3 

1.04 

14.3 

373 

1.07 

14.0 

1.02 

14.7 

1.07 

13.7 

1.06 

13.5 

371 

1.10 

14.7 

1.05 

14. B 

1.08 

15.0 

1.08 

13.7 

370 

1.05 

17.1 

0.98 

15.4 

1.00 

15.8 

1.04 

14.6 

J64 

0.94 

17.0 

0.89 

17.4 

0.92 

15.5 

0.92 

17.8 

L-N 

0.94 

17.7 

0.86 

16.5 

0.86 

15.8 

0.94 

14.9 

JWB 

0.99 

19.5 

0.99 

18.8 

1.01 

19.4 

0.99 

19.6 

US66 

0.95 

16.5 

0.89 

15.6 

0.92 

15.5 

0.95 

14.4 

US62 

0.84 

30.0 

0.70 

32.4 

0.72 

30.3 

0.73 

32.7 

DENS 

1.45 

21.3 

1  .00 

21.0 

0.94 

23.0 

1.31 

19.5 

Number 

Points 

of 

56908 

28273 

33835 

25825 

S3- 

4 

S  STA¬ 

■1 

SETA- 2 

R 

** 

FF  c 

R 

R 

Or 

MSIS86B 

1.04 

11.1 

MSI S83B 

0.98 

11.8 

MSIS83A 

0.99 

12.1 

MSIS79 

0.98 

11.5 

MSIS77 

0.98 

11.2 

J77 

0.94 

13.7 

J73 

0.96 

11.6 

J71 

0.99 

12.1 

J70 

0.97 

12.0 

J64 

0.90 

11.6 

L-N 

0.93 

14.6 

JWB 

0.86 

11.0 

US66 

0.90 

11.5 

US62 

0.93 

17.6 

DENS 

0.79 

19.9 

1.01 

9.8 

0.94 

11.0 

0.92 

9.7 

0.87 

11.6 

0.93 

10.1 

0.88 

11.8 

0.96 

11.7 

0.92 

11.7 

0.96 

11.5 

0.92 

11.2 

0.88 

12.6 

0.89 

13.9 

0.92 

9.8 

0.92 

10.2 

0.94 

9.9 

0.95 

10.1 

0.99 

9.3 

0.93 

10.4 

0.99 

11.1 

0.88 

11.3 

0.99 

9.9 

0.91 

11.4 

0.90 

10.4 

0.82 

10.9 

0.99 

11.0 

0.88 

11.2 

1 . 13 

12.0 

0.96 

15.0 

0.87 

14.9 

0.99 

15.2 

Number  of 
Points  38215 


56530  277442 
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More  detailed  analyses  of  the  statistical  properties  of  the  SETA 
accelerometer  data  have  been  made  by  partitioning  data  (to  240  km) 
into  b°  geographic  latitude  bins  and  day/night  local  time  bins  as 
well  as  four  kp  four  bins.  Local  time  bins  of  08-16  hr  and  16-24 
hr  were  used  to  separate  SETA  dayside  and  nightside  northern  hemis¬ 
phere  data  occurring  at  the  same  latitude. 

Mean  values  relative  to  the  J71  MSIS83B  and  MSIS86  models  for 
SETA-1  data  (Mar-Apr  79)  are  shown  in  Figs.  2-4  respectively  as  a 
function  of  geographic  latitude.  Part  a  of  each  Figure  consists  of 
daytime  measurements  (  ~  1000  hours  local  time)  and  part  b  consists 
of  nighttime  measurements  (  ~  2200  hours  local  time).  The  horizontal 
axis  from  right  to  left,  follows  the  satellite  trajectory.  The 
satellite  altitude  as  a  function  of  latitude  typically  decreases 
from  240  km  at  low  northern  hemisphere  latitudes  on  the  nightside 
to  200  km  near  the  pole  down  to  about  170  km  near  30°N  on  the 
dayside.  It  then  increases  to  about  180  km  at  the  equator.  Since 
the  orbit  is  not  quite  polar,  no  data  are  available  for  the  85-90° 
bins.  Model  ratios  indicate  dayside  relative  density  enhancements 
near  20°N  and  65°N.  Nighttime  densities  are  underestimated,  partic¬ 
ularly  near  30°N,  oy  10-15%.  J71  tends  to  overestimate  the  geo¬ 

magnetic  activity  response  at  all  latitudes  and  particularly  at 
night.  The  high  Kp  bin  data  are  about  8%  (dayside)  and  15%  (night¬ 
side)  lower  than  those  of  the  low  Kp  bin.  The  MSIS  models  show 
good  agreement  of  the  low  and  high  Kp  bins  on  the  dayside  and  about 
a  7%  overestimation  on  the  nightside. 

Mean  values  of  SETA-2  data  (Jul-Nov  82)  are  compared  to  J71,  MS83B 
and  MSS86B  in  Figs.  5-7  respectively.  Peaks  in  ratios  near  65°N  on 
the  dayside  and  increasing  ratios  with  decreasing  latitude  on  the 
nightside  are  again  prominent  features  of  all  data  sets.  The 
geomagnetic  activity  representation  of  MSIS  is  less  accurate  than 
for  SETA-1  data.  Maximum  density  response  overestimations  occur 
for  the  highest  Kp  bin  at  high  latitudes  on  the  nightside.  J71 
provides  a  better  estimation  of  the  nightside  density  response  than 
it  did  for  SETA-1  data  and  than  MSIS83B  for  the  SETA-2  data. 

Figs.  8  and  9  show  the  standard  deviations  for  J71  and  MSIS83 
respectively.  These  Figures  show  a  striking  latitude  dependence. 

The  major  errors  occur  at  auroral  and  polar  latitudes,  even  during 
geomagnetical ly  quiet  conditions.  The  maximum  error  of  about  15% 
occurs  for  the  highest  Kp  bin  in  both  cases.  Similar  results 
were  shown  by  Marcos  [13J  for  the  SETA-2  and  S3-4  data  relative  to  the 
J71  model.  The  MSIS86  data  (not  shown)  give  a  similar  result. 

The  previous  analyses  show  that  lower  mean  values  and  standard 
deviations  are  obtained  with  the  SETA-1,  -2  and  -3  and  S3-4  satel¬ 
lite  data  set.  However,  systematic  differences  between  the  data 
from  these  satellites  and  the  data  from  the  AE/S3-1  satellites  are 
to  be  expected  114]. 
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Mean  values  of  SETA-1  data  to  J71  model  plotted  as  a  function 
of  geographic  latitude  (four  Kp  bin6). 
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(b) 


Fig.  3 


Mean  valus  of 
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SETA  1  data  to  MS1S  83  model  plotted  as  a  function 
latitude  (four  Kp  bins). 
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Fig.  U  Mean  values  of  SETA  1  data  to  MS1S  86  model  plotted  as  a  function 
of  geographic  latitude  (four  Kp  bins). 
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Fig.  5  Mean  values  of  SETA-2  data  to  J71  model  plotted  as  a  function 
of  geographic  latitude  (four  Kp  bins). 
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Fig.  6  Mean  values  of  SETA-2  data  to  MSIS  83  model  plotted  as  a  function 
of  geographic  latitude  (four  Kp  bins). 
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Fig.  7  Mean  values  of  SETA-2  data  to  MSTS  86  model  plotted  as  a  function 
of  geographic  latitude  (four  Kp  bins). 
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Fig.  8 


Standard  deviations  of  ratios  of  SETA- 1  data  to  J 7 1  model 
plotted  as  a  function  of  geographic  latitude  (four  Kp  bins). 
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Fig- 


9  Standard  deviations  of  ratios  of  SETA-1  dati  to  MSIS  83  model 
plotted  as  a  function  of  geographic  latitude  (four  Kp  bins). 
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The  Table  2  AE/S3-1  ratios  imply  densities  generally  above  model 
values  while  the  SETA  ratios  imply  densities  below  model  values. 

This  difference  between  the  two  satellite  data  sets  is  believed  due 
mainly  to  uncertainties  in  the  ballistic  coefficients  used  to 
convert  orbital  drag  measurements  to  atmospheric  density  values. 

Drag  coefficient  (Cq)  estimates  contribute  the  major  uncertainty  to 
ballistic  coefficient  errors.  Our  estimated  errors  in  Cp  are  _+ 

10%.  It  is  possible  that  the  SETA  C(j  errors  exceed  +_  15%  due  to  the 
spacecraft  shape  121 J, Therefore,  systematic  errors  in  drag  coefficient 
may  completely  account  for  the  absolute  value  differences  between 
AE/S3-1  and  SETA  total  mass  densities.  We  conclude  that,  on  average, 
recent  models  reflect  density  to  within  the  jh  10%  accuracy  of  the 
Atmosphere  Explorer  accelerometer  measurements. 

The  standard  deviation  difference  has  been  attributed  to  be  the 
manner  in  which  altitude,  latitude  and  local  time  variations  were 
sampled  for  each  data  set.  The  AE/S3-1  data  were  sampled  over  a  mix 
of  altitudes,  latitudes  and  local  times.  For  SETA  data,  these 
parameters  were  highly  correlated;  e.g.,  mid-latitudes  were  sampled 
near  170  km  and  1030  hrs  LT  on  the  dayside.  AE-E  and  S3-1  data 
have  been  further  analyzed  at  low  latitudes  to  show  quant i tat i vely 
that  reducing  the  range  of  variables  over  which  a  model  must  predict 
density  leads  to  a  reduced  standard  deviation  in  agreement  with  the 
SETA  values. 

The  AE-E  satellite  acquired  data  between  _+  20°  latitude  over  a 
period  of  about  one  year.  During  this  time  data  were  sampled  over 
four  local  time  cycles.  Data  obtained  in  the  spinning  mode  were 
analyzed  in  altitude,  latitude,  Kp,  local  time  and  seasonal  bins. 
Results  presented  in  Table  3,  show  reductions  in  o  associated  with 
these  bin  limitations.  The  AE-E  170  -  175  km  data,  grouped  by 

season  show  values  of  between  5%  and  7%  in  good  agreement  with 

SETA  low  latitude  values.  A  similar  analysis  was  performed  with 
the  S3-1  data.  Results  of  binning  are  given  in  Table  4  for  170  - 

175  km  0-20°N  and  for  170  -  200  km  for  0-80°N.  Local  time  was  near 

0730  hours  for  this  data  set.  Note  that  for  0-10°N,  170  -  200  km, 
o  is  6.5%  in  excellent  agreement  with  the  Table  4  AE-E  data. 

Therefore  the  typically  15%  standard  deviations  obtained  from  the 
AE/S3-1  data,  not  the  lower  values  obtained  with  the  S3-4/SETA 
flights,  are  representati ve  of  model  errors  over  "all"  conditions. 

DISCUSSION 

Accelerometer  data  have  revealed  unmodeled  features  of  lower  thermo¬ 
spheric  structure.  Particularly  significant  are  observed  variations 
at  high  latitudes  as  a  function  of  latitude,  local  time  and  level 
of  geomagnetic  activity.  The  importance  of  new  knowledge  at  auroral 
and  polar  latitudes  is  borne  out  by  the  maximum  standard  deviaations 
of  models  in  those  regions.  However,  additional  density  data  alone 
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Table  3.  Atmosphere  Explorer  -  E  Statlatlcal  Result*  with  Local  Tine,  Latitude,  Season, 
Altitude  and  Kp  Blna  Spinning  Orbit  Data 


J71  «IS  83 


Altitude 

Lat 

Kp 

LT 

Us 

N. 

R 

0 

R 

0 

170-175 

120 

All 

All 

All 

1929 

0.088 

0.136 

0.028 

0.108 

L 

0+10 

0-3+ 

All 

All 

280 

0.104 

0. 142 

0.025 

0.100 

▼ 

120 

0-3+ 

10-11 

DJF 

47 

0.058 

0.078 

-0.041 

0.071 

Table  4.  S3-1  Statistical  Result*  va.  Latitude;  K.p  0  to  3+ 


J71 

MS  IS 

83 

ttean 

Altitude 

Lat 

LT 

N. 

R 

a 

R 

a 

170- 

175 

0-10 

0736 

51 

0.214 

0.066 

0.058 

0.068 

170- 

175 

10-20 

0737 

66 

0.206 

0.059 

0.076 

0.066 

170- 

200 

0-10 

0736 

136 

0. 194 

0.065 

0.016 

0.073 

10-20 

737 

244 

0.180 

0.061 

0.025 

0.073 

20-30 

0736 

188 

0.138 

0.075 

0.026 

0.068 

30-40 

0733 

134 

0.132 

0.084 

0.068 

0.086 

40-50 

0734 

101 

0. 153 

0.104 

0.122 

0.088 

50-60 

0748 

226 

0.095 

0.117 

0.106 

0.098 

60-70 

0805 

306 

0.038 

0.081 

0.090 

0.076 

1 

r 

70-80 

0846 

232 

0.131 

0.124 

0.172 

0. 109 
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will  not  be  sufficient  to  improve  models.  More  recent  models, 
utilizing  extensive  satellite  measurements  of  atmospheric  composi¬ 
tion  have  been  shown  to  specify  total  density  with  about  the  same 
accuracy  as  older  models  based  on  satellite  orbital  decay  measure¬ 
ments  . 

Problems  in  accurately  describing  the  spat i a  1 /tempora 1  character¬ 
istics  of  density  variability  are  indicated  by  the  diagram  of  Fig. 

10.  The  complex  phenomena  contributing  to  atmospheric  heating  and 
dynamics  are  dominated  by  two  sources  of  solar  energy:  electromag¬ 
netic  radiation  and  particles.  Solar  EUV,  the  most  important 
thermospheric  heat  source,  is  deposited  mainly  at  low  latitudes  and 
in  the  summer  hemisphere.  The  circulation  and  structure  of  the 
thermosphere  are  controlled  primarily  by  this  heating.  The  processes 
involved  are  reasonably  well  understood.  At  auroral  latitudes,  the 
situation  is  different.  The  thermosphere  is  driven  by  el ectrodynaini c 
coupling  between  the  solar  wind,  magnetosphere,  ionosphere  and 
neutral  atmosphere.  Energetic  charged-particle  precipitation  from 
the  magnetosphere  and  Joule  (resistance)  heating  due  to  electric 
current  Systems  within  che  ionosphere  are  the  major  energy  sources. 
Also,  ionospheric  plasma,  driven  by  the  magnetospheri c  convection 
electric  field,  transfers  momentum  to  the  neutral  particles  and  can 
produce  winds  of  about  one  km  sec-^  .  Roth  the  'Magnitude  and  spa¬ 
tial  extent  of  the  auroral  processes  are  extremely  variable  and 
related  to  the  level  of  geomagnetic  activity.  During  large  storms, 
the  energy  dumped  into  the  relatively  narrow  auroral  region  can 
exceed  the  global  EUV  input.  These  loc  lized  heat  sources  can 
cause  large  and  sudden  density  variations,  perturb  circulation 
patterns  and  launch  gravity  waves.  High  latitude  processes  can 
therefore  significantly  alter  thermospheri c  density,  composition 
and  temperature  on  a  global  scale  as  a  result  of  large  geomagnetic 
disturbances.  Hence,  density  observed  at  a  given  location  can 
depend  on  both  the  past  history  of  density  at  that  location  and  on 
events  that  occurred  at  varying  distances  and  times. 

Dynamic  models  are  now  being  developed  [22  and  zJJ.  These  are 
theoretical  three-dimensional  and  t ime- dependent  models  that  numeri¬ 
cally  solve  the  primitive  energy  «nd  momentum  equations.  They 
attempt  to  understand  the  physics,  chemistry  and  dynamics  respon¬ 
sible  for  thermospheric  structure.  Parameter! zat ions  of  various 
input  fields  are  necessary  for  the  specification  of  toe  physical 
processes  described  in  these  models.  These  inputs  include  solar  UV 
and  EUV  energy  and  their  heatin'-  efficiencies  and  energy  and 
momentum  sources  resulting  from  the  solar  wind-magnetosphere  inter¬ 
action.  The  output  is  a  giooal  simulation  of  the  structure  and 
dynamics  of  the  thermosphere.  !n  their  present,  form,  the  models 
would  not  be  feasible  for  operational  use  since  they  require  a  CRAY 
computer.  However,  these  humeri.,  a !  '  1  s  provide  the  means  for 

guiding  development  o'  new  generation  of  operational  models 
utilizing  realistic  Physical  concepts. 


SETA-1  density  and  winds  obtained  during  the  March  22,  1979 
geomagnetic  storm  have  provided  a  unique  opportunity  to  test  and 
evaluate  our  knowledge  and  modeling  capability  of  the 
magnetosphere/ionosphere/thermosphere  system.  This  period  was  one 
of  intense  solar-terrestrial  observations  and  analysis.  As  a 
result,  the  time-dependent  variations  of  auroral  and  magnetospheric 
convection  model  parameterizations  incorporated  within  the  National 
Center  for  Atmospheric  Research  thermospheric  general  cirri' *tlc 
model  (TGCM)  were  prescribed  in  detail  as  input  to  the  storm  simula¬ 
tion  L 24 J .  The  TGCM  predictions  of  density  variability  for  this 
storm  period  are  compared  L 25 J  with  the  satellite  data  and  the 
MSIS83  model  in  Fig  11.  Three-dimensional  depictions  are  gSen  for 
latitudes  above  45°  in  the  northern  hemisphere  for  loc;l  times 
~  2200  hrs  (left  side)  and  ~  1100  hrs  (right  side).  In  the  day 
sector,  the  SETA  measurements  clearly  confirm  the  TGCM  predictions 
of  a  high-latitude  (60°-85°)  localized  region  of  intense  heating 
and  density  enhancement  (  ~  60%)  between  1400  and  1500  UT  (Universal 
Time),  the  period  of  maximum  magnetic  disturbance.  The  night  sector 
is  characteri zed  by  more  wavelike  structures  and  density  perturbations 
of  20-40%  in  both  the  SETA  and  TGCM  data.  The  extension  of  the 
TGCM  "hot  spot"  farther  into  the  nighttime  sector  than  is  consistent 
with  SETA  data  is  believed  due  to  inaccuracies  in  the  TGCM  convection 
model  position,  orientation  or  latitudinal  extent.  Over  the  complete 
range  of  latitudes,  Ref.  25  shows  that  the  TGCM  and  SETA  data  are 
characteri zed  by  generally  elevated  densities  (10-20%)  after  the 
storm  (2200  UT)  as  compared  with  pre-storm  conditions  (0800  hrs 
GMT).  The  MSIS  model  predicts  density  enhancements  of  only  10-20% 
between  1100  UT  and  1800  UT,  without  the  pronounced  latitudinal  and 
temporal  structure  indicated  by  the  TGCM  and  the  SETA  data.  It  is 
evident  that  the  TGCM  accounts  for  time-dependent  localized  processes 
with  a  much  higher  time  resolution  than  the  MSIS83  model. 

CONCLUSIONS 

Satellite  accelerometer  measurements  have  provided  the  first  exten¬ 
sive  description  of  density  variability  at  low  satellite  altitudes 
and  high  latitudes.  The  data  have  been  obtained  over  a  wide  range 
of  solar  and  geomagnetic  conditions.  Statistical  evaluations  of 
density  results  from  seven  satellites  demonstrate  that  model  accura¬ 
cies  have  not  progressed  significantly  in  the  past  20  yeans.  Mean 
values  are  given  within  *  10%  with  standard  deviations  of  *  15%. 

Errors  encountered  for  specific  applications  have  been  shown  to  be 
sensitive  to  geographic  location,  local  time  and  level  of  geomag¬ 
netic  activity  as  well  as  to  the  model  used.  Development  of  better 
models  will  result  from  improvements  in  our  knowledge  of  the  magni¬ 
tude,  distribution  and  transport  of  thermospheric  heating.  This 
requires  coordinated  theor^t  ical  ,/observat  i  jnal  programs.  Specific 
recommendat i ons  for  the  lower  thermosphere  include: 


Extensive  interpretation  and  analyses  of  the  accel erometer  data 
base.  Emphasis  should  be  placed  on  case  studies  of  spati al /temporal 
variations  related  to  yeomaynetic  activity  as  a  function  of  solar 
flux,  season,  latitude,  lonyitude  and  local  time. 

Continued  development  of  dynamic  ylobal  circulation  models. 

Accel erometer  density  and  winds  data  as  well  as  Dynamics  Explorer 
satellite  data  are  needed  to  validate  and  improve  input  assumptions. 

More  realistic  indices  to  describe  energy  inputs  related  to 
auroral  processes  and  their  interaction  with  the  maynetosphere/ 
ionosphere  system.  Auroral  electrojet  parameters  and  new  hiyh 
latitude  indices  derived  from  yround-based  maynetometers  can  be 
evaluated  with  accel erometer  density  data  for  the  same  periods. 

An  extensive  iower  thermosphere  measurements  proyram  that 
provides  a  basis  for  developiny  accurate  relationships  between 
energy  inputs  and  ylobal  density  variability.  An  initial  payload 
measuring  density,  composition,  winds,  temperature  and  ionospheric 
properties  is  required.  A  low  eccentricity  sun-synchronous  orbit 
(1400/0200  LT)  would  provide  data  in  local  time-latitude  regions 
not  covered  by  the  present  data  set.  Coordination  with  ground- 
based  thermosphere  dynamics  measurements  programs  (such  as  the 
Coupled  Energetics  and  Dynamics  of  Atmospheric  Regions  program) 
will  greatly  enhance  data  i nterpretat 1  on  . 

The  above  recommendations  have  focussed  on  applications  of  accelero¬ 
meter  data  to  improve  density  models.  Since  these  data  measure  the 
parameter  to  be  modeled,  atmospheric  density,  they  must  be  an 
integral  part  of  future  programs  that  attempt  to  improve  density 
models.  Continued  progress  in  numerical  modeling  and  in  geomagnetic 
index  development  together  with  1 nterpretat i ons  of  present  data  and 
a  critical  new  measurements  program  can  provide  important  advances  in 
thermospheric  models. 
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MSIS  Empirical  Model  Status  and  Directions  for  Improvement 
by  A.  E.  Hedin  (NASA/GSFC) 

for  AFSPACECOM  Atmospheric  Neutral  Density  Specialist  Conference 
Colorado  Springs,  CO  March  22-23,  1988 

I.  Introduction 

With  the  detection  of  atmospheric  drag  effects  in  early 
artificial  satellites,  data  became  available  which  led  to  the  well  known 
series  of  Jacchia  models.  The  J65  (also  called  J64)  model  was  the 
earliest  comprehensive  global  model  based  on  satellite  drag  (orbital 
decay).  This  model  was  the  first  to  include  the  principal  types  of 
thermospheric  density  variations  and,  unlike  the  standard  atmospheres, 
provides  an  approximation  to  the  true  spatial  and  temporal  variations  of 
the  thermosphere.  While  drag  based  models  use  temperature  and 
composition  as  intermediate  parameters  for  the  calculation  of  total 
density,  these  intermediate  parameters  may  be  in  error  since  drag  data 
provide  no  direct  information  on  composition  and  temperature.  The  J77 
model  Introduced  separate  pseudo-temperatures  for  each  constituent ,  but 
this  significant  artificial  complication  was  only  partially  successful 
in  reproducing  composition  variations. 

The  OGO-6  satellite  mass  spectrometer  launched  in  1969  provided 
the  first  extensive  measurements  of  the  densities  of  molecular  nitrogen, 
atomic  oxygen,  and  helium  in  the  thermosphere.  Their  sum  provides  an 
independent  determination  of  total  density.  The  observed  variations  in 
composition  could  be  quite  different  from  the  Jacchia  model  predictions 
and  led  to  a  new  approach  to  represent  the  observed  variability.  The 
ground  based  Incoherent  scatter  measurements  of  temperature  were 
subsequently  combined  with  the  in  situ  composition  measurements  in  the 
MSIS  (Mass  Spectrometer  and  Incoherent  Scatter)  models  to  provide 
composition  and  temperature  as  well  as  total  density  predictions  of 
equivalent  accuracy  for  various  geographical,  temporal,  and  solar 
conditions.  The  latest  model,  MSIS-86,  was  chosen  for  CIRA  1986,  but 
publication  of  CIRA  1986  has  been  delayed. 

II .  Present  accuracy 

Comparisons  of  measured  total  densities  with  various  models 
(Tables  1  H  2  show  mean  absolute  differences  and  standard  deviations 
for  magnetically  quiet  conditions,  and  Table  3  shows  standard  deviations 
for  magnetically  active  conditions,  all  in  terms  of  log  to  base  e) 
indicate  a  general  accuracy  of  15%  to  20%  in  the  200  to  400  km  range  and 
25  to  30%  at  higher  altitudes.  Data  is  scarce  at  lower  altitudes,  but 
what  is  available  suggests  model  accuracies  in  the  15  to  20%  range. 
Ranking  the  models  from  1  to  5  (1  best)  for  each  data  set  in  selected 
altitude  ranges  and  plotting  (Figs.  1  to  3)  for  each  model  how  many  data 
sets  (counting  each  altitude  range  separately)  are  best  (or  second  best, 
etc.)  represented  by  that  model  shows  that  there  is  little  difference 
between  models  with  respect  to  absolute  differences  between  data  and 
model.  With  respect  to  standard  deviations  between  data  and  model,  the 
later  (MSIS)  models  have  a  better  ranking  with  more  data  sets  than  the 
earlier  Jacchia  models.  However,  the  difference  between  the  worst  and 


best  model  is  often  only  a  few  percentage  points. 

The  data  set  labeled  Jacchia  drag,  for  instance,  consists  of 
densities  determined  from  orbital  decay  by  the  Smithsonian  group  and 
were  used  in  generating  the  Jacchia  models  but  not  the  MSIS  models.  The 
comparison  results  point  to  both  consistency  (within  15%)  between  two 
quite  different  techniques  (mass  spectrometer  and  drag)  and  between  two 
different  time  periods,  since  the  drag  data  generally  come  from  an  older 
time  period  than  the  mass  spectrometer  measurements. 

While  total  density  predictions  have  improved  relatively  little, 
the  predictions  of  composition  and  temperature  have  improved 
considerably  (Table  4  and  Fig.  4).  The  most  dramatic  improvements  are 
for  the  minor  constituents  in  each  altitude  range,  but  there  are 
improvements  at  altitudes  where  a  species  is  not  minor  also.  For 
example,  in  the  200  to  400  km  range  predictions  by  J70  for  both  N2  and  0 
density  are  worse  than  predictions  by  MSIS-86,  yet  the  differences  in 
total  density  are  small  because  the  errors  in  J70  model  predictions  for 
these  two  major  gases  approximately  compensate.  The  J70  model  helium 
densities  are  not  as  good  as  MSIS-66  above  800  km  (where  He  is  a  major 
constituent),  while  J70  compares  better  to  drag  densities  than  MSIS. 
Solving  this  puzzle,  which  may  be  due  to  errors  in  the  drag  coefficient 
which  would  appear  as  errors  in  the  Jacchia  densities  or  the  possible 
existence  of  a  hot  oxygen  population  which  has  not  been  taken  into 
account,  would  certainly  improve  calculation  of  drag  force  at  the  higher 
altitudes.  Since  the  drag  coefficient  depends  on  atmospheric  composition 
and  temperature  it  is  important  to  predict  these  quantities  correctly  in 
order  to  calculate  drag  force  or  infer  density.  The  accuracy  of  the 
composition  at  one  altitude  also  determines  the  accuracy  of  extrapolated 
densities  using  diffusive  equilibrium. 

Winds  at  high  latitudes  can  change  drag  fo^ce  by  10  to  20%  but 
this  effect  has  not  been  included  in  deductions  of  density  nor  in 
tracking  errors  because  no  suitable  wind  model  existed.  The  Dynamics 
Explorer  satellite  measured  winds  over  the  whole  globe  and  these 
together  with  Atmosphere  Explorer  wind  measurements  have  been  r  «d  to 
create  the  first  empirical  model  of  global  winds.  An  example  for 
solstice  conditions  is  shown  in  Fig.  5  where  the  solid  contours  show 
temperature  and  the  dashed  contours  pressure.  The  winds  blow 
approximately  from  high  to  low  pressure  areas  but  are  higher  in  the 
polar  region  due  to  magnetospheric  forcing.  Note  that  temperature  and 
pressure  (as  well  as  density)  do  not  maximize  at  the  same  location  and 
this  is  contrary  to  a  basic  assumption  of  the  early  Jacchia  models.  The 
winds  shown  are  in  addition  to  the  rigid  body  rotation  of  the ' atmosphere 
with  the  earth.  Super-rotation  '■ffects  (a  net  or  zonally  averaged  flow 
beyond  the  rigid  body  rotation  i  3ed)  are  not  of  great  practical 
importance  for  drag.  Because  the  effect  on  drag  force  depends  on  the 
direction  the  satellite  is  moving  with  respect  to  the  wind,  the  wind 
effects  may  not  have  seriously  corrupted  a  density  model  if  there  were  a 
sufficient  variety  of  orbits  involved  in  the  model  database. 

Waves  provide  a  basic  floor  to  achievable  prediction  accuracy 
which  available  measurements  indicate  is  in  the  2  to  10%  range  depending 
on  latitude  (larger  at  high  latitudes).  A  model  predicting  the 
likelihood  of  waves  of  a  certain  amplitude  in  the  upper  thermosphere  is 


currently  feasible.  Unless  the  lower  thermosphere  boundary  can  be 
described  in  some  detail  with  observations,  the  lower  atmosphere  weather 
will  also  provide  some  noise  level  in  the  thermosphere  (particularly  the 
lower  thermosphere)  which  is  not  well  established. 

III.  Calculation  efficiency 

The  following  times  were  found  for  total  density  calculation  on 
a  MicroVax  computer: 

MSIS-86  47  millisec 
J70  40  millisec 

J77  123  millisec 

These  calculations  did  not  use  any  optimizing  algorithms.  Order  of 
magnitude  timing  improvements  are  possible  for  the  Jacchia  models 
because  of  the  unique  density-temperature  relationship  which  is  assumed 
in  their  formulation.  However,  it  is  this  very  assumption  which  is  a 
basic  cause  of  errors  in  composition  and  temperature  in  these  models. 

The  MSIS  timing  can  be  Improved  by  about  a  factor  of  2  for  total  density 
by  skipping  the  calculation  of  constituents  which  are  of  minor 
Importance  for  the  particular  altitude. 

It  is  likely  that  the  timing  of  any  of  the  models  can  be 
improved  for  particular  applications  where  certain  types  of  variations 
can  be  disregarded. 

IV.  Directions  for  accuracy  Improvement 
IVa.  General 

Although  engineering  and  operational  interests  frequently  focus  on 
total  density,  an  effort  to  improve  predictions  of  total  density  without 
regard  to  composition,  temperature,  and  wind  is  not  particularly 
interesting  to  the  wider  scientific  community  because  it  Is  not  likely 
to  improve  our  understanding  of  the  atmosphere  and  its  processes.  In  any 
case,  more  accurate  density  by  itself  may  not  lead  to  higher  accuracy 
drag  force  determinations  since  the  drag  coefficient  depends  on 
composition  of  the  atmosphere  and  drag  depends  on  wind  as  well. 
Spacecraft  design  also  depends  more  and  more  on  accurate  composition  for 
problems  like  material  degradation  and  surface  glow. 

Improvements  are  likely  to  come  in  small  steps  rather  than  dramatic 
breakthroughs  since  the  current  accuracy  limit  seems  to  have  a  number  of 
diverse  causes  such  as  measurement  accuracy  of  historical  data, 
differences  between  techniques,  undiscovered  systematic  variations, 
unpredictable  natural  variability,  and  difficult  to  predict  (with 
limited  ancillary  data)  effects  of  magnetic  storms.  It  should  also  not 
be  forgotten  that  current  models  already  account  for  a  large  proportion 
of  the  observed  density  variance.  Improved  models  of  total  density  and 
other  atmospheric  parameters  will  require  a  sustained  program  involving 
a  wide  range  of  experimental  and  theoretical  approaches  and  must  include 
new  measurements. 
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IVb.  Existing  Data 

Existing  data  (from  previous  scientific  satellites,  ground  based 
techniques,  and  possibly  unexploited  satellite  tracking  data)  should  be 
analyzed  for  a  better  understanding  of: 

magnetic  storm  effects; 
semiannual  variations; 

EUV  (F10.7)  variations; 
interactions  between  known  effects; 
differences  between  solar  cycles; 
differences  between  techniques. 

This  is  the  easiest  approach  and  shows  promise  of  making  some 
improvement  in  the  near  term.  These  studies  will  also  benefit  from 
ongoing  scientific  programs  involving  ground  based  measurements  of 
parameters  other  than  density. 

IVc.  New  Measurements 

A  model  can  only  be  as  good  as  the  available  data.  No  theory  is 
good  enough  to  calculate  atmospheric  parameters  without  reference  tc 
measurements.  Furthermore,  we  have  no  certainty  that  the  relations  we 
derive  (such  as  between  density  and  the  10.7  cm  solar  flux)  hold  for  the 
indefinite  future. 

New  density  and  other  measurements  must  be  designed  to  surpass 
the  required  accuracy.  Densities  must  be  corrected  for  wind  and  be 
cross-checked  by  different  techniques  to  check  basic  assumptions.  If 
global  predictions  are  required  then  measurements  must  be  global  in 
scope.  Since  we  know  the  semiannual  variation  is  different  from  year  to 
year,  we  assume  the  lower  boundary  will  be  different  from  year  to  year 
(if  not  day  to  day),  and  we  are  not  sure  what  happens  in  each  new  solar 
cycle,  measurements  must  be  taken  reasonably  near  the  required  time. 


IVd.  Theoretical  Support 

Theoretical  analysis  and  Global  Circulation  Models  (GCM) : 
provide  guidance  for  organizing  and  interpreting  measurements  and  fill 
in  where  measurements  are  missing  (such  as  the  lower  thermosphere). 
Current  models  should  be  evaluated  for  their  current  ability  to  predict 
data.  Just  as  empirical  models.  Theoretical  models  may  be  particularly 
valuable  for  highly  dynamic  effects  and  may  possibly  be  combined  with 
empirical  models.  An  important  question  is  how  accurate  and  detailed 
our  knowledge  must  be  of  the  auroral  energy  inputs  in  order  that  the 
additional  detail  predicted  by  theoretical  models  is  sufficiently 
correct  in  phase  and  amplitude  as  to  be  more  useful  than  a  smoothed 
representation. 
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IVe.  Ancillary  Measurements 

To  provide  energy  input  for  theoretical  models  and  correlation 
parameters  for  empirical  models  we  must  continue  or  improve  measurements 
such  as : 

solar  radio  flux; 
solar  EUV  flux; 

IMF; 

magnetic  indices; 
auroral  morphology; 
particle  precipitation. 

The  present  plans  to  measure  an  EUV  band  on  a  long  term  basis 
from  the  GOES  satellite  is  an  encouraging  step.  While  there  may 
be  some  prediction  improvement  possible  by  correlating  this  new 
measurement  with  the  10.7  cm  flux  and  thus  introducing  it  into  the 
existing  models,  it  will  only  be  through  direct  correlation  with  new 
measurements  of  density  and  temperature  that  the  full  usefulness 
of  this  measurement  can  be  realized  and  validated. 

To  both  interpret  drag  data  in  terms  of  density  and  make  better 
predictions  of  drag  from  density  models  there  is  a  need  for  improved 
knowledge  of  the  drag  coefficient  under  a  wide  variety  of  conditions 
through  laboratory  and  space  experiments. 
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Overall  mean  (magnetically  quiet)  (preliminary) 


Alt 

Pts  MSIS 

-86 

MSIS 

-83 

MSIS 

-77 

J77 

J70 

Jacchia 

200-400 

3197 

.05 

(3) 

-.05 

(3) 

-.21 

(4) 

-.03 

(1) 

-.04 

(2) 

drag 

400-800 

6516 

.05 

(1) 

-.06 

(2) 

-.23 

(3) 

-.05 

(1) 

-.06 

(2) 

800-1200 

3386 

.04 

(2) 

.04 

(2) 

-.04 

(2) 

-.03 

(1) 

-.05 

(3) 

Barlier 

120-000 

1050 

.  02 

(1) 

-.03 

(2) 

-.02 

(1) 

-.04 

(3) 

-.07 

(4) 

drag 

200-400 

4761 

.04 

(3) 

.  02 

(1) 

-.06 

(4) 

-.04 

(3) 

.03 

(2) 

400-800 

1447 

.  10 

(2) 

.  10 

(2) 

-.04 

(1) 

.  10 

(2) 

.  15 

(3) 

AE-C  MESA 

120-200 

5746( .25) 

.  11 

(3) 

.  11 

(3) 

.05 

(2) 

.05 

(2) 

.00 

(1) 

accel 

200-250 

6101( . 55) 

.02 

(1) 

.04 

(3) 

-.08 

(4) 

.03 

(2) 

-.03 

(2) 

AE-C  OSS 

120-200 

6447( .65)- 

.01 

(1) 

-.01 

(1) 

-.05 

(3) 

-.02 

(2) 

-.07 

(4) 

ms 

200-890 

6279( .2) 

.02 

(2) 

.06 

(3) 

-.08 

(4) 

.09 

(5) 

.01 

(1) 

AE-D  MESA 

120-200 

1 1024( . 7 )  - 

.03 

(3) 

.00 

(1) 

.00 

(1) 

.02 

(2) 

-.02 

(2) 

accel 

200-250 

4399 

.08 

(3) 

-.04 

(2) 

-.  11 

(4) 

.01 

(1) 

-.04 

(2) 

AE-D  OSS 

120-200 

11787 

.01 

(1) 

.01 

(1) 

.  01 

(1) 

.04 

(2) 

.01 

(1) 

ms 

200-390 

10923( .75)- 

.06 

(3) 

-.02 

(1) 

-.20 

(4) 

.03 

(2) 

-.03 

(2) 

AE-E  MESA 

120-200 

11455( .45) 

.02 

(2) 

.02 

(2) 

-.01 

(1) 

-.01 

(1) 

-.03 

(3) 

accel 

200-250 

8363 

.01 

(2) 

.00 

(1) 

-.06 

(4) 

.05 

(3) 

-.01 

(2) 

AE-E  OSS 

120-200 

1 1812( .62)- 

.01 

(2) 

.00 

(1) 

-.07 

(4) 

-.05 

(3) 

-.07 

(4) 

ms 

200-400 

9741 ( . 12)- 

.03 

(2) 

-.03 

(2) 

-.  15 

(3) 

-.01 

(1) 

-.03 

(2) 

AE-E  NACE 

120-200 

7533 

.05 

(1) 

-.05 

(1) 

-.05 

(1) 

-.05 

(1) 

-.07 

(2) 

ms 

200-400 

10158( . 14)- 

.09 

(3) 

-.09 

(3) 

-.22 

(4) 

-.02 

(2) 

-.07 

(1) 

400-600 

10243( .20)- 

.  17 

(4) 

-.14 

(3) 

-  .32 

(5) 

-.12 

(2) 

-.08 

(1) 

l-< 

5 

l-l 

B 

l-i 

6 

1-7 

1-! 

5 

2-' 

7 

2-' 

7 

2 

2 

2-9 

l 

2- 

10 

3-' 

7 

3-1 

S 

3- 

3 

3-4 

3-! 

3 

4- 

1 

4-1 

D 

4- 

9 

4-0 

1 

4-; 

3 

5- 

1  , 

5-1 
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Table  1 


Overall  Standard  Deviation  (magnetically  quiet)  (preliminary) 


Alt  Pts  MSIS-86  MSIS-83  MSIS-77  J77  J70 


Jacchia 

200-400 

3197 

.  15 

(1) 

.  16 

(2) 

.  18 

(3) 

.  15 

(1) 

.16 

(2) 

drag 

400-800 

6516 

.26 

(2) 

.26 

(2) 

.29 

(3) 

.25 

(1) 

.26 

(2) 

800-120G 

1  3386 

.27 

(3) 

.26 

(2) 

.30 

(4) 

.29 

(3) 

.22 

(1) 

Barlier 

120-200 

1050 

.22 

(1) 

.23 

(2) 

.22 

(1) 

.22 

(1) 

.23 

(2) 

drag 

200-400 

4761 

.20 

(1) 

.20 

(1) 

.21 

(2) 

.21 

(2) 

.21 

(2) 

400-800 

1447 

.31 

(1) 

.31 

(1) 

.33 

(3) 

.32 

(2) 

.31 

(1) 

AE-C  MESA 

120-200 

5746( .25) 

.  16 

(2) 

.  16 

(2) 

.  15 

(1) 

.  18 

(3) 

.  18 

(3) 

accel 

200-250 

6101( .55) 

.21 

(1) 

.21 

(1) 

.21 

(1) 

.21 

(1) 

.21 

(1) 

AE-C  OSS 

120-200 

6435( .65) 

.  12 

(2) 

.  12 

(27 

.  10 

(1) 

.  13 

(3) 

.  13 

(3) 

ms 

200-390 

6279( .2) 

.  14 

(1) 

.  14 

(1) 

.  14 

(1) 

.  15 

(2) 

.  16 

(3) 

AE-D  MESA 

120-200 

11024( .7) 

.  15 

(2) 

.  15 

(2) 

.  14 

(1) 

.  15 

(2) 

.  14 

(1) 

accel 

200-250 

4399 

.  18 

(1) 

.  19 

(2) 

.  18 

(1) 

.  18 

(1) 

.  18 

(1) 

AE-D  OSS 

120-200 

11787 

.11 

CD 

.  11 

(1) 

.  11 

(1) 

.  13 

(2) 

.  11 

(1) 

ms 

200-390 

10923( .75) 

.  17 

(1) 

.18 

(2) 

.18 

(2) 

.  17 

(1) 

.  17 

(1) 

AE-E  MESA 

120-200 

11455( .45) 

.  12 

(1) 

.  12 

(1) 

.  12 

(1) 

.  14 

(2) 

.  14 

(2) 

accel 

200-250 

8363 

.  18 

(1) 

.  18 

(1) 

.  19 

(2) 

.20 

(3) 

.20 

(3) 

AE-E  OSS 

120-200 

11812( .62) 

.15 

(3) 

.  14 

(2) 

.  13 

(1) 

.  16 

(4) 

.  17 

(5) 

ms 

200-400 

9741 ( . 12) 

.19 

(2) 

.  18 

(1) 

.21 

(4) 

.20 

(3) 

.21 

(4) 

AE-E  NACE 

120-200 

7533 

.  11 

(1) 

.11 

(1) 

.  15 

(2) 

.  17 

(4) 

.  16 

(3) 

ms 

200-400 

10158( . 14) 

.20 

(2) 

.  19 

(1) 

.20 

(3) 

.20 

(2) 

.20 

(2) 

400-600 

10243( .20) 

.  19 

(2) 

.  19 

(2) 

.22 

(4) 

.  18 

(1) 

.20 

(3) 

1-12 

1-10 

1-10 

1-7 

1-7 

2-7 

2-11 

2-4 

2-7 

2-6 

3-2 

3-0 

3-4 

3-5 

3-6 

4-0 

4-0 

4-3 

4-2 

4-1 

5-1 
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Table  2 


Overall  Standard  Deviation  (magnetically  active)  (preliminary) 


Alt  Pts  MSIS-86  MSIS-83  MSIS-77  J77  J70 


Jacchia 

200-400 

2464 

.  18 

(1) 

drag 

400-800 

4007 

.29 

(2) 

800-1200 

2549 

.29 

(2) 

Barlier 

120-200 

1050 

.22 

(2) 

drag 

200-400 

2767 

.  20 

(1' 

400-800 

824 

.  30 

(2) 

AE-C  MESA 

120-200 

12359( .25) 

.  16 

(2) 

accel 

200-250 

1 1585 ( . 55) 

.21 

(2) 

AE-C  OSS 

120-200 

12370( .65) 

.  13 

(2) 

ms 

200-390 

6279( .2) 

.  15 

(1) 

AE-D  MESA 

120-200 

8907 ( .7) 

.  15 

(1) 

accel 

200-250 

3723 

.  19 

(1) 

AE-D  OSS 

120-200 

7608 

.  13 

(1) 

ms 

200-390 

7476 ( . 75) 

.22 

(2) 

AE-E  MESA 

120-200 

8289( .45) 

.  13 

(1) 

accel 

200-250 

5582 

.  17 

(1) 

AE-E  OSS 

120-200 

7868( .62) 

.  16 

(2) 

ms 

200-400 

7S54( . 12) 

.  i.9 

(1) 

AE-E  NACE 

120-200 

6766 

.  13 

(1) 

ms 

200-400 

9079( . 14) 

.  23 

(2) 

400-600 

7987 ( .20) 

.  23 

(1) 

1-11 

2-10 

3- 0 

4- 0 


.19  (2) 

.20  (3) 

.21  (4) 

.  18 

(1) 

.30  (3) 

.31  (4) 

.29  (2) 

.28 

(1) 

.30  (3) 

.33  (5) 

.32  (4) 

.23 

(1) 

.22  (2) 

.22  (2) 

.23  (3) 

.21 

(1) 

.21  (2) 

.22  (3) 

.23  (4) 

.21 

(2) 

.29  (1) 

.32  (4) 

.31  (3) 

.30 

(2) 

.16  (2) 

.15  (1) 

.21  (4) 

.  18 

(3) 

.21  (2) 

.20  (1) 

.23  (3) 

.21 

(2) 

.13  (2) 

.11  (1) 

.16  (4) 

.  14 

(3) 

.16  (2) 

.17  (3) 

.19  (4) 

.  16 

(2) 

.15  (1) 

.15  (1) 

.17  (2) 

.  15 

(1) 

.19  (1) 

.19  (1) 

.20  (2) 

.  19 

(1) 

.13  (1) 

.13  (1) 

.15  (3) 

.  14 

(2) 

.21  (1) 

.21  (1) 

.23  (3) 

.21 

(1) 

.13  (1) 

.13  (1) 

.16  (3) 

.  15 

(2) 

.18  (2) 

.18  (2) 

.20  (4) 

.  19 

(3) 

.15  ( 1  ) 

.15  (1) 

.19  (4) 

.  18 

(3) 

.20  (2) 

.22  (4) 

.23  (5) 

.21 

(3) 

.13  ( 1 N 

.17  (3) 

.19  (4) 

.  16 

(2) 

.22  (1) 

.24  (3) 

.23  (2) 

.  23 

(2) 

.23  (1) 

.27  (2) 

.23  (1) 

.23 

(1) 

1-10 

1-9 

1-1 

1-1 

8 

2-9 

2-3 

2-4 

2-1 

8 

3-2 

3-5 

3-6 

3- 

5 

4-0 

4-3 

4-9 

4-1 

0 

5-1 

5-1 
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Table  3 


Ovt,iall  Standard  Deviation  (magnetically  quiet) 
MSIS  selected  temperature  and  composition  data 


Alt 

Pts 

MSIS- 

-86 

MSIS- 

-83 

MSIS- 

-77 

J77 

J70 

Tempera- 

200-400 

7608 

92 

(1) 

99 

(3) 

97 

(2) 

102 

(4) 

110 

(6) 

ture 

400-800 

1363 

96 

(2) 

106 

(4) 

95 

CD 

102 

(3) 

122 

(5) 

He 

120-200 

804 

.28 

(1) 

.  29 

(2) 

.32 

(3) 

.47 

(4) 

.85 

(5) 

Density 

200-400 

7906 

.27 

(1) 

.28 

(2) 

.36 

(3) 

.42 

(4) 

.87 

(5) 

400-800 

5344 

.26 

(1) 

.27 

(2) 

.  34 

(3) 

.40 

(4) 

.62 

(6) 

800-999 

313 

.  17 

(1) 

.  19 

(2) 

.  39 

(4) 

.44 

(5) 

.37 

(3) 

0  +  02 

120-200 

854 

.  16 

(1) 

.  17 

(2) 

.  18 

(3) 

.23 

(4) 

.23 

(4) 

Density 

200-400 

8067 

.21 

(1) 

.  22 

(2) 

.25 

(4) 

.24 

(3) 

.27 

(5) 

400-800 

4864 

.  29 

(1) 

.29 

(1) 

.31 

(3) 

.30 

(2) 

.30 

(2) 

N2 

120-200 

929 

.  16 

(1) 

.  16 

(1) 

.  19 

(2) 

.21 

(4) 

.20 

(3) 

Density 

200-400 

8247 

.  37 

(1) 

.40 

(2) 

.40 

(2) 

.47 

(3) 

.56 

(4) 

400-800 

3910 

.43 

(1) 

.44 

(2) 

.50 

(3) 

.51 

(4) 

.76 

(5) 

Ar 

120-200 

867 

.  34 

(1) 

.36 

(2) 

.41 

(4) 

.39 

(3) 

.44 

(5) 

Density 

200-400  _ 

3178 

.47 

(1) 

-  5? 

(2) 

.62 

(4) 

.58 

(4) 

.80 

(5) 

1-13 

1-2 

1-1 

1-0 

1-0 

2-1 

2-10 

2-3 

2-1 

2-1 

3-0 

3-1 

3-6 

3-4 

3-2 

4-0 

4-1 

4-4 

4- 8 

5- 1 

4- 2 

5- 9 
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-  Density  and  pressure  were  derived  from  the  hydrostatic  equation 
with  the  aids  of  balloon  radiosonde  at  the  lower  boundary. 
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Abstract.  The  NCAR  thermospheric- general  circulation  model  (TGCM)  is  extended 
to  include  a  self-consistent  aeronomic  scheme  of  the  thermosphere  and  ionosphere.  The 
model  now  calculates  total  temperature,  instead  of  perturbation  temperature  about  some 
specified  global  mean,  global  distributions  of  N(2D  ),  N(4S)  and  NO  ,  and  a  global  iono¬ 
sphere  with  distributions  of  O +,  NO+,  02+  .  AG+  ,  N+  electron  density  and  ion  tem¬ 
perature  as  well  as  the  usual  fields  of  winds,  temperature  and  major  composition. 
Mutual  couplings  between  the  thermospheric  neutral  gas  and  ionospheric  plasma  occur 
at  each  model  time  step.  This  is  the  first  Eulerian  model  of  the  ionosphere  and  it  is 
solved  on  the  TGCM  geographic  grid.  Steady  state  results  for  solar  minimum  equinox 
conditions  are  presented. 


Introduction 

The  development  of  global  numerical  models  of  the  thermosphere  (Fullcr-Rowell  and 
Rees,  1980,  Dickinson  et  ai,  1981)  and  global  numerical  models  of  the  ionosphere  (Sojka 
and  Schunk,  1985)  have  proceeded  generally  independently  of  each  other  Global  empiri¬ 
cal  models  of  the  ionosphere  (e  g.,  Chiu,  1975)  were  used  to  specify  ion  drag  in  thermos¬ 
pheric  general  circulation  models  and  global  empirical  atmospheric  models,  such  as  the 
mass  spectrometer/incoherent  scatter  models  (e  g.,  Hedin,  1987),  were  used  to  specify 
atmospheric  properties  for  ionospheric  models.  Recently.  Fuller-Rowcll  et  til-  (1987)  con¬ 
structed  a  coupled  model  of  the  polar  ionosphere  and  the  global  thermosphere  to  studv 
the  dynamical  and  chemical  interactions  between  the  neutral  and  ionospheric  consti¬ 
tuents  ot  the  polar  regions.  They  adapted  the  Lagrangian  ionospheric  model  of  Quegan 
el  al.  (1982)  to  the  Eulerian  thermospheric  model  of  Fuller-Rowell  and  Rees  (1980)  using 
an  interpolative  scheme  at  each  model  time  step.  Their  calculation  was  successful  in 


reproducing  a  number  of  observed  features  in  the  high  latitude  thermosphere  and 


ionosphere. 

In  this  paper,  we  present  a  self-consistent  coupled  Eulerian  model  of  the  thermosphere 
and  ionosphere  using  an  extension  of  the  aeronomic  scheme  developed  by  Roble  et  al. 
(1987)  and  Roble  and  Ridley  (1987).  The  ionosphere  and  thermosphere  are  both  solved 
on  the  TGCM  geographic  grid  and  includes  consideration  of  displaced  geomagnetic  and 
geographic  poles.  Ion  drift  for  the  ionospheric  calculation  is  obtained  from  the  empirical 
model  of  Richmond  et  al.  (1980)  for  low-  and  mid-latitudes  and  the  empirical  model  of 
Heelis  et  al.  (1982)  for  high  latitudes.  Results  for  solar  minimum  equinox  conditions  are 
presented  that  show  good  agreement  with  MSIS-86  (Hedin,  1987).  The  self-consistent 
model  requires  only  specifications  of  external  sources  as  solar  EUV  and  UV  fluxes,  aurora 
particle  precipitation,  ionospheric  convection  pattern,  and  the  amplitudes  and  phases  of 
semi-diurnal  tides  from  the  lower  atmosphere. 

TGCM  Modifications 

The  NCAR-TGCM  solves  the  primitive  equations  of  dynamic  meteorology  adapted  to 
the  physics  appropriate  to  thermospheric  heights.  The  basic  model  and  developments 
have  been  described  in  detail  by  Dickinson  et  al.  (1981;  1984),  Roble  et  al.  (1982),  Fesen 
et  al.  (1986)  and  Roble  and  Ridley  (1987).  The  TGCM  has  been  used  for  numerous  stu¬ 
dies  involving  ground-based  and  satellite  measurements.  All  of  these  studies  were  per¬ 
formed  with  a  TGCM  that  solved  for  perturbation  temperatures  about  a  global  mean 
temperature  profile  that  was  obtained  from  an  empirical  model  such  as  MSIS  (Dickinson 
et  al.,  1981).  The  previous  procedure  had  been  adopted  because  physical  processes,  such 
as  non-local  thermodynamic  equilibrium  (NLTE)  cooling  emissions  from  CO 2  and  NO, 
and  cooling  by  eddy  transport,  were  not  well  known  at  the  time  of  initial  TGCM 
development.  But'  recently,  Roble  et  al.  (1987)  designed  a  self-consistent  global  mean 
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model  of  the  basic  structure  of  the  thermosphere  and  ionosphere  which  were  successful 
in  calculating  various  thermospheric  and  ionospheric  properties  for  both  solar  minimum 
and  maximum  conditions.  The  major  heating  and  cooling  processes  that  appear  to  he 
responsible  for  maintaining  the  global  mean  structure  were  identified  and  sufficient 
confidence  in  the  aeronomic  scheme  and  adapted  rate  coefficients  was  developed  to  war¬ 
rant  inclusion  in  a  new  version  of  the  TGCM. 

The  thermodynamic  equation  in  the  original  NCAR-TGCM  (equation  (1)  of  Dickinson 
et  al.,  1981)  is  modified  to  calculate  total  temperature, 


dT  _  _ge 
dt  P„ 


PoC, 


d  f  Kt_  dT_ 
dZ  y  H  dZ 


H  ‘  C. 
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-  V  •  VT  -u> 


+(Q  -L)/C, 
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where  T  is  temperature,  t  is  time,  g  gravitational  acceleration,  Cp  is  specific  heat  per 
unit  mass,  Z  is  our  vertical  coordinate  defined  a s  Z  —  In  P0  /P  .  where  P  is  pressure 
and  P0  is  a  reference  pressure  (  5  X  1(1 *  n  b  ),  KT  is  molecular  thermal  conductivity . 
H  is  the  scale  height,  KE  is  the  eddy  thermal  conductivity,  p  is  total  density.  V  is  the 
horizontal  velocity  vector,  w  is  TGCM  “vertical  velocity.”  —  dZ  '  dt  .  R  is  the  imver- 
sal  gas  constant,  m  is  the  mean  mass  and  Q  and  L  are  heating  and  cooling  rates, 
respectively.  The  eddy  diffusion  coefficient  that  is  used  for  these  calculations  is  (Dickin¬ 
son  et  al.,  1981)  KE  (Z )  —  5  X  10"6exp(-7.  Z)s  1  The  eddy  diffusion  coefficient, 
multiplied  by  scale  height  squared  to  be  in  dimensional  units,  has  a  maximum  value  of 
160  m2  s'1  at  the  lower  boundary  (97  km).  The  various  heating  and  cooling  rates  calcu¬ 
lated  by  the  TGCM  are  described  by  Roble  et  al.  (1987). 


The  TGCM  also  has  been  modified  to  include  minor  species  transport  by  solving  addi- 


tional  composition  equations  of  the  form: 
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The  first  three  terms  in  En  represent  gravitational  force,  thermal  diffusion  and  the 
frictional  interaction  with  major  species.  Also  rpn  is  the  mass  mixing  ratio,  An  is  the 
molecular  diffusion  coefficient  and  mn  is  the  mass  of  the  minor  species,  n  ,  fn  is  the 
average  mass  of  the  major  species,  a„  is  a  thermal  diffusion  coefficient,  £  is  a  matrix 
operator  that  represents  the  frictional  interaction  of  a  minor  species  with  the  three 
major  species  O  ,  02  and  N 2,  and  Sn  and  Rn  represents  the  mass  sources  and  sinks  of 
minor  species  n  respectively. 


The  upper  boundary  condition  is  diffusive  equilibrium,  — —  -  En  ipn  =  0  and  the 

dZ 


lower  boundary  condition  can  be  either  photochemical  equilibrium,  a  specified  mass  mix¬ 
ing  ratio,  or  a  specified  mass  flux.  The  odd  nitrogen  chemistry  embedded  in  the  TGCM 
is  the  same  as  that  discussed  by  Roble  et  al.  (1987).  N(~D)  is  assumed  to  be  in  photo¬ 
chemical  equilibrium  throughout  the  thermosphere  and  transport  equations  are  solved 
for  N(*S)  and  NO.  We  consider  N(AS)  in  photochemical  equilibrium  and  specify  a 
constant  number  density  of  8  X  106  for  NO  at  the  lower  boundary. 

The  ion-chemistry  scheme,  rate  coefficients,  and  calculation  procedure  for  the  iono¬ 
sphere  have  been  discussed  by  Roble  et  al.  (1987)  and  Roble  and  Ridley  (1987)  We  con- 


sider  both  magnetic  field  aligned  diffusion  and  E  X  B  transport  in  solving  the  0  + 
transport  equation, 


^  -  Q  +  Ln  =  -  S.  •  nV  (3) 

where  V  —  Vp  +  V±  and, 

£„  =  (£  •  U  -  j-V  Pi)  +  k  -all 

^=w£xi 

n  is  0+  number  density,  t  is  time,  Q  and  Ln  are  production  and  loss  rates  of  0+  and 
Vff  and  are  velocities  parallel  and  perpendicular  to  the  magnetic  field  line,  b_  is  a 
unit  vector  along  the  magnetic  field,  v  is  the  0+  ion-neutral  collision  frequency,  g  is 
gravity,  p{  and  Pj  are  the  ion  mass  density  and  pressure  respectively,  u  is  the  neutral 
wind  vector,  |  B  |  is  magnetic  field  strength  and  £  is  the  electric  field  vector.  Equa¬ 
tion  (3)  is  solved  on  the  TGCM  geographic  latitude  and  longitude  grid  using  the 
mathematical  framework  for  calculating  neutral  dynamics  (Washington  and  Williamson, 
1977).  The  ionosphere  is  thus  solved  using  an  Eulerian,  as  compared  to  the  Lagrangian 
(Sojka  and  Schunk,  1985;  Fuller-Rowell  et  al.,  1987)  approach.  The  upper  boundary  con¬ 
dition  is  nVp  =  <h,  where  is  a  plasma  flux  from  the  magnetosphere  and  is  assumed  to 
be  ^ero  for  these  calculations.  The  lower  boundary  condition  assumes  photochemical 
equilibrium,  n  —  Q/L  Once  the  0+  density  distribution  is  determined,  the  NO +  , 
O  ^  ,  N 2  and  N  +  distributions  are  determined  by  assuming  photochemical  equilibrium 
and  solving  a  quartic  equation  as  discussed  by  Roble  and  Ridely  (1987).  TGCM  tem¬ 
perature,  dynamics,  composition  and  the  global  ionosphere  are  calculated  at  each  time 


step  ( ~  150  s)  and  the  heating  and  cooling  rates  from  chemical  reactions  and  1R 
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emissions  etc.  are  coupled  into  the  thermodynamic  equation  at  the  next  time  step,  as 
described  by  Roble  et  al.  (1987).  Likewise,  neutral  winds,  temperature,  composition, 
and  ionization  rates  are  updated  at  each  time  step  for  the  ionospheric  calculation.  The 
electron  energy  equation  is  not  solved  at  the  present  time  and  the  electron  temperature 
is  arbitrarily  set  equal  to  the  neutral  temperature  for  these  initial  calculations. 

Global  Structure  for  Equinox  Solar  Minimum 

The  new  TGCM  is  run  to  steady  state,  a  diurnally  reproducible  solution,  for  equinox 
conditions  during  solar  cycle  minimum.  The  Hinteregger  (1980)  solar  EUV  and  UV 
measurements  for  July  1976  are  used  to  calculate  photoionization  and  dissociation  rates. 
The  aurora  particle  precipitation  model  of  Foster  et  al.  (1986),  assuming  a  power  input 
of  11  GW,  and  a  Heelis  et  al.  (1982)  ionospheric  convection  pattern,  with  a  cross-polar 
cap  potential  drop  of  60  kV,  are  used  in  the  calculation.  Tidal  amplitudes  are  identical 
to  those  described  by  Fesen  et  al.  (1986).  The  calculated  global  ionospheric  structure  at 
1900  UT  is  shown  in  Figure  1.  The  F-region  peak  density  structure  is  illustrated  by  con¬ 
tours  of  / 0  F 2  plasma  frequency  (MHz)  as  shown  in  Figure  la.  The  prominent  feature  is 
the  density  enhancement  aligned  along  the  magnetic  equator  on  the  dayside  of  the  Earth 
with  a  long  tail  extending  into  the  nightside.  Aurora  particle  precipitation  causes 
enhancements  at  high  latitudes.  Minimum  f  0  F  2  occur  in  mid  latitudes  at  night.  The 
very  low  values  in  that  region  are  attributed  to  F-region  decay  and  the  lack  of  a 
specified  plasma  flow  from  the  magnetosphere  for  this  calculation. 

Contours  of  electron  density  along  the  TGCM  Z  =  -4  constant  possure  surface  near 
120  km  are  shown  in  Figure  lb.  Photochemical  equilibrium  prevails  at  this  altitude  with 
maximum  values  occurring  at  the  subsolar  point  and  in  the  curoral  oval. 


TGCM  calculated  contours  of  neutral  gas  temperate  re  and  wind  vectors  at  1900  UT 
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along  the  Z  —  +2  constant  pressure  surface  near  280  km  are  shown  in  Figure  2a.  Con¬ 
tours  of  the  neutral  gas  temperature  along  the  same  constant  pressure  surface  obtained 
from  MSIS-86  are  shown  in  Figure  2b  for  comparison.  There  is  excellent  agreement 
between  the  TGCM  calculated  temperature  and  MSIS-86  on  the  dayside  of  the  Earth. 
At  night,  the  midnight  temperature  bulge  calculated  in  the  TGCM  appears  too  large, 
compared  with  MSIS-86,  indicating  that  the  tidal  specification  may  be  too  large  or  the 
data  averaging  in  MSIS-86  eliminates  this  variable  feature.  The  neutral  wind  pattern  is 
similar  to  that  calculated  in  previous  studies  with  a  diurnal  component  at  low-  and  mid¬ 
latitudes  and  strong  enhancements  at  high  latitudes  caused  by  ion  convection. 

A  latitudinal  cross-section  of  the  ionospheric  density  structure  at  1900  UT  along  the 
45°  W  longitude  (16  solar  local  time)  is  shown  in  Figure  3.  The  prominent  features  are 
the  equatorial  structure  in  the  F-region  and  auroral  zone  enhancements  superimposed  on 
a  general  dayside  solar  background.  The  low  electron  density  values  in  the  upper  F- 
region  at  mid-latitudes  are  attributed  to  a  zero  flux  plasma  exchange  between  the  iono¬ 
sphere  and  magnetosphere  and  the  low  plasma  temperature. 

A  longitudinal  cross-section  along  the  17.5°  N  latitude  circle  is  shown  in  Figure  1.  In 
the  lower  ionosphere,  photochemical  equilibrium  controls  the  density  distribution. 
Above  about  250  km,  however,  plasma  motion  caused  bv  ionospheric  E  X  B  drifts  and 
neutral  winds  forcing  ionization  up  and  down  magnetic  field  lines  control  the  general 
shape.  In  particular,  the  post-sunset  elevation  of  the  F-layer  is  clearly  shown. 

Finally,  the  TGCM  calculated  nitric  oxide  distribution  at  1900  FT  along  the  Z  ~  1 

constant  pressure  surface  is  shown  in  Figure  5.  The  NO  densities  maximize  at  high  lati¬ 
tudes  with  roughly  a  factor  of  3  difference  between  high  and  low  latitudes.  There  is  also 
a  diurnal  NO  density  variation  in  the  equatorial  region  with  maximum  values  in  the  late 
afternoon.  These  features  are  similar  to  the  Atmosphere  Explorer-C  measurements 
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described  by  Stewart  and  Cravens  (1978)  and  Cravens  and  Stewart  (1978).  There  is  also 
a  magnetic  field  control  at  high  latitudes  associated  with  aurora  production  similar  to 
that  described  by  Cravens  and  Killeen  (1988). 

The  overall  fields  calculated  by  fhe  TGCM  include,  Tn  ,  u  ,  v  .  w  .  h  .  and  mass  mixing 
ratios  for  02,  O  ,  N2,  N('D  ),  N  {'S ),  NO  ,  as  well  as  Of  Of  ,  \0+,  A",  Nf  ,  nf 
and  T,  .  These  fields  are  obtained  at  each  time  step  and  only  a  few  are  presented  in  this 
paper  to  illustrate  the  basic  TGCM  calculated  structure. 

Summary  and  Conclusions 

The  self-consistent  aeronomic  scheme  developed  by  Roble  et  al.  (1987)  has  been  incor¬ 
porated  into  the  NCAR-TGCM  and  used  to  calculate  the  global  structure  of  e  thermo¬ 
sphere  and  ionosphere  for  solar  minimum  equinox  conditions.  The  Eulerian  model  gives 
a  reasonable  representation  of  the  global  ionospheric  structure,  and  it  is  numerically 
efficient  since  thermospheric  and  ionospheric  properties  are  calculated  on  the  same  grid. 
The  TGCM  time  step  is  2  1/2  minutes  for  both  the  thermosphere  and  ionosphere  calcu¬ 
lation  and  it  takes  approximately  1  hour  of  Cray  XMP-18  CPU  time  to  simulate  1  day. 
The  numerical  scheme  is  stable  and  has  been  tested  for  solar  minimum  and  maximum 
with  equinox  and  solstice  conditions  respectively  and  for  aurora  power  inputs  of  100  GW 
and  a  cross-polar  cap  potential  drop  of  150  kV  for  ionospheric  convection.  The  results 
also  indicate  that  a  better  specification  of  magnetosphere-ionosphere  plasma  transport 
including  a  calculation  of  electron  temperature  is  needed  to  improve  the  overall  structure 
of  the  topside  ionosphere  above  the  F-region  peak. 

Concluding  Remarks 


Considerable  progress  has  been  made  in  recent  years  ;n  improving  our  understanding 


of  the  main  physical  processes  operating  in  the  thermosphere  and  ionosphere.  This  pro¬ 
gress  was  possible  because  of  the  availability  of  satellite  data  on  the  structure  and 
dynamics  of  the  thermosphere  and  ionosphere  as  well  as  in  the  development  of  new  glo¬ 
bal  dynamic  models  for  use  in  interpreting  measurements  and  for  predictions  of  global 
phenomena.  It  appears  that  it  may  be  possible,  in  the  near  future,  to  develop  a 
sufficiently  realistic  TGCM  for  use  in  the  prediction  of  thermospheric  weather  that  can 
be  used  to  improve  forecasts  of  satellite  drag  and  at  the  same  time  give  a  prediction  of 
ionospheric  phenomena  for  communication  purposes. 
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EUV  MONITORING  ANO  GEOMAGNETIC  STORM  FORECASTS 
USING  THE  USAF-NOAA  SOLAR  X-RAY  IMAGERS  (SXI) 


W.  J.  Wagner 

Space  Environment  Laboratory,  NOAA/Environmental  Research  Laboratories, 

Boulder,  Colorado  80303 


Abstract 

Neutral  atmosphere  density  models  will  require  continuous  real  tine 
information.  The  histories  and  prognoses  from  data  on  radiant  extreme 
ultraviolet  (EUV)  solar  flux  and  geomagnetic  storm  heating  are  needed  for 
advancing  these  operational  models  beyond  mere  climatology.  The  USAF-NOAA 
Solar  X-Ray  Imagers  (SXI)  will  be  providing  these  data  in  the  1990's. 

Introduction 

The  Sun  is  the  source  of  disturbances  to  the  space  environment.  An 
upgrading  of  the  present  solar  remote  sensing  systems  is  planned  through  the 
implementation  of  X-ray  and  EUV  imaging  of  the  Sun  (Wagner  et  al.  1987).  The 
use  of  Solar  X-Ray  Imagers  (SXI)  will  give  important  improvements  in  the 
prediction  of  proton  flare  events,  geomagnetic  storms,  and  heating  of  the 
Earth's  neutral  atmosphere.  The  SXI  will  permit  better  performance  of  USAF 
missions  by  providing  forecasts  of  disruptions  due  to  these  effects. 
Forecasts  give  USAF  decision-makers  the  planning  time  which  allows  for 
shifting  to  backup  systems  or  for  attaching  more  credibility  to  other  warning 
systems.  The  improved  space  environment  services  resulting  from  the  SXI  will 
also  act  as  a  force  multiplier.  When  forecasts  are  available  for  the  support 
of  advanced  systems,  they  provide  economies  in  resources  such  as  design 
complexity,  launch  weight,  and,  appropriate  to  this  workshop,  satellite  orbit 
1 ifetime. 

The  design  of  the  SXI  copies  that  of  a  telescope  on  NASA's  1973  Skylab 
spacecraft.  Anomalous  orbital  drag  caused  premature  re-entry  of  the  Skylab  in 
1979  before  a  space  shuttle  rescue  mission  could  be  mounted  to  save  this,  the 
nation's  first  space  station.  In  a  sense,  as  we  will  see  below,  SXI  may  be 
considered  to  be  a  spin-off  from  Skylab,  whose  pathfinding  has  returned  to 
serve  against  such  future  calamities.  Already  today,  early  in  solar  activity 
cycle  22,  scientific  spacecraft  are  being  endangered  (Speich  1988)  and 
missions  are  being  delayed  (Industry  Observer  1988)  at  the  same  time  that  DOD 
launch  schedules  (Covault  1988)  are  imperiled  by  the  unpredicted  decay  of 
other  satellite  orbits.  I  will  briefly  discuss  the  SXI  instrument  and  its 
expected  contributions  to  the  neutral  density  atmospheric  models  which  are 
intended  to  mitigate  losses  from  satellite  drag. 
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The  Solar  X-Ray  Imager  and  its  Contributions 

The  Solar  X-Ray  Imagers  will  be  flown  on  the  NOAA  GOES  geostationary 
satellites  at  22000  mile  altitude  with  no  day/night  cycles  (Figure  1).  With 
the  SXI,  grazing  incidence  optics  will  give  continuous  real  time  data  (Figure 
2)  consisting  of  both  extreme  ultraviolet  and  soft  X-ray  images  of  the 
complete  solar  disk.  Images  will  arrive  once  per  ten  minutes  (at  fastest 
rates,  up  to  one  per  minute)  and  will  have  pixel  size  of  5  arc-seconds.  The 
passbands  to  be  used  include  10  -  20  and  20  -  60  Angstroms  in  the  soft  X-ray 
region,  and  255  -  300  Angstroms  in  the  EUV. 

To  select  the  optimal  EUV  band  for  monitoring  the  radiant  energy  input  to 
the  thermosphere,  previous  investigations  by  AFGL  were  used  (Hinteregger 
1981).  Additional  studies  of  the  AFGL  data  (Hedin  1984;  Hedin  and  Mayr  1987; 
Roble  1987)  seem  to  indicate  that  solar  flux  in  the  255  -  300  Angstrom  band 
best  correlates  with  the  density  of  thermospheric  nitrogen  (Figure  3).  A  heat 
rejection  prefilter  covers  the  aperture  of  the  SXI  thus  preventing  passage  of 
wavelengths  above  approximately  350  Angstroms.  Unclear  at  this  time  is 
whether  to  include  the  strong  but  relatively  constant  contribution  of  He  II 
304  Angstroms.  However,  it  appears  (Figure  3,  from  Hedin  1984)  that  this 
spectral  line  would  serve  only  to  degrade  the  data  if  the  intention  were  to 
monitor  the  changing  thermospheric  densities. 

The  SXI  will  measure  the  EUV  flux  over  an  entire  solar  cycle.  With  a 
series  of  SXIs  installed  on  the  GOES  satellites,  it  is  expected  that  long-term 
variations  of  these  255  -  300  Angstrom  EUV  fluxes  will  be  available  for 
comparison  with  other  longer-established  solar  indices  (Figure  4).  Although, 
no  overall  absolute  radiometric  calibration  will  be  provided  in  orbit,  any 
sensitivity  changes  are  expected  to  be  gradual  and  uniform  with  time  and  thus 
amenable  to  periodic  checks  by  other  space  instruments.  The  SXI  detector,  a 
charge-coupled  device,  will  be  regularly  calibrated  onboard  in  a  visible 
wavelength  band. 

In  addition  to  monitoring  the  solar  EUV  radiant  flux  for  input  into 
thermospheric  density  models,  the  SXI  will  also  aid  in  the  prediction  of  the 
heating  of  the  Earth's  atmosphere  that  is  caused  by  geomagnetic  storms.  The 
location  of  X-ray  coronal  holes  and  the  details  of  their  evolution  and  shape 
will  be  used  to  infer  high  speed  streams  in  the  solar  wind  structure.  Long- 
duration  X-ray  events,  believed  to  mark  the  launch  of  coronal  mass  ejections, 
will  also  be  detected  (by  the  soft  X-ray  channel)  as  they  leave  the  Sun  . 

A  Memorandum  of  Agreement  was  signed  by  the  USAF  and  NOAA  in  June  1987 
providing  funds  M  October  1989.  At  the  same  time,  NASA  declined  to  join  the 
.  -"■  GOES  Munches  are  now  posted  for  May  1992,  February  1997,  and 
!■  e  r*9\  f  funding  earlier  than  Fy  °C  can  be  identified,  the  1992 
'  *  -u  '  r  r-y  the  fir'-.  JXI  into  oroit. 


Conclusions 


The  SXI  program  will  deliver  10  to  15  years  of  data  and  prediction 
services  through  the  USAF  and  NOAA  forecast  centers  in  Colorado  (Mulligan  and 
Wagner  1987).  Accurate  forecasts  of  satellite  drag  can  save  the  nation  real 
money.  They  will  permit  the  coherent  and  consistent  space  shuttle  mission 
planning  which  was  advocated  by  the  Rogers  Commission  on  the  Challenger 
accident.  Today,  circumstances  still  resemble  the  situation  in  1979  at  the 
time  of  the  inadvertant  Skylab  re-entry.  However,  with  workshops  such  as  this 
one  and  using  neutral  atmosphere  density  models  supplied  with  on-going  real 
time  EUV  input  data  and  with  geomagnetic  storm  forecasts  from  Solar  X-ray 
Imager,  we  are  likely  to  minimize  losses  due  to  unanticipated  orbital  decay  in 
the  future. 
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Figure  Captions 


Figure  1.  Real  time  SXI  EUV  and  X-ray  data  will  be  delivered  directly  to  Air 
Force  and  NOAA  forecast  centers  for  dissemination  to  users. 

Figure  2.  X-ray  image  of  the  Sun  from  the  NASA  Skylab  mission.  Data  from 

the  SXI  will  resemble  this.  Currently,  no  image  exists  which  was 
taken  in  the  255  -  300  Angstrom  band  designated  for  the  SXI. 

Figure  3.  This  work  by  Hedin  (1984)  shows  that  nitrogen  densities  (error  bars 

in  the  top  panel)  correlate  better  with  255  -  300  Angstrom  fluxes 
(second  panel)  than  they  do  with  (a  model  driven  by)  10-cm  radio 
fluxes  which  are  shown  as  the  solid  curve  in  the  top  panel.  The 
inclusion  of  He  II  304  Angstroms  (third  panel)  within  this  EUV  band 
would  actually  lower  the  correlation  with  atmospheric  densities. 

Figure  4.  Solar  cycle  variations  of  various  indices  that  have  been  used  as 
surrogates  for  the  heating  by  EUV  flux  (according  to  Donnelly 
1988).  A  comparison  of  their  dissimilar  solar  cycle  variations 
demonstrates  the  differences  among  these  surrogates. 
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Abstract 


The  Remote  Atmospheric  and  Ionospheric  Detection  System  (RAIDS)  experiment,  to  fly  on  a 
TIROS  spacecraft  in  the  late  1980's,  consists  of  a  comprehensive  set  of  one  limb  imaging 
and  seven  limb  scanning  optical  sensors.  These  eight  instruments  span  the  spectral  range 
from  the  extreme  ultraviolet  to  the  near  infrared,  allowing  simultaneous  observations  of 
the  neutral  and  ion  composition  on  the  day  and  night  side  as  well  as  in  the  auroral 
region.  The  primary  objective  of  RAIDS  is  to  demonstrate  a  system  for  remote  sensing  of 
the  ionosphere  to  produce  global  maps  of  the  electron  density,  peak  altitude  and  critical 
frequency . 


Introduct ion 


The  Remote  Atmospheric  and  Ionospheric  Detection  System  (RAIDS)  is  a  satellite 
experiment  designed  to  perform  a  comprehensive  study  of  the  upper  atmospheric  airglow 
emissions.  RAIDS  observations  will  be  used  to  develop  and  test  techniques  for  remote 
sensing  of  the  neutral  atmosphere  and  ionosphere  on  a  global  scale.  The  RAIDS  experiment 
consists  of  a  set  of  one  limb- imaging  and  seven  limb-scanning  optical  sensors  which  sample 
the  altitude  interval  75  -  750  km  and  cover  the  wavelength  range  500  -  8700  A. 

The  goal  of  the  RAIDS  experiment  is  to  obtain  a  set  of  simultaneous  airglow  profiles  at 
a  number  of  wavelengths  which  will  be  used  to  develop  and  evaluate  techniques  for  neutral 
atmospheric  and  ionospheric  remote  sensing.  The  RAIDS  instrumentation  will  acquire  a 
global  database  of  airglow  intensities  which  will  be  used  in  conjunction  with,  and 
compared  to,  theoretical  models  of  radiation  transport,  photochemistry  and  dynamics  to 
examine  in  detail  the  relationships  between  atmospheric  composition  and  airglow.  The 
primary  focus  of  RAIDS  will  be  on  remote  sensing  of  the  ionosphere  since  there  is 
considerable  interest  by  the  ionospheric  and  high  frequency  propagation  communities  in 
monitoring  the  ionosphere  in  real-time  on  a  global  basis. 

RAIDS  has  been  proposed  to  the  U.S.  Air  Force  Space  Test  Program  and  is  currently  being 
considered  for  flight  on  a  NOAA  TIROS-N  series  satellite.  The  planned  orbit  will  be 
sun-synchronous  in  the  early  afternoon  at  an  altitude  of  870  km. 

Experiment  Objectives 

The  ionospheric  electron  density  on  the  dayside  will  be  determined  from  measurements  of 
the  extreme  ultraviole.  .d'JV)  emissions  ionized  oxygen  (primarily  at  834  A)  in  the  F2 

region  of  the  ionosphere1'  .  The  ionized  oxygen  dayglow  is  produced  from  lonization- 
excitation  if  neutral  atomic  oxygen  and,  as  the  EUV  radiation  undergoes  multiple 
resonant  scattering,  the  radiation  fjeld  acquires  the  signature  of  the  O  ion 
distribution.  Measurements  of  the  0  EUV  radiation  by  limb  scanning  instruments  can  be 
used  to  infer  the  electron  density  profile  throughout  most  of  the  F2  region  where  the  0* 
density  is  largely  equal  to  the  electron  density. 

The  nightside  electron  density  will  be  determined  by+measurement  of  the  nightglow  from 
neutral  oxygen  atoms  produced  by  the  recombination  of  0  and  0,  ions  and  electrons. 

These  emissions  have  been  observed  at  a  number  of  wavelengths  Including  the  EUV  911  8 
continuum,  the  far  ultraviolet  (FUV).1304  S,  1356  A  multiplets,  and  in  the  visible  and 
near  infrared  at  6300  X  and  7774  ' 5 ' 6 ' 7 .  scanning  through  the  limb  gives  the  shape  of 

the  emission  profile  which  can  be  used  to  determine  the  altitude  distribution  of  the 
electron  density.  The  magnitude  of  the  nightglow  recombination  emission  has  been  shown  to 
be  proportional  to  the  square  of  the  peak  electron  density  in  the  F2  region*'  . 
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Another  scheme  for  ior.osc^er i  -  remote  sensing  to  he  explored  by  RAIDS 
experiment  involves  the  s  im„i tuner.-s  measurement  of  dayglow  emissions  from  neutral  oxygen 
and  molecular  nitrogen  in  the  ..art  of  the  spectrum.  The  atomic  jxygen/nolec  jl  ar 

nitrogen  ratio  has  been  shown  to  :e  related  to  tne  distribution  of  0  ions’ . 

In  addition  to  the  EL'V  and  FL'V  sensors,  the  RAIDS  experiment  will  contain  sensors  to 
measure  airglow  emissions  in  the  middle  and  rear  ultraviolet  ;  M'JV,  NUV)  and  the  near 
infrared  (NIR)  which  will  oe  used  to  determine  the  neutral  atmospheric  composition  (0,  N0 
and  Oj),  temperature,  energetic  electro^  flux  and  the  densities  of  .minor  atmospheric  ‘ 
speciSs  such  as  NO,  N,  Ma ,  H,  He  and  Mg  .  The  magor  and  minor  neutral  species  will 
provide  important  inputs  into  pnotcchem i ? a  1  and  dynamical  models  of  the  thermosphere  and 
ionosphere.  By  making  airglow  observations  over  an  extended  spectral  range,  the  RAIDS 
experiment  will  provide  enouoh.  mbits  to  adequately  test  the  theoretical  models.  Table  1 
summarizes  the  magor  scientific  chgectives  of  RAIDS  and  the  emission  features  that  will  be 
observed  to  satisfy  the  objectives. 

TAB LL  1.  RAIDS  SCIENTIFIC  OBJECTIVES 


SCIENCE  OBJECTIVES 

METHOD 

WAVELENGTHS  (X) 

Dayside  Electron  Density 

0  Day glow 

540-834  (0+) 

Nightside  Electron 

0  Recombination 

911,  1304,  1356,  7774  (0) 

Density 

0,  Recombination 

6100  (0) 

Neutral  Composition 

N .  Airglow 

1250-2000  (N-  LBH 1 

And  Temperature 

0  Airglow 

989,  1304,  1356,  1641  (C) 

j  Absorption 

l’S0-53S0  (0^  Schunann-Runge ! 
7400-10000  !0,  Atmospheric' 

0*0  Association 

2500-5000  (0.,  Herzberg  I) 

Energetic  Electrons 

N . ,  N.  A 1 n low 

1 

1250-20 CO  (N,  LBH' 

3100-4000  !N.  Second  Positive 
i 

N j  First  Negat ive ) 

Minor  Species 

NO 

NO  Airoiow 

2000-2  SCO  (NO  v  1 

N 

NO 

1900-2400  (NO  f  i 

Na 

Na  ‘ 

5890  ;  Na  "D") 

Mg 

M-rn  * 

2 8  CO  (Mg  ) 

H 

u  * 

1027  (H  Ly  a  ) 

He 

He 

564  ( He  ) 

Observing  Scheme 

Figure  1  is  an  artist's  cent  - ;  - 1  or.  of  the  RAIDS  experiment  ih  orbit  aboar  ‘  =  ~I7.-'D 

satellite.  The  experiment  w,'. '.  be  mounted  to  the  underside  of  the  spacecraft  a "  tve 

instruments  will  view  albn  i  th.  an  1 1  -  1  oc  1 1  y  direction.  The  instrument  r  - *  i  -  .3 

divided  into  three  sec  : :  a  fixed  central  *T"  Shaped  structure  housm.  ~  '■  '  '  " 

mechanism  and  the  1  i-b~  mag  mg  .instrument ,  and  two  rotatable  side  compart -cm.',  s  -  1 .  c  ‘  '-.c  ;S 
the  limb-scanning  instruments.  The  primary  instrument  package  has  dmens:""  '  :  I .  i  m 
(X)  by  63.5  cm  (Yl  bv  83.8  cm  ' 1 i  where  X  is  n  the  nadir  direction,  i  is  ocuU'  . -  •  ;e 

velocity  vector  and  Z  is  along  the  long  axis  of  the  satellite. 
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scanning  gratings  and  single  photomultiplier  detectors,  and  (3)  single  channel 
photometers.  The  spectrographs  record  emissions  from  several  wavelengths  simultaneously 
over  a  fixed  spectral  interval  while  the  spectrometers  can  be  commanded  either  to  scan 
across  wavelength  range,  obtaining  data  at  each  wavelength  m  series  or  stop  and  record 
data  at  a  single  wavelength.  The  photometers  have  fixed  ban-passes  which  isolate  a  single 
airglow  emission  feature.  A  summary  of  some  of  the  more  important  instrument  parameters 
is  contained  in  Table  2. 


TABLE  2.  RAIDS  INSTRUMENT  PARAMETERS 


Instrument  Wavelength  Bandpass 
Range 


Field* 

of 

Vie  v 


Integration# 

Time 


Aperture 

Area 


Peak 
Sens . 


Telemetry 

Rate 


<X, 

(X) 

(deg2 ) 

EUV 

Spectrograph 

500-1100 

12 

0. 1x2.4 

FUV 

Spectrograph 

1300-1700 

8 

4.0x0. 1 

<1UV 

Spectrometer 

1900-3200 

10 

0. 1x2. 1 

NUV 

Spectrometer 

2950-4000 

7 

0. 1x2. 1 

NIR 

Spectrometer 

7250-8700 

7 

0. 1x2.1 

5890 

Photometer 

5890 

20 

0. 2x2.  1 

6300 

Photometer 

6300 

20 

0. 2x2.1 

7774 

Photometer 

7774 

20 

0. 1x2.  1 

(sec) 

( mm2  ) 

(C/s/R) 

(Bits/s) 

0.5 

70x70 

0.5 

2560 

87.0 

25x25 

2.5 

4096 

0.1  or  . 

025 

42x50 

7.  5 

800 

0.1  or  . 

025 

21x25 

1.8 

800 

O 

b- 

0 

025 

21x25 

2.8 

800 

0.1 

50x50 

6.  2 

200 

0.1 

50x50 

8.0 

200 

0.1 

25x25 

4.6 

200 

*  The  first  and  second  values  of  the  f ield-of-v lew  correspond  to  the  vertical  and 
horizontal  dimensions  respectively. 


@  The  integration  times  of  the  MUV,  NUV,  and  NIR  spectrometers  are  given  for  both 
the  f  ixed-waveleng tn/al 1 1 tude-scan  and  the  f ixed-al 1 1 1 ude/waveleng t h-scan  modes 
respect i vel y . 


El'V  and  FL’V  Spectrographs 

The  ELV  instrument  is  an  f/S,  Wadsworth  concave  grating  spectrograph  (Figure  2). 

Inc. n;m  l.gnt  enters  through  a  mechanical  modulation  collimator  and  is  focussed  on  a 
.sensitive  MC?  array  detector  by  a  spherical  reflectance  grating  used  m  first 
,  -j  r.  r  -  ,  ►  .  r.  ■,  nas  ,3  a  '2  x  "2  mm,  a  ruling  Jensity  of  l".  0  line  s/mm  with  a 

:.:s  irdet  o'.azo  at  700  8  and  a  355  mm  focal  length.  The  grating  can  be  rotated  mt 
two  txasif.  ons  focussing,  350  8  jegments  between  500  and  1000  A  on  the  detector.  T 
two  segment-:  overlap  so  that  the  0  834  A  feature  is  measured  m  either  grating  position. 

;-.e  spectra  1  msoijtion  is  12  A.  The  field  of  view  defined  by  the  collimator  and  grating 
. x  2.1  neg.  The  instrument  is  mounted  on  the  scan  platform  so  that  the  long  axis  of 
--a  fw>ii  ot  view  is  parallel  to  the  horizon  on  orbit. 
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Figure  2.  Layout  diagram  of  the  RAIDS  EUV  spectrograph. 

The  EUV  detector  consists  of  a  stack  of  three  MCP's  in  a  *Z"  stack  configuration  with  a 
wedg e-a nd— s tr 1 p  type  anode.  The  MCP's  are  36  mm  diameter  rimless  plates  with  a  27  mm 

active  area  and  80:1  leng th/d lameter  channel  ratio.  The  detector  is  windowless  and  a  Mr,  p 

phJtacatnode  is  deposited  directly  on  the  top  MCP.  The  detector  can  resolve  256  x  256  2 
pixels  but  on  orbit  it  will  be  operated  using  128  x  1  resolution  elements  aligned  along 
the  wavelength  dispersion  direction.  The  detector  is  enclosed  in  a  vacuum  housing  which 
is  sealed  with  a  DC  powered  cover  assembly  which  will  be  opened  m  flight. 

Tho  design  of  the  FUV  spectrograph  uses  an  f/S  125  mm  focal  length  imaging  Wadsworth 
configuration.  The  optical  elements  in  the  spectrograph  consist  of  a  spherical  telescope 
mirror  used  off-axis,  a  collimating  mirror,  a  concave  grating  and  a  MCP  array  detector 
(Fioure  3).  The  two  mirrors  and  the  grating  have  focal  lengths  of  125  mm.  The  0.25  mm  x 

9.0  mm  entrance  slit  is  at  the  focus  of  the  telescope  mirror  defines  a  field  of  view  of 

0.1  x  4.0  deg.  The  grating  has  a  ruled  area  31  mm  x  31  mm,  a  ruling  density  of  3600 
lines/mm  and  is  blared  at  1400  ft.  The  FUV  spectrograph  contains  no  moving  parts  and  it 
win  be  .lxed  to  the  mounting  surface  of  the  spacecraft.  The  long  axis  of  the  entrance 
slit  wi  l  be  perpendicular  to  the  horizon  during  flight  and  the  spectrograph  will  image 
vertically  along  the  length  of  the  slit.  The  wavelength  dispersion  direction  is 
perpend  ;c ui. a r  to  the  slit.  The  wavelength  coverage  is  1300  -  1700  8  and  the  spectral 
resolution  is  6-8  X.  r 
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FAR  ULTRAVIOLET  IMAGING  SPECTROGRAPH 
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Figure  3.  Layout  diagran  of  the  RAIDS  Ft'V  imaging  spectrograph. 

The  FUV  spectrograph  detector  consists  of  a  per~anentiy  sealed  evacuated  housing 
enclosing  a  stack  of  four  MCP'n  and  a  wedge  and  strip  anode  array.  The  detector  window  is 
MgF_  and  is  coated  on  the  inside  with  a  CsX  photocathode.'  The  MCP’s  are  25  mm  in  diameter 
and  nave  a  lengtn/diameter  channel  ratio  of  40:1.  The  detector  will  be  operated  in  a  two 
dimensional  mode  using  256  pixels  along  the  dispersion  direction  and  64  pixels  al_oq  the 
altitude  direction.  The  altitude  resolution  is  3.5  km  and  the  minimum  ray  height  ooint 
altitude  range  imaged  is  75  kn  to  285  km. 

MU V ,  Nil V  and  N IR  Spectrometers 

The  MUV,  NUV  and  NIR  spectrometers  consist  of  f/5,  125  mm  focal  length  Ebert-Fastie 

monochromators  with  off-axis  telescopes  containing  spherical  mirrors.  The  grating 
positions  are  determined  by  stepper-driven  cam  mechanisms  vncn  can  be  controlled  by  the 
flight  microprocessor  either  to  scan  continuously  m  wavelemth.  or  to  step  to  a  particular 
wavelength  position  and  record  data  at  that  fixed  wavelength.  The  grating  cams  and  aear 
u'.ios  were  closer,  ,  .  =.  • ;  .  e  t  .•  scan  of  earn  spec  .  *eter  consists  nf  1'5"4  motor 

steps  and  takes  approximately  25  seconds.  The  entrance  arcl  exit  silts  of  the 
monochromators  are  matched  and  fields  of  view  of  the  three  spectrometers  are  0.1  x  2.1 

deg.  As  with  the  EL'V  spectrograph  .  tr.e  spectrometers  are  to  me  meun.-ed  to  the  scan 
platform  with  the  long  axis  of  the  slit  parallel  to  the  tor : ton  in.  flioht. 

The  telescope  of  the  MUV  spectrometer  is  f/5  with  a, 250  mm  focal  length .  V  t  r- 
instrument  uses  a  3600  iine/mm  grating  blazed  at  240C  K.  ■nth  an  entrance  slit  sic;  ~.c 
0.5  mm  x  9.0  mm  the  spectrometer  -as  a  resolution  of  9.6  A.  The  spectral  cover  a  a z  ..  the 
MUV  spectrometer  extends  from  »  9LC~  3iCC  K.  The  .MUV  s  pec  t  rem  e  te  r  detect. r  a  2C  ■“ 

diameter  photomultiplier  with  a  CsTe  pnotcca  th.cde  ana  MoF.  window  and  ...  -  •  "d  at  t  h . 

exit  slit  of  the  monochr  .r.ator . 

The  MUV  spectrometer  nas  an  f,  5 ,  ...  mr.  te.es  o  .u  a  no  a  1.25  :.ht.  -/  •  - 
grating  has  2400  lines /mm  and  is  blazed  at  3CQ0  A.  The  spectral  cover. me  • 

2950-4000  A  and  the  spectral  resolution  is  7  X.  The  NUV  spectrometer  .  2 

diameter  photomultiplier  with  a  oi-aikaii  photocathode  and  aiass  wmdc-  Fiiv.;  :  -  \  . 
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Figure  4.  Layout  diagram  of  the  NUV  spectrometer. 

The  telescope  and  slit  dimensions  of  the  NIR  spectrometer  ar^  identical  to  the  NUV 
instrument.  The  NIR  grating  has  1800  line/mm  and  is  blazed  at  7600  a.  The  wavelength 
scan  range  is  7400-8700  A  and  the  spectral  resolution  is  7  2.  Tho  detector  is  a 
photomultiplier  with  a  GaAs  photocathode  and  glass  window.  Because^detector  must  be 
cooled  to  reduce  the  dark  count  to  an  acceptable  level,  the  photomultiplier  tube^TSTs 
housed  in  a  separate  box  and  is  radiatively  cooled  to  approximately  -20  C.  A  3  nm 
J: anetor  fiber  optical  bundle  is  used  c-  deliver  the  light  from  the  monochromator  to  the 
detector.  A  lens  assembly  couples  light  from  the  monochromator  into  the  ‘iber  bundle. 

Photometers 

The  RAIDS  experiment  includes  three  photometers  with  bandpasses  centered  at  the  sodium. 
5890  A  line,  and  the  oxygen  6300  8  and  7774  8  lines.  The  photometers  use  off-axis 
telescopes  identical  to  those  on  the  spectrometers,  and  field  stops  at  the  foci  of  the 
telescopes  limit  the  f leld -of-view  to  0 . 1  x  2.1  deg.  At  the  field  stop  an  optical 
coupling  assembly  collects  light,  renders  it  parallel,  passes  it  through  an  interference 
filter  and  refocusses  it  into  a  3  mm  diameter  fiber  optical  bundle.  The  bundle  transnits 
the  light  to  the  radiatively  cooled  (separate)  detector  box.  Each  photometer  uses  a  GaAs 
photomultiplier  as  a  detector. 

In  addition  to  providing  information  on  the  distribution  of  the  Na  D  line  airglow,  the 
5890  A  photometer  will  be  used  as  an  independent  absolute  calibration  of  the  scan  platform 
angle.  The  height  of  the  peak  of  the  Na  airglow  is  90  ♦  5  km,  and  limb  scans  of  this 

2  36 

SPl£  Voi  M7  UPrtviolH  T tchnologr  O  966) 


emission  feature  will  provide  a  check  on  the  accuracy  of  the  instrument  look 
direction.  The  Na  airglow  will  be  measured  on  both  the  dayside  and  mghtside. 


Flight  Microprocessor  and  Scan  Platform 

The  RAIDS  flight  microprocessor  (FMP),  which  is  based  on  a  Harris  HS-8085RH 
microprocessor,  functions  primarily  as  a  spacecraft  command  interpreter  and  data  manager 
for  telemetry.  The  FMP  accepts  8-bit  coded  commands  and  executes  the  associated  stored 
command  macros  after  checking  for  conflicts.  The  FMP  manipulates  the  large  quantity  of 
data  received  from  the  instruments  and  multiplexes  it  into  the  satellite  telemetry  stream. 
Other  operations  performed  by  the  FMP  include  power  sequencing  for  the  instruments  and 
heaters,  setting  the  wavelength  positions  of  the  spectrometers,  controlling  the  scan 
platform,  performing  analog-to-digital  conversions,  and  controlling  data  buffers  for  the 
instruments . 

The  electronics  have  been  designed  to  minimize  the  microprocessor  intervention  and  to 
prevent  over-tasking.  Most  subsystems  are  autonomous  except  for  inputs  such  as 
wavelength,  power  on/off  commands,  and  cross-strapping  the  various  data  channels.  Analog 
to  digital  conversions,  discrete  digital  status,  and  the  positions  of  the  various  stepper 
motors  must  be  read  and  stored  by  the  microprocessor. 

The  scan  platform  is  the  principal  structural  unit  of  the  RAIDS  instrument.  All  the 
instruments  except  the  FUV  spectrograph  are  mounted  to  the  movable  section  of  the 
platform.  The  mechanically  scanned  instruments  are  located  on  one  of  four  vertical  fins. 
There  are  two  sets  of  fins  on  each  side  of  the  platform  along  with  the  stationary  PUV 
spectrograph  in  the  center.  The  fins  are  coupled  in  pairs  into  two  integral  scanning 
assemblies  that  are  tilted  independently  by  a  pinion  and  sector  gear  drive.  The  scanning 
assemblies  rotate  about  a  point.  ftJW  to  the  chassis  base. 

A  single  stepper  motor  rotates  the  scanning  assemblies.  The  direction  of  the  motor 
rotation  is  determined  by  the  state  of  microswitches  accuated  by  the  platform  at  the  top 
and  bottom  of  the  scan  cycle.  A  counter  is  used  to  indicate  the  position  of  the  platform 
and  is  zeroed  after  each  cycle.  A  linear  potentiometer  is  also  attached  to  the  platform 
to  monitor  the  platform  motion. 
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FROM  DMSP-F4  UV  DATA 


RAIDS  DATABASE 


•  NOAA  data  tapes 
stored 


•  Arrival  and  timespan  (rev*, 
date,  etc.)  of  NOAA  tapes* 
cataloged  on  a  PC  AT  usino 
Rbase.  NOAA  tapes 
converted  from  IBM  format 
to  VAX  format. 

•  Data  blocks: 
Reconstructed 
time-ordered  data 

•  Averaged,  calibrated  data 
(intensity  units) 


•  Maps  of  observed  intensities  ft# 
selected  features:  auroral 

and  airglow  backgrounds 

•  ZA'may  be  further  condensed 
to  minimize  storage  soace 

or  act  as  a  browse  file 

•  Emission  rates, 
density  profiles, 
atmospheric  composition 


•  Scientific  analysis 


•  "Final"  product  as  maps 
TEC,  f„F2,  Hmixt  nmgx  , 

electron,  ion  and 


neutral  temperatures. 
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BEGIN  TO  CREATE  AN  UPPER  ATMOSPHERIC  DATA  BASE  INCLUDING 
MEASURED  INPUTS  OF  THE  MAJOR  ENERGY  SOURCES 
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CHAIRMAN 


SUMMARY  OF  THE  NEUTRAL  DENSITY  CONFERENCE 
by  Dr  H.  Beat  Wackernagel 

Hq  Space  Command/Di rectorate  of  Operations  Analysis 

This  afternoon  Lieutenant  Colonel  George  Davenport  and  I  will 
summarize  what  we  have  heard  at  this  conference.  Colonel 
Davenports  summary  will  be  from  the  perspective  of  the  purveyor 
of  atmospheric  data,  while  I  will  deal  with  the  same  matter  from 
the  perspective  of  the  consumer  (or  user)  of  atmospheric  data. 

In  addition  there  will  be  a  recapitulation  of  action-items  as  we 
heard  them  -  corrections  from  the  floor  are  welcome! 

The  user  has  the  task  of  calculating  satellite  orbits,  both 
in  terms  of  determining  the  orbit  on  the  basis  of  observations  and 
environmental  data  (establishing  past  positions),  as  well  as  the 
task  of  calculating  predictions  of  future  positions  of  the 
satellite.  For  doing  this  the  user  needs  those  data  which  are 
called  for  by  the  model  (or  models)  for  the  caluclation  of  both, 
the  past,  and  future  satellite  positions. 

Computational  Efficiency. 

In  the  presentation  by  Capt  de  la  Barre  computational 
efficiency  was  mentioned  as  one  of  the  requirements.  Since  this 
aspect  is  sometimes  not  fully  appreciated  by  the  model-building 
community,  I  will  take  this  occasion  to  elaborate  and  to  provide 
some  numbers. 

The  total  catalog  contains  some  19,000  satellites  of  which 
7,300  are  in  orbit  today,  while  the  others  have  decayed.  Hence, 
the  number  of  satellites  the  computing  center  has  to  handle  -  on 
a  real-time  basis  -  is  a  thousandfold  more  than  the  handful  of 
satellites  used  in  a  research  project. 

The  data  influx  is  presently  in  the  vicinity  of  50,000  to 
60,000  observations  per  day,  seven  days  per  week.  For  planning 
purposes  use  100,000  per  day,  which  means  an  average  of  about  one 
observation  per  second  -  24  hours  a  day.  In  a  realtime  system 
you  may  not  fall  behind  -  the  data  must  be  processed  as  they  flow 
into  the  system.  In  a  favorable  case  each  observation  is 
processed  twice:  a  first  time  to  verify  the  identification,  so 
that  it  can  be  placed  into  the  correct  data  file  and  a  second  time 
for  calculating  an  improved  orbit.  In  the  case  of  an  untagged  or 
mistagged  observation,  the  computational  effort  is  vastly  higher. 
Because  of  computational  limitations  only  about  one  or  two 
percent  of  the  satellites  are  maintained  with  Special 
Perturbations  (including  a  full  atmospheric  density  model). 

Density  information  is  required  at  each  integration  step  for  each 
iteration  in  the  improvement  process  for  each  satellite. 
Conjectures  for  such  numbers  are: 


2000  x  5  x  100 


a  million. 


Because  many  other  activities  must  go  on  concurrently  in  a 
realtime  system,  only  a  small  fraction,  perhaps  two  percent,  of 
the  computing  resources  can  be  devoted  to  neutral  atmospheric 
density  calculations;  that  comes  to  perhaps  2000  seconds  per  day 
during  which  a  million  requests  must  be  serviced;  that  divides 
out  to  two  milliseconds  per  request.  If  the  entire  catalog  is  to 
be  serviced,  that  number  is  hundred  times  smaller,  namely  only 
20  microseconds. 

Hence,  in  addition  to  better  accuracy,  the  software  of  the 
next  model  must  be  computationally  efficient  in  the  way  it  is 
operationally  implemented.  An  available  option  might  be  to 
evaluate  the  model  once  every  few  hours  and  store  the  results  in 
tables  from  which  the  data  can  be  retrieved  more  repidly;  however 
there  are  important  limitations  as  to  the  size  of  available 
storage  areas. 

Additional,  data  retrieval  algorithms  assessing  large  data  bases 
require  non-trival  amounts  of  machine  time.  This  is  clearly  an 
area  where  some  trade-offs  have  to  be  studied  quite  carefully. 
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The  existing  procedures  for  specifying  and  predicting  the 
atmospheric  neutral  density  at  satellite  altitudes  leave  a  great 
deal  to  be  desired.  Air  Weather  Service  just  reports  the 
observations  that  describe  the  environment,  and  the  various 
customers  enter  these  input  data  into  their  own  atmospheric  neutral 
density  models.  Each  customer  uses  these  observed  data  to  execute 
the  model  that  has  been  installed  on  its  computer  system;  the  model 
may  or  may  not  generate  output  that  really  satisfies  the  customer's 
needs . 

Air  Weather  Service  should  be  providing  customers  with  the 
output  from  the  best  possible  atmospheric  neutral  density  model 
rather  than  input  observations.  Rather  than  merely  relaying 
observed  environmental  parameters  for  the  customers  to  input  into 
their  own  models,  the  appropriate  procedure  would  be  for  Air 
Weather  Service  to  provide  the  various  customers  with  the 
parameters  that  they  really  need,  the  actual  specification  or 
prediction  of  the  atmospheric  neutral  density  at  satellite 
altitudes . 

For  Air  Weather  Service  to  assume  responsibility  for  executing 
the  optimum  atmospheric  neutral  density  model  will  require  a 
fundamental  change  in  the  way  that  Air  Weather  Service  and  the 
customer,  in  this  case  Air  Force  Space  Command,  conduct  business. 
Currently  Air  Weather  Service  forwards  to  Air  Force  Space  Command  a 
standardized  message,  including  the  observed  and  predicted  ?800 
MHz,  10.7  centimeter  wavelength  solar  radio  flux  (the  Ottawa  Flux 
or  Fio)^  the  90-day  mean  F]0,  and  the  observed  and  predicted 
geomagnetic  index  Ap;  Air  Force  Space  Command  extracts 
appropriate  parameters  from  the  Air  Weather  Service  message  for 
input  into  the  Jaccbia  6^  atmospheric  neutral  density  model 
resident  on  the  computer  system  inside  Cheyenne  Mountain;  finally 
Air  Force  Space  Command  uses  the  output  from  the  Jacchia  6^  model 
to  carry  out  drag  computations.  In  the  future.  Air  Weather  Service 
will  accumulate  all  of  the  best  possible  observed  and  forecast 
environmental  parameters;  these  data  will  be  input  into  a 
state-of-the-art  atmospheric  neutral  density  model  that  Air  Weather 
Service  will  execute  at  the  new'  Space  Forecast  Center,  which  is  to 
be  established  at  Falcon  AFS;  Air  Weather  Service  will  transmit 
relevant  model  output  parameters  from  the  Space  Forecast  Center  to 
Air  Force  Space  Command  in  Cheyenne  Mountain;  and  Air  Force  Space 
Command  will  use  the  atmospheric  neutral  density  output  from  the 
Air  Weather  Service  model  to  compute  the  drag  on  its  highest 
interest  satellites. 

Air  Force  Space  Command  is  the  primary  driver  for  i  mi  t  ovement  s 
in  the  capability  to  r  p.c-e  i  f  \  ind  predict  the  atmerpheri  r  neutral 


density  and  for  the  development  of  an  improved  atmospheric  neutral 
density  model.  The  Air  Force  Space  Command  Directorate  of 
Astrodynamics  initially  stated  the  command's  requirements  in  a 
letter  to  the  ^th  Weather  Wing  Aerospace  Sciences  Division  on  21 
May  1986.  On  25  May  1988,  the  Air  Force  Space  Command  Deputy  Chief 
of  Staff  for  Plans  elevated  and  reinforced  the  command's 
requirements  in  correspondence  to  his  counterparts  at  Military 
Airlift  Command  and  Air  Force  Systems  Command;  he  forwarded  a 
letter  from  the  United  States  Space  Command  Director  of  Operations 
that  stated  refined  requirements  for  improved  density  models.  Such 
specific  articulated  requirements  motivated  Air  Weather  Service  to 
develop  an  improved  capability  to  specify  and  predict  the 
atmospheric  neutral  density  at  satellite  altitudes. 

The  Air  Force  Geophysics  laboratory  is  about  to  award  a 
contract  for  advanced  development  of  the  Vector  Spherical  Harmonics 
model,  a  derivative  of  the  highly  sophisticated  Thermospheric  - 
Ionospheric  General  Circulation  Model  that  was  developed  at  the 
National  Center  for  Atmospheric  Research.  Military  Airlift  Command 
will  award  a  contract  for  operational  software  development  soon 
after  the  advanced  development  effort  begins.  This  operational 
software  development  will  ensure  that  the  Space  Forecast  Center  can 
execute  the  Vector  Spherical  Harmonics  model  efficiently  and  that 
it  will  generate  products  that  customers,  such  as  Air  Force  Space 
Command,  can  use  effectively. 

At  this  Atmospheric  Neutral  Density  Specialist  Conference, 
representatives  from  Air  Weather  Service  and  Air  Force  Space 
Command  agreed  that  they  should  begin  a  dialogue  to  work  out  the 
interface  between  the  new  model's  output  and  the  customer's  input. 
This  dialogue  has  begun  and  must  continue  throughout  operational 
software  development  to  ensure  that  the  new  model  offers  a 
quantifiable  improvement  over  the  older  Jacchia  models  and  that  it 
is  computationally  efficient. 

Other  related  studies  continue  as  the  new  model  advanced 
development  and  operational  software  development  proceed.  For 
example,  Air  Weather  Service  continues  to  search  for  an  improved 
capability  to  predict  the  solar  radio  flux,  the  Ottawa  Flux 
( F i o ) >  and  t>e  geomagnetic  index,  Ap,  that  remain  the  primary 
inputs  to  currently  operational  atmospheric  neutral  density  models. 
Air  Weather  Service  also  continues  to  monitor  the  continued 
availability  of  the  Ottawa  Flux  in  light  of  repeated  threats  by  the 
Canadian  National  Research  Council  to  close  down  the  Algonquin 
Observatory  whore  the  world  standard  Fip  is  measured  each  day. 

The  conference  participants  seemed  pleased  with  the 
presentations,  discussion,  and  progress  made  during  this 
conference.  Tn  fact ,  during  the  closing  session  it  was  suggested 
that  a  meeting  or  a  conference  such  as  this  one  should  he  convened 
on  a  periodic  and  recurring  basis,  perhaps  annually.  One  specific 
goal  for  a  future  meeting  of  this  nature  could  ho  to  investigate 
the  availability  and  the  possiMity  of  sharing  databases. 
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